Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 554,036 A→G 62.2% G140G (GGT→GGC ybiY ← predicted pyruvate formate lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE554,0360AG62.2% 12.1 / 14.1 29G140G (GGT→GGCybiYpredicted pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base A (7/4);  new base G (8/10);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATG  >  minE/553968‑554102
                                                                    |                                                                  
tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                                                    >  1:1912952/1‑69 (MQ=255)
tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                                                    >  1:2774543/1‑69 (MQ=255)
tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                                                    >  1:2412750/1‑69 (MQ=255)
tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                                                    >  1:328482/1‑69 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                         <  1:390170/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                         <  1:1625375/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                         <  1:3314223/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                         <  1:3118498/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                         <  1:3078550/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                         <  1:3054992/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                         <  1:647061/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                         <  1:2895660/69‑1 (MQ=255)
           ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                         <  1:915567/69‑1 (MQ=255)
            tCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                                                         <  1:1694716/68‑1 (MQ=255)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                                >  1:3410239/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                                >  1:3384847/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                                >  1:2664145/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                                >  1:499818/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                                >  1:289767/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                                >  1:1342487/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                                >  1:2590993/1‑69 (MQ=25)
                    ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc                                                >  1:2630402/1‑69 (MQ=25)
                               gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt                                     >  1:1336350/1‑69 (MQ=255)
                               gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt                                     >  1:1921113/1‑69 (MQ=255)
                               gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt                                     >  1:1978631/1‑69 (MQ=255)
                               gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt                                     >  1:541103/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:1664293/68‑1 (MQ=255)
                                                           gggCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACgg          <  1:3240785/68‑1 (MQ=255)
                                                           gggCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACgg          <  1:2729036/68‑1 (MQ=255)
                                                           gggCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACgg          <  1:683021/68‑1 (MQ=255)
                                                                  ccaccCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCAt   >  1:1438154/1‑68 (MQ=255)
                                                                   caccCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATg  <  1:2215333/68‑1 (MQ=255)
                                                                   caccCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATg  <  1:1444967/68‑1 (MQ=255)
                                                                   caccCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATg  <  1:2239754/68‑1 (MQ=255)
                                                                   caccCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATg  <  1:2252242/68‑1 (MQ=255)
                                                                   caccCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATg  <  1:2386954/68‑1 (MQ=255)
                                                                   caccCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATg  <  1:2823906/68‑1 (MQ=255)
                                                                   caccCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATg  <  1:2628291/68‑1 (MQ=255)
                                                                   caccCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATg  <  1:1330183/68‑1 (MQ=255)
                                                                   caccCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATg  <  1:2663002/68‑1 (MQ=255)
                                                                    |                                                                  
TGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATCGTAAAACGGTTTATCGCGCAGAACGGTCGTCATG  >  minE/553968‑554102

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: