Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 552,054 C→T 53.0% S490S (TCG→TCA ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE552,0540CT53.0% 1.6 / 56.3 66S490S (TCG→TCAybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base C (9/22);  new base T (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 = 1.00e+00

TAGTAACGGATTTGCGTATCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAGCGCC  >  minE/551992‑552120
                                                              |                                                                  
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:1107556/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:2012645/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:1039616/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:2722008/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:3096165/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:3110946/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:3153484/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:382792/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:458059/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:1443442/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:713207/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:1132488/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:1176087/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:732849/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:753394/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:1251486/69‑1 (MQ=255)
tagtaACGGATTTGCGTATCCCATGCGTCCATTACTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:868397/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:1558840/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                            <  1:2852789/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                            <  1:564404/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                            <  1:2630095/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                            <  1:1986199/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                            <  1:574615/69‑1 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                            >  1:2579782/1‑69 (MQ=255)
                  tCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAgg                                            <  1:887608/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:2302426/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                              <  1:3219900/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                              <  1:1686848/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                              <  1:3324127/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                              <  1:336945/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                              <  1:2981978/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                              <  1:2014983/68‑1 (MQ=255)
                                 aCTTCATCGAAGTTGTTGAAGTTACCCGCTGACAAc                                                              <  1:464646/36‑1 (MQ=37)
                                  cTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCt                              <  1:1187773/67‑1 (MQ=255)
                                        cGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcat                      <  1:1858553/69‑1 (MQ=255)
                                        cGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcat                      <  1:1597308/69‑1 (MQ=255)
                                        cGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcat                      <  1:467102/69‑1 (MQ=255)
                                        cGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcat                      <  1:701708/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:3406655/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga                 >  1:185710/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga                 >  1:2484537/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga                 >  1:492587/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga                 >  1:1477340/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga                 >  1:1435472/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga                 >  1:2006787/1‑66 (MQ=255)
                                             ttgttgAAGTTACCCGCTGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGcatcacga                 >  1:2539899/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:1278404/1‑69 (MQ=255)
                                                        aCCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAg      >  1:1291757/1‑69 (MQ=255)
                                                        aCCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAg      >  1:2668658/1‑69 (MQ=255)
                                                             ccGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAGcgcc  >  1:2319051/1‑68 (MQ=255)
                                                              |                                                                  
TAGTAACGGATTTGCGTATCCCACGCGTCCATCACTTCATCGAAGTTGTTGAAGTTACCTGCCGACAACGCTTTTTCTTGTGGCAGGAACACTTTGCCGCTGGTGGCATCATGCCCGCCTTCCAGCGCC  >  minE/551992‑552120

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: