Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,080,000 C→G 31.1% K22N (AAG→AAC cbpA ← curved DNA‑binding protein, DnaJ homologue that functions as a co‑chaperone of DnaK

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,080,0000C→G31.1% 2.3 ‑4.5 0.02066cbpAcurved DNA‑binding protein, DnaJ homologue that functions as a co‑chaperone of DnaK
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (4/9):  total (11/9)
Fisher's exact test strand distribution p-value = 1.08e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGC  >  REL606/1079965‑1080034
                                   |                                  
ggATGGTATTTGCGGGCAAGTCGACGATAGGCGGTc                                    >  1:5871450/1‑36 (MQ=255)
  aTGGTATTTGCGGGCAAGTCGACGATAGGCGGTGtt                                  >  1:3321513/1‑36 (MQ=255)
     gTATTTGCGGGCAAGTCGACGATAGGCGGTTttttt                               >  1:3581141/1‑33 (MQ=37)
     gTATTTGCGGGCAAGTCGACGATAGGCGGGGTTGAt                               >  1:15222/1‑36 (MQ=37)
      tATTTGCGGGCAAGTCGACGATAGGCGGTCTTGAtt                              <  1:861393/36‑1 (MQ=255)
      tATTTGCGGGCAAGTCGACGATAGGCGGTCTTGAtt                              <  1:2913646/36‑1 (MQ=255)
        tttGCGGGCAAGTCGACGATAGGCGGTGTTGCtttt                            >  1:4256885/1‑34 (MQ=25)
           gCGGGCAAGTCGACGATCGGCGGTCTTGATTGTCtt                         <  1:4238418/36‑1 (MQ=255)
           gCGGGCAAGTCGACGATAGGCGGGGTTGGTTGTTtt                         >  1:3547137/1‑36 (MQ=19)
            cgcgCAAGTCGACGATAGGCGGTCTTGATTGTCTTg                        <  1:3602957/33‑1 (MQ=255)
            cGGGCAAGTCGACGATAGTCGGTCTTGATTGTCTTg                        <  1:2168356/36‑1 (MQ=255)
            cGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTg                        <  1:4736438/36‑1 (MQ=255)
              ggCAAGTCGACGATAGGCGGTCTTGATTGTCTTgag                      <  1:4204998/36‑1 (MQ=255)
               gCAAGTCGACGATAGGCGGTCTTGATTGTCTTgaga                     <  1:1912255/36‑1 (MQ=255)
                   gTCGACGATAAGCGGGGTTATTTTTCTTGAGATtgt                 >  1:1776614/1‑36 (MQ=2)
                    tCGACGATAGGCGGAGTTGATTGTTTTTTGATcgtg                >  1:1424925/1‑35 (MQ=12)
                      gacgaTAGGCGGTCTTGTTTTTCTTTAGATcgtccg              >  1:2305144/1‑36 (MQ=37)
                      gacgaTAGGCGGTCTTGATTTTCTTTTGATcgtccg              >  1:5039936/1‑36 (MQ=25)
                             ggCGGTCTTGATTGTCTTGAGATCGTCCGTCGGttt       <  1:4052650/36‑1 (MQ=255)
                                ggTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCAc    >  1:2359877/1‑36 (MQ=255)
                                  tCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGc  >  1:5863834/1‑36 (MQ=255)
                                   |                                  
GGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGC  >  REL606/1079965‑1080034

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG