Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,840,168 C→A 20.0% D105Y (GAC→TAC)  gyrB ← DNA gyrase subunit B

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,840,1680C→A20.0% 2.9 ‑3.0 0.039313gyrBDNA gyrase subunit B
Reads supporting (aligned to +/- strand):  new base (1/8):  ref base (16/14):  total (17/22)
Fisher's exact test strand distribution p-value = 5.18e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATC  >  REL606/3840133‑3840202
                                   |                                  
gTGCAGACCGCCGGACACTTTATAGGAGTTATCGTc                                    <  1:5658828/36‑1 (MQ=255)
   cAGACCGCCGGACACTTTATAGGAGTTATCGTCaaa                                 <  1:313561/36‑1 (MQ=255)
      aCCGCCGGACACTTTATAGGAGTTATCGTCAAAttt                              >  1:5556973/1‑36 (MQ=255)
       ccgccgGACACTTTATAGGAGTTATCGTCAAATTTa                             <  1:540838/36‑1 (MQ=255)
       ccgccgGACACTTTATAGGAGTTATCGTCAAATTTa                             <  1:4013055/36‑1 (MQ=255)
        cgccgGACACTTTATAGGAGTTATCGTCAAATTTAc                            <  1:4338495/36‑1 (MQ=255)
        cgccgGACACTTTATAGGAGTTATCGTCAAATTTAc                            <  1:4247195/36‑1 (MQ=255)
        cgccgGACACTTTATAGGAGTTATCGTCAAATTTAc                            >  1:3721418/1‑36 (MQ=255)
         gccgGACACTTTATAGGAGTTATCGTCAAATTTAcc                           <  1:4121250/36‑1 (MQ=255)
            ggACACTTTATAGGAGTTATCGTCAAATTTACCTcc                        >  1:1559764/1‑36 (MQ=255)
            ggACACTTTATAGGAGTTATCGTCAAATTTACCGcc                        <  1:961970/36‑1 (MQ=255)
            ggACACTTTATAGGAGTTATCGTCAAATTTACCGcc                        <  1:2313119/36‑1 (MQ=255)
             gACACTTTATAGGAGTTATCGTCAAATTTACCGCCt                       <  1:4372315/36‑1 (MQ=255)
             gACACTTTATAGGAGTTATCGTCAAATTTACCGCCt                       >  1:4193544/1‑36 (MQ=255)
              acacTTTATAGGAGTTATCGTCAAATTTACCGCCTg                      >  1:227145/1‑36 (MQ=255)
              acacTTTATAGGAGTTATCGTCAAATTTACCGCCTg                      >  1:3983908/1‑36 (MQ=255)
               ctttttGTAGGATTTACCGAAAATTTTACCGCCTGc                     <  1:3253093/33‑1 (MQ=2)
                 cTTTATGGGAGTTATCGTAAAATTTCCCGCCTGCGt                   <  1:5029345/36‑1 (MQ=25)
                  tttATAGGAGTTATCGTCAAATTTACCGCCTCCGTg                  >  1:2479412/1‑36 (MQ=255)
                   ttATAGGAGTTATCGTAAAATTTACCGCCTGCGTGc                 <  1:3863062/36‑1 (MQ=255)
                    tataGGAGTTATCGTCAAATTTACCGCCTGCGTGCa                >  1:533151/1‑36 (MQ=255)
                    tataGGAGTTATCGTCAAATTTACCGCCTGCGTGCa                >  1:475860/1‑36 (MQ=255)
                    tataGGAGTTATCGTAAAATTTACCGCCTGCGTGCa                <  1:5126228/36‑1 (MQ=255)
                     ataGGAGTTATCGTCAAATTTACCGCCTGCGTGCAg               >  1:2463214/1‑36 (MQ=255)
                     ataGGAGTTATAGTCAAATTTACCGCCTGAGTGCAg               >  1:4423027/1‑36 (MQ=37)
                      taGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGa              >  1:3462441/1‑36 (MQ=255)
                      taGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGa              >  1:3278822/1‑36 (MQ=255)
                        ggAGTTATCGTCAAATTTACCGCCTGCGTGCAGAAc            <  1:5045659/36‑1 (MQ=255)
                        ggAGTTATCGTAAAATTTACCGCCTGCGTGCAGaaa            >  1:4361465/1‑35 (MQ=255)
                           tttATCGAAAAATTTACCGCCTGCGTGCAGAACGGt         <  1:1318561/35‑1 (MQ=25)
                             tATCGTAAAATTTACCTCCTGCGTGCAGAACGGTCa       <  1:535038/36‑1 (MQ=37)
                             tATCGTAAAATTTACCGCCTGCGTGCAGAACGGTCa       <  1:806625/36‑1 (MQ=255)
                              aTCGTCAAATTTACCGCCTGCGTGCAGAACGGTCAt      <  1:1890229/36‑1 (MQ=255)
                               tCGTAAAATTTACCGCCTGCGTGCAGAACGGTCATg     <  1:5193847/36‑1 (MQ=255)
                                cttCAATTTTCCCGCCTGCGTGCAGAACGGTCATGa    <  1:1352280/34‑1 (MQ=37)
                                 gTCAAATTTACCGCCTGCGTGCAGAACGGTCATGAt   >  1:5336537/1‑36 (MQ=255)
                                 gTCAAATTTACCGCCTGCGTGCAGAACGGTCATGAt   <  1:5456511/36‑1 (MQ=255)
                                 gTCAAATTTACCGCCTGCGTGCAGAACGGTCATGAt   >  1:1565894/1‑36 (MQ=255)
                                 gTCAAATTTACCGCCTGCGTGCAGAACGGTCATGAt   <  1:5081677/36‑1 (MQ=255)
                                 gTCAAATTTACCGCCTGCATGCACAACGGGCTTGAt   >  1:2909955/1‑36 (MQ=16)
                                  acAAATTCACCGCCTGCGTTCAGAACGGTCATGATc  <  1:5588907/35‑1 (MQ=37)
                                   |                                  
GTGCAGACCGCCGGACACTTTATAGGAGTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATC  >  REL606/3840133‑3840202

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 14 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG