Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,130,566 G→A 25.2% intergenic (‑305/‑219) thrS ← / → arpB threonyl‑tRNA synthetase/ECK1718:JW5278:b1720; hypothetical protein, N‑ter fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,130,5660GA25.2% 92.2 / 24.8 75intergenic (‑305/‑219)thrS/arpBthreonyl‑tRNA synthetase/ECK1718:JW5278:b1720; hypothetical protein, N‑ter fragment
Reads supporting (aligned to +/- strand):  ref base G (43/13);  new base A (12/7);  total (55/20)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAGTTCGTTGCGTGATTTGGCGGTAGT  >  minE/1130503‑1130629
                                                               |                                                               
aGTACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAg                                                             >  1:500422/1‑68 (MQ=255)
aGTACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAg                                                             >  1:1781430/1‑68 (MQ=255)
aGTACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAg                                                             >  1:1861403/1‑68 (MQ=255)
aGTACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAg                                                             >  1:1468584/1‑68 (MQ=255)
aGTACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAg                                                             >  1:2110232/1‑68 (MQ=255)
aGTACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAg                                                             >  1:2324302/1‑68 (MQ=255)
aGTACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAg                                                             >  1:841083/1‑68 (MQ=255)
aGTACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAg                                                             >  1:845061/1‑68 (MQ=255)
aGTACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAg                                                             >  1:988057/1‑68 (MQ=255)
  tACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCgg                                                          >  1:697940/1‑69 (MQ=255)
  tACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCgg                                                          <  1:657392/69‑1 (MQ=255)
  tACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCgg                                                          <  1:1012977/69‑1 (MQ=255)
  tACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCgg                                                          >  1:420541/1‑69 (MQ=255)
  tACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCgg                                                          >  1:2743468/1‑69 (MQ=255)
  tACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCgg                                                          >  1:249556/1‑69 (MQ=255)
  tACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCgg                                                          <  1:906568/69‑1 (MQ=255)
  tACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCgg                                                          <  1:931654/69‑1 (MQ=255)
  tACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCgg                                                          >  1:994881/1‑69 (MQ=255)
  tACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCgg                                                          >  1:1088658/1‑69 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGATCGGATTTACAGCAGCATAAACTGaa                                    >  1:1114907/1‑67 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGaa                                    >  1:1012397/1‑67 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGaa                                    >  1:2413165/1‑67 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGaa                                    >  1:141753/1‑67 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGaa                                    >  1:3201131/1‑67 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGaa                                    >  1:3215493/1‑67 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGaa                                    >  1:3261624/1‑67 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGaa                                    >  1:2081962/1‑67 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGaa                                    >  1:430619/1‑67 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGaa                                    >  1:1521064/1‑67 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGaa                                    >  1:818348/1‑67 (MQ=255)
                          caacTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGa                                     >  1:3416245/1‑66 (MQ=255)
                               gATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCa                                >  1:1742684/1‑66 (MQ=255)
                               gATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAg                               >  1:1378977/1‑67 (MQ=255)
                               gATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAg                               >  1:379810/1‑67 (MQ=255)
                               gATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAg                               >  1:53356/1‑67 (MQ=255)
                               gATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAg                               >  1:1191181/1‑67 (MQ=255)
                               gATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAg                               >  1:3038494/1‑67 (MQ=255)
                               gATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAg                               >  1:2393296/1‑67 (MQ=255)
                               gATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAg                               >  1:1526996/1‑67 (MQ=255)
                               gATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAg                               >  1:1445635/1‑67 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:142670/1‑69 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:2339889/1‑69 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:1735598/1‑69 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:2977756/1‑69 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:2991461/1‑69 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:3001247/1‑69 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:3427490/1‑69 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:206607/1‑69 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:314389/1‑69 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:1629596/1‑69 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:1097067/1‑69 (MQ=255)
                                 tAACATCACAGAGTGTCTGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCaaa                           >  1:198474/1‑69 (MQ=255)
                                       cACTGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTg                                      <  1:89573/52‑1 (MQ=255)
                                       cACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTg                                      <  1:1968784/52‑1 (MQ=255)
                                       cACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTg                                      <  1:1185232/52‑1 (MQ=255)
                                       cACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTg                                      <  1:1868625/52‑1 (MQ=255)
                                       cACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTg                                      <  1:1271627/52‑1 (MQ=255)
                                       cACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTg                                      <  1:1581811/52‑1 (MQ=255)
                                       cACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTg                                      <  1:2410799/52‑1 (MQ=255)
                                        aCGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAGTTCGtt                    >  1:3396531/1‑69 (MQ=255)
                                        aCGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAGTTCGtt                    >  1:1397119/1‑69 (MQ=255)
                                        aCGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAGTTCGtt                    >  1:211800/1‑69 (MQ=255)
                                        aCGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAGTTCGtt                    >  1:1614490/1‑69 (MQ=255)
                                        aCGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAGTTCGtt                    >  1:1806005/1‑69 (MQ=255)
                                        aCGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAGTTCGtt                    >  1:295362/1‑69 (MQ=255)
                                             gCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAGTTCGTTGCGTg               <  1:65828/69‑1 (MQ=255)
                                                  tGAAAGGTTTACTACGAGCGGATTTACAGCAGCATTAACTGAATCCAGCAAAg                          <  1:1599079/53‑1 (MQ=255)
                                                  tGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAg                          <  1:688542/53‑1 (MQ=255)
                                                  tGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAg                          <  1:472002/53‑1 (MQ=255)
                                                  tGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAg                          <  1:3401713/53‑1 (MQ=255)
                                                  tGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAg                          <  1:1465459/53‑1 (MQ=255)
                                                  tGAAAGGTTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAg                          <  1:2734828/53‑1 (MQ=255)
                                                  tGAAAGGGTTACTACGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAg                          <  1:3057294/53‑1 (MQ=255)
                                                          ttACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAGTTCGTTGCGTGATTTGGCGgtagt  >  1:1371543/1‑69 (MQ=255)
                                                          ttACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAGAAAAGTTCGTTGCGTGATTTGGCGgtagt  <  1:358683/69‑1 (MQ=255)
                                                               |                                                               
AGTACACCAGCAAAACGCTTAAAGAACAACTGATAACATCACGGTGCGTCTGAAAGGTTTACTGCGAGCGGATTTACAGCAGCATAAACTGAATCCAGCAAAGTTCGTTGCGTGATTTGGCGGTAGT  >  minE/1130503‑1130629

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