Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,093,109 +C 8.0% coding (2210/2220 nt) ygiQ ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,093,1091.→C 8.0% 2.5 ‑3.0 ‑0.125coding (2210/2220 nt)ygiQhypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/2):  ref base (19/4):  total (19/6)
Fisher's exact test strand distribution p-value = 5.00e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.40e-01

GCGTGTCAGTTGCCGGGGATATTTTTTAACGT.CCAA.CCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTC  >  REL606/3093074‑3093144
                                     |                                   
gCGTGTCAGTTGCCGGGGATATTTTTTAACGT.CCaa                                      >  1:1231043/1‑36 (MQ=255)
 cGTGTCAGTTGCCGGGGATATTTTTTAACGT.CCAA.c                                    >  1:5219736/1‑36 (MQ=255)
  gtgtCAGTTGCCGGGGATATTTTTTAACGT.CCAA.cc                                   >  1:7590492/1‑36 (MQ=255)
  gtgtCAGTTGCCGGGGATATTTTTTAACGT.CCAA.cc                                   >  1:5261529/1‑36 (MQ=255)
   tgtCAGTTGCCGGGGATATTTTTTAACGT.CCAA.ccg                                  >  1:3767568/1‑36 (MQ=255)
    gtCAGTTGCCGGGGATATTTTTTAACGT.CCAA.ccgc                                 >  1:5834574/1‑36 (MQ=255)
      cAGTTGCCGGGGATATTTTTTAACGT.CCAA.ccgccg                               >  1:7238745/1‑36 (MQ=255)
        gTTGCCGGGGATATTTTTTAACGT.CCAA.CCGCCGCt                             >  1:7357804/1‑36 (MQ=255)
         ttGCCGGGGATATTTTTTAACGT.CCAA.CCGCCGCtt                            >  1:353945/1‑36 (MQ=255)
              ggggATATTTTTTAACGC.ACAA.CCTCCGGTTTATGc                       >  1:6481390/1‑36 (MQ=16)
                 gATTTTTTTAAACCC.CCCC.CCGCCGCTTTAGGCCGt                    <  1:1044174/36‑1 (MQ=12)
                 gATATTTTTTAACGT.CCAA.CCGCCGCTTTAGGCCGt                    >  1:6715427/1‑36 (MQ=255)
                 gATATTTTTTAACGT.CCAA.CCGCCGCTTTAGGCCGt                    >  1:2318219/1‑36 (MQ=255)
                  tttttttttACCGC.CCAC.CCGCCGCTTTAGGCCGtt                   <  1:6788201/33‑1 (MQ=19)
                  atatTTTTTAACGT.CCAA.CCGCCGCTTTAGGCCGtt                   >  1:4157477/1‑36 (MQ=255)
                  atatTTTTTAACGT.CCAA.CAGCCGCTTTAGGCCGtt                   >  1:2710558/1‑36 (MQ=255)
                   tatTTTTTAACGT.CCAA.CCGCCGCTTTATGCCGttt                  >  1:5045821/1‑36 (MQ=255)
                   tatTTTTTAACGT.CCAA.CCGCCGCTTTAGGCCGttt                  >  1:193959/1‑36 (MQ=255)
                   tatTTTTTAACGT.CCAA.CCGCCGCTTTAGGCCGttt                  >  1:6263972/1‑36 (MQ=255)
                   atttttttTACCC.CCAACCCGCCGCTTTAGGCCGtt                   <  1:7869574/34‑1 (MQ=12)
                   atttttttAACCCCCCAC.CCGCCGCTTTAGGCCGtt                   <  1:6295539/34‑1 (MQ=12)
                      tttttAACGT.CCAA.CCGCCGCTTTAGGCCGTttctt               >  1:6833056/1‑36 (MQ=255)
                       ttttAACGT.CCAA.CCGCCGCTTTAGGCCGTttcttc              >  1:4881058/1‑36 (MQ=255)
                       ttttAACCC.CCAACCCGCCGCTTTAGGCCGTttctt               <  1:4003080/36‑1 (MQ=19)
                        tttAACGT.CCAA.CCGCCGCTTTAGGCCGTTTCTTCg             >  1:1874424/1‑36 (MQ=255)
                         ttAACGT.CCAA.CCGCCGCTTTAGGCCGCTTCTTCga            >  1:4693569/1‑35 (MQ=255)
                            aCGT.CCAA.CCGCCGCTTTAGGCCGTTTCTTCGCAcc         >  1:6864191/1‑36 (MQ=255)
                             cGT.CCAA.CCGCCGCTTTAGGCCGTTTCTTCGCACCa        >  1:5299553/1‑36 (MQ=255)
                              gT.CCAA.CCGCCGCTTTAGGCCGTTTCTTCGCACCAg       >  1:545064/1‑36 (MQ=255)
                              gT.CCAA.CCGCCGCTTTAGGCCGTTTCTTCGCACCAg       >  1:3860172/1‑36 (MQ=255)
                              cccccAC.CCGCCGCTTTAGGCCGTTTCTTCGCACCa        <  1:559664/33‑1 (MQ=37)
                                 ccAA.CCGCCGCTTTAGGCCGTTTCTTCGCACCAGCa     >  1:7009475/1‑36 (MQ=255)
                                   aa.CCGCCGCTTTAGGCCGTTTCTTCGCACCAGCAtt   >  1:6866922/1‑36 (MQ=255)
                                   aa.CCGCCGCTTTAGGCCGTTTCTTCGCACCAGCAtt   >  1:4318912/1‑36 (MQ=255)
                                    a.CCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTc  >  1:2469026/1‑36 (MQ=255)
                                     |                                   
GCGTGTCAGTTGCCGGGGATATTTTTTAACGT.CCAA.CCGCCGCTTTAGGCCGTTTCTTCGCACCAGCATTC  >  REL606/3093074‑3093144

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG