Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 554,030 C→T 54.5% E142E (GAG→GAA ybiY ← predicted pyruvate formate lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE554,0300CT54.5% ‑1.6 / 24.5 33E142E (GAG→GAAybiYpredicted pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base C (11/4);  new base T (8/10);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

TGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGG  >  minE/553968‑554094
                                                              |                                                                
tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGc                                                                  >  1:317294/1‑63 (MQ=255)
tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                                            >  1:2774543/1‑69 (MQ=255)
tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                                            >  1:1912952/1‑69 (MQ=255)
tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                                            >  1:328482/1‑69 (MQ=255)
tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                                            >  1:2412750/1‑69 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                 <  1:3054992/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                 <  1:3078550/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                 <  1:3118498/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                 <  1:2895660/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                 <  1:3314223/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                 <  1:390170/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                 <  1:647061/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                 <  1:1625375/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                 <  1:915567/69‑1 (MQ=255)
            tCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                 <  1:1694716/68‑1 (MQ=255)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                        >  1:3384847/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                        >  1:3410239/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                        >  1:499818/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                        >  1:1342487/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                        >  1:289767/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                        >  1:2664145/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                        >  1:2630402/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                        >  1:2590993/1‑69 (MQ=25)
                               gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt                             >  1:1978631/1‑69 (MQ=255)
                               gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt                             >  1:1921113/1‑69 (MQ=255)
                               gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt                             >  1:541103/1‑69 (MQ=255)
                               gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt                             >  1:1336350/1‑69 (MQ=255)
                                 gccatcgccatTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt                             >  1:2970276/1‑67 (MQ=255)
                                 gccatcgccatTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATc                            >  1:3022680/1‑68 (MQ=255)
                                 gccatcgccatTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATc                            >  1:1775444/1‑68 (MQ=255)
                                                           gggCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACgg  <  1:2729036/68‑1 (MQ=255)
                                                           gggCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACgg  <  1:3240785/68‑1 (MQ=255)
                                                           gggCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACgg  <  1:1664293/68‑1 (MQ=255)
                                                           gggCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACgg  <  1:683021/68‑1 (MQ=255)
                                                              |                                                                
TGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGG  >  minE/553968‑554094

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: