Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 552,051 T→C 53.1% A491A (GCA→GCG ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE552,0510TC53.1% 6.3 / 54.0 66A491A (GCA→GCGybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base T (9/22);  new base C (15/20);  total (24/42)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TAGTAACGGATTTGCGTATCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAG  >  minE/551992‑552116
                                                           |                                                                 
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:1107556/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:2012645/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:458059/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:382792/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:713207/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:732849/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:2722008/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:753394/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:868397/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:1443442/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:3096165/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:1039616/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:1251486/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:3110946/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:1176087/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:1132488/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:3153484/69‑1 (MQ=255)
 agtaACGTATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:1567453/68‑1 (MQ=255)
 agtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:2004455/68‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCTACAACGCTTTTTCTTGTGGCAgg                                        <  1:1426253/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                        <  1:2852789/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                        <  1:2630095/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                        >  1:2579782/1‑69 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                        <  1:564404/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                        <  1:887608/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                        <  1:1558840/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                        <  1:574615/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                        <  1:1986199/69‑1 (MQ=255)
                             cATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGcc                             >  1:2484640/1‑69 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGTTTTTTCTTGTGGCAGGAACACTTTGCCGCt                          <  1:1116937/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCtg                         >  1:89362/1‑69 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCtg                         >  1:2946470/1‑69 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                          <  1:3324127/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                          <  1:2981978/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                          <  1:336945/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                          <  1:1686848/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                          <  1:2302426/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                          <  1:2014983/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                          <  1:3219900/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                          <  1:464646/36‑1 (MQ=37)
                                  cTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                          <  1:1187773/67‑1 (MQ=255)
                                        cGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcat                  <  1:1597308/69‑1 (MQ=255)
                                        cGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcat                  <  1:701708/69‑1 (MQ=255)
                                        cGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcat                  <  1:1858553/69‑1 (MQ=255)
                                        cGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcat                  <  1:467102/69‑1 (MQ=255)
                                         gAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcat                  <  1:3144609/68‑1 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACt                                  >  1:923717/1‑48 (MQ=39)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACt                                  >  1:2005248/1‑48 (MQ=39)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga             >  1:3110072/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga             >  1:1435472/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga             >  1:3406655/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga             >  1:2539899/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga             >  1:2484537/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga             >  1:2006787/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga             >  1:492587/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga             >  1:1477340/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga             >  1:185710/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacg              >  1:1687779/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacg              >  1:158013/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacg              >  1:1512812/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCACGc             >  1:1637841/1‑69 (MQ=255)
                                                        aCCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCa   >  1:1900352/1‑68 (MQ=255)
                                                        aCCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAg  >  1:326349/1‑69 (MQ=255)
                                                        aCCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAg  >  1:2668658/1‑69 (MQ=255)
                                                        aCCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAg  >  1:1291757/1‑69 (MQ=255)
                                                        aCCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAg  >  1:1278404/1‑69 (MQ=255)
                                                           |                                                                 
TAGTAACGGATTTGCGTATCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAG  >  minE/551992‑552116

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: