Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,408,853 G→A 56.8% S374S (TCG→TCA yejF → fused predicted oligopeptide transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,408,8530GA56.8% 16.3 / 57.3 74S374S (TCG→TCAyejFfused predicted oligopeptide transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (25/7);  new base A (20/22);  total (45/29)
Fisher's exact test for biased strand distribution p-value = 9.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

ACCTATTCGTCATCGCATTCAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGAAGGCTTACGGGTTCACCAGCCGACGCTT  >  minE/1408805‑1408919
                                                |                                                                  
aCCTATTCGTCATCGCATTCAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGt                                                <  1:75821/69‑1 (MQ=255)
       cGTCATCGCATTCAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAg                                         >  1:825932/1‑69 (MQ=255)
       cGTCATCGCATTCAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAg                                         >  1:1273598/1‑69 (MQ=255)
       cGTCATCGCATTCAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAg                                         >  1:1739434/1‑69 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAgg                               >  1:587727/1‑67 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGaa                             >  1:1427221/1‑69 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGaa                             >  1:2044217/1‑69 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGaa                             >  1:2300406/1‑69 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGaa                             >  1:2348701/1‑69 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGaa                             >  1:1832839/1‑69 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGaa                             >  1:2508829/1‑69 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGaa                             >  1:2772031/1‑69 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGaa                             >  1:3170432/1‑69 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGaa                             >  1:1410465/1‑69 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGaa                             >  1:3257367/1‑69 (MQ=255)
                   cAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGaa                             >  1:965942/1‑69 (MQ=255)
                       tggtATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGAAGGCt                         <  1:883064/69‑1 (MQ=255)
                        ggtATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGatta                                     >  1:1334037/1‑56 (MQ=255)
                        ggtATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAgg                               >  1:3223536/1‑62 (MQ=255)
                        ggtATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAgg                               >  1:1062076/1‑62 (MQ=255)
                        ggtATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAgg                               >  1:3092460/1‑62 (MQ=255)
                        ggtATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAgg                               >  1:3273285/1‑62 (MQ=255)
                        ggtATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAgg                               >  1:2474975/1‑62 (MQ=255)
                        ggtATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAgg                               >  1:3375834/1‑62 (MQ=255)
                        ggtATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAgg                               >  1:1946895/1‑62 (MQ=255)
                              tCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGAAGGCTTACggg                   >  1:1405789/1‑68 (MQ=255)
                              tCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGAAGGCTTACggg                   >  1:2759525/1‑68 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGc                                           >  1:879497/1‑40 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:1619824/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:3348195/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:544618/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:1437909/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:3118334/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:653661/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:1184858/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:2822653/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:789781/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:1666785/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:1986311/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:2125/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:215192/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:2307618/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:2400483/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:2489099/1‑41 (MQ=255)
                                  gATCCAAACTCCTCACTCAACCCACGACTCAACGTTTTGCa                                          >  1:1437325/1‑41 (MQ=255)
                                  gATCCAAACTCCTCAATCAACCCACGACTCAACGTTTTGCaa                                         >  1:2965504/1‑41 (MQ=38)
                                  gATCCAAACTCCTAACTCAACCCACGACTCAACGTTTTGCaa                                         >  1:2861282/1‑41 (MQ=255)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:692188/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:831908/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:599503/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:111358/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:1218723/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:1003120/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:562449/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:2026477/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:1318819/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:1136403/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:1399196/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:3223607/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:3211707/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:3125032/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:1909784/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:1512883/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:3016709/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:3331712/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:2842717/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:2590679/37‑2 (MQ=39)
                                       aaaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:2556835/37‑2 (MQ=39)
                                        aaCTCCTCACTCAACCCACGACTCAACGTTTTGCaa                                         <  1:2849399/36‑2 (MQ=39)
                                               cGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGAAGGCTTACGGGTTCACCAGCCGACGCtt  <  1:822693/68‑1 (MQ=255)
                                               cGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGAAGGCTTACGGGTTCACCAGCCGACGCtt  <  1:1860072/68‑1 (MQ=255)
                                               cGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGAAGGCTTACGGGTTCACCAGCCGACGCtt  <  1:592083/68‑1 (MQ=255)
                                               cGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGAAGGCTTACGGGTTCACCAGCCGACGCtt  <  1:2199124/68‑1 (MQ=255)
                                               cGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGAAGGCTTACGGGTTCACCAGCCGACGCtt  <  1:2070539/68‑1 (MQ=255)
                                                |                                                                  
ACCTATTCGTCATCGCATTCAGGTGGTATTTCAGGATCCAAACTCCTCGCTCAACCCACGACTCAACGTTTTGCAGATTATTGAGGAAGGCTTACGGGTTCACCAGCCGACGCTT  >  minE/1408805‑1408919

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: