Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,493,676 T→C 19.0% intergenic (+58/‑317) gyrB → / → ibpA DNA gyrase, subunit B/heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,493,6760TC19.0% 192.4 / 33.1 84intergenic (+58/‑317)gyrB/ibpADNA gyrase, subunit B/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base T (27/41);  new base C (9/7);  total (36/48)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

TCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAG  >  Exported/2493608‑2493743
                                                                    |                                                                   
tCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCtt                                                                    <  1:1940002/70‑1 (MQ=255)
tCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCtt                                                                    >  1:1967521/1‑70 (MQ=255)
tCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTg                                                                   <  1:4865322/71‑1 (MQ=255)
 cGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCtt                                                                    >  1:1070304/1‑69 (MQ=255)
             gcgcTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                       <  1:4091137/70‑1 (MQ=255)
             gcgcTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                       <  1:880253/70‑1 (MQ=255)
             gcgcTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAAt                                                      <  1:227000/71‑1 (MQ=255)
              cgcTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAAtt                                                     <  1:4010695/71‑1 (MQ=255)
               gcTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAAtt                                                     <  1:4099130/70‑1 (MQ=255)
                   aCCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGt                                                             <  1:4542828/58‑1 (MQ=255)
                      aTGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAt                                             <  1:2327127/71‑1 (MQ=255)
                       tGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                  <  1:3570730/65‑1 (MQ=255)
                       tGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                  <  1:633648/65‑1 (MQ=255)
                          gAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                         <  1:4002687/71‑1 (MQ=255)
                          gAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                         <  1:4009331/71‑1 (MQ=255)
                            gCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTa                                                            >  1:2053843/1‑50 (MQ=255)
                            gCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTa                                                            >  1:2254195/1‑50 (MQ=255)
                               tGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                  <  1:5127138/57‑1 (MQ=255)
                               tGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                  <  1:3892662/57‑1 (MQ=255)
                                 ccTGATGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATg                                                              >  1:2928634/1‑43 (MQ=14)
                                 ccTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTg                                         <  1:1247188/64‑1 (MQ=255)
                                  cTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                       <  1:2609100/49‑1 (MQ=255)
                                  cTGATGCGCTACGCTTATCAGGACTACGAAAATCTTGCAATGTATTGaa                                                       <  1:2547853/49‑1 (MQ=255)
                                   tGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATg                                                              <  1:5612276/41‑1 (MQ=255)
                                   tGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                       <  1:4711895/48‑1 (MQ=255)
                                    gATGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                  <  1:3508759/36‑1 (MQ=255)
                                    gATGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                  <  1:3413032/36‑1 (MQ=255)
                                    gATGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                  <  1:1947356/36‑1 (MQ=255)
                                    gATGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                  <  1:5310180/36‑1 (MQ=255)
                                    gATGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                  <  1:5587305/36‑1 (MQ=255)
                                     aTGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGAt                               <  1:254560/70‑1 (MQ=14)
                                     aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                  >  1:880369/1‑35 (MQ=255)
                                     aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTg                                                         >  1:4193514/1‑44 (MQ=255)
                                     aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTg                                                         >  1:3670612/1‑44 (MQ=255)
                                     aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTg                                                         >  1:1971126/1‑44 (MQ=255)
                                     aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                               <  1:4306832/70‑1 (MQ=255)
                                      tGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGAtt                                            >  1:1173598/1‑56 (MQ=14)
                                      tGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGAtt                                            >  1:4838595/1‑56 (MQ=14)
                                      tGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                          >  1:3730010/1‑58 (MQ=255)
                                      tGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                          >  1:3610521/1‑58 (MQ=255)
                                       gcgcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATaa                             <  1:3365579/70‑1 (MQ=255)
                                        cgcTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTa                                                            >  1:2611704/1‑38 (MQ=17)
                                          cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                         >  1:3665076/1‑71 (MQ=255)
                                          cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                         <  1:2182818/71‑1 (MQ=255)
                                          cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                         <  1:2305677/71‑1 (MQ=255)
                                           tACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                       <  1:4754536/40‑1 (MQ=255)
                                           tACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGt                        >  1:3218808/1‑71 (MQ=255)
                                            ccGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGc                                    >  1:2676404/2‑58 (MQ=14)
                                             cGCTTATCAGGCCTACGAAAATTCTGCAATTTATTGAATTTGCACGATTTTGTa                                       <  1:983921/54‑1 (MQ=11)
                                             cGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCa                                                 >  1:3049601/1‑44 (MQ=11)
                                             cGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTa                                       <  1:1981738/54‑1 (MQ=14)
                                              gCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                       <  1:1402563/37‑1 (MQ=255)
                                              gCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                       <  1:698211/37‑1 (MQ=255)
                                              gCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                       <  1:5158473/37‑1 (MQ=255)
                                              gCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                       <  1:4811310/37‑1 (MQ=255)
                                                ttATCAGGCCTACGAAAATCTTGCAATGTATTGAAtt                                                     >  1:3932463/1‑37 (MQ=255)
                                                  aTCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAAcg                   <  1:3397947/69‑1 (MQ=14)
                                                    cAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTg                                         >  1:3333938/1‑45 (MQ=14)
                                                    cAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTg                                         >  1:5193231/1‑45 (MQ=14)
                                                     aGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAACGCCGCa              <  1:4731642/71‑1 (MQ=14)
                                                     aGGCCTACGAAAATCTTGCAATGTATTGAATTTTTATGATTTTGTAGGc                                    >  1:4527492/1‑49 (MQ=255)
                                                         cTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGtt                       <  1:3876547/58‑1 (MQ=11)
                                                         cTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAACGCCGCATcc           <  1:1039345/70‑1 (MQ=11)
                                                         cTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCCg          >  1:2040621/1‑71 (MQ=11)
                                                         cTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTa                                       >  1:3301491/1‑42 (MQ=255)
                                                         cTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc           <  1:5466208/70‑1 (MQ=255)
                                                          tACGAAAATCTTGCAATGTATTGAATTTTCATGAttt                                           <  1:679215/37‑1 (MQ=255)
                                                          tACGAAAATCTTGCAATGTATTGAATTTTCATGAttt                                           <  1:4736070/37‑1 (MQ=255)
                                                          tACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgccgc               <  1:4336483/65‑1 (MQ=255)
                                                          tACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAc                    <  1:3905775/60‑1 (MQ=255)
                                                            cGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                          >  1:1183731/1‑36 (MQ=255)
                                                            cGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                          >  1:87395/1‑36 (MQ=255)
                                                            cGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                          >  1:2569587/1‑36 (MQ=255)
                                                            cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgc                  >  1:4880723/1‑60 (MQ=255)
                                                            cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgc                  >  1:4238418/1‑60 (MQ=255)
                                                            cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCa              <  1:2333479/64‑1 (MQ=255)
                                                            cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc        >  1:4910070/1‑70 (MQ=255)
                                                            cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc        >  1:2170843/1‑70 (MQ=255)
                                                            cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc        >  1:3010509/1‑70 (MQ=255)
                                                            cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc        >  1:3104643/1‑70 (MQ=255)
                                                              aaaaTCTTGCAATGTATTGAATTTTCATGATTTTGTa                                       >  1:2625088/1‑37 (MQ=255)
                                                              aaaaTCTTGCAATGTATTGAATTTTCATGATTTTGTa                                       >  1:1563655/1‑37 (MQ=255)
                                                              aaaaTCTTGCAATGTATTGAATTTTCATGATTTTGTa                                       >  1:2374485/1‑37 (MQ=255)
                                                               aaaTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaaga    <  1:1621509/71‑1 (MQ=255)
                                                                 aTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                         <  1:3741076/48‑1 (MQ=255)
                                                                 aTCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCAcgcc                 >  1:4676899/1‑56 (MQ=255)
                                                                   cTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGtt                       >  1:2264722/1‑48 (MQ=255)
                                                                    ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaagaag  >  1:2946024/1‑68 (MQ=255)
                                                                    ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaagaag  >  1:4330642/1‑68 (MQ=255)
                                                                    ttGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaagaa   >  1:462593/1‑67 (MQ=255)
                                                                    |                                                                   
TCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGAAG  >  Exported/2493608‑2493743

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: