Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,493,675 C→T 18.6% intergenic (+57/‑318) 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
*minE2,493,6750CT18.6% 204.8 / 31.6 86intergenic (+57/‑318)gyrB/ibpADNA gyrase, subunit B/heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base C (28/42);  new base T (9/7);  total (37/49)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

ATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGA  >  minE/2493605‑2493741
                                                                      |                                                                  
atatCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATc                                                                    >  1:3449625/1‑71 (MQ=255)
   tCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCtt                                                                  >  1:1967521/1‑70 (MQ=255)
   tCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCtt                                                                  <  1:1940002/70‑1 (MQ=255)
   tCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTg                                                                 <  1:4865322/71‑1 (MQ=255)
    cGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCtt                                                                  >  1:1070304/1‑69 (MQ=255)
                gcgcTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                     <  1:880253/70‑1 (MQ=255)
                gcgcTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGaa                                                     <  1:4091137/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:633648/65‑1 (MQ=255)
                          tGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                <  1:3570730/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:3892662/57‑1 (MQ=255)
                                  tGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTc                                                <  1:5127138/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:5310180/36‑1 (MQ=255)
                                       gATGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                <  1:5587305/36‑1 (MQ=255)
                                       gATGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                <  1:1947356/36‑1 (MQ=255)
                                       gATGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                <  1:3508759/36‑1 (MQ=255)
                                       gATGCGCTACGCTTATCAGGCCTACGAAAATCTTGc                                                                <  1:3413032/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:1971126/1‑44 (MQ=255)
                                        aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTg                                                       >  1:3670612/1‑44 (MQ=255)
                                        aTGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGAt                             <  1:4306832/70‑1 (MQ=255)
                                         tGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGAtt                                          >  1:4838595/1‑56 (MQ=14)
                                         tGCGCTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGAtt                                          >  1:1173598/1‑56 (MQ=14)
                                         tGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                        >  1:3610521/1‑58 (MQ=255)
                                         tGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGAtttt                                        >  1:3730010/1‑58 (MQ=255)
                                          gcgcTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATaa                           <  1:3365579/70‑1 (MQ=255)
                                           cgcTCCGCTTATCAGGCCTACGAAAATTCTGCAATGTa                                                          >  1:2611704/1‑38 (MQ=17)
                                             cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                       <  1:2182818/71‑1 (MQ=255)
                                             cTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCg                       >  1:3665076/1‑71 (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:5193231/1‑45 (MQ=14)
                                                       cAGGCCTACGAAAATTCTGCAATGTATTGAATTTGCACGATTTTg                                       >  1:3333938/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:4736070/37‑1 (MQ=255)
                                                             tACGAAAATCTTGCAATGTATTGAATTTTCATGAttt                                         <  1:679215/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:3010509/1‑70 (MQ=255)
                                                               cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc      >  1:3104643/1‑70 (MQ=255)
                                                               cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc      >  1:2170843/1‑70 (MQ=255)
                                                               cGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc      >  1:4910070/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)
                                                                      |                                                                  
ATATCGATATTTAATGGCGCTAACCATGCGAGCGTGCCTGATGCGCTACGCTTATCAGGCCTACGAAAATCTTGCAATGTATTGAATTTTCATGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGA  >  minE/2493605‑2493741

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: