Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 554,004 A→G 47.4% M151T (ATG→ACG)  ybiY ← predicted pyruvate formate lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE554,0040AG47.4% ‑2.3 / 31.0 38M151T (ATG→ACG) ybiYpredicted pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base A (13/7);  new base G (8/10);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

AGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATC  >  minE/553951‑554068
                                                     |                                                                
aGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCgg                                                     >  1:1570467/1‑67 (MQ=255)
aGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGc                                                    >  1:1040649/1‑68 (MQ=255)
aGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGc                                                    >  1:517335/1‑68 (MQ=255)
aGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGc                                                    >  1:692743/1‑68 (MQ=255)
   tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                                  <  1:116816/67‑1 (MQ=255)
   tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                                  <  1:2799231/67‑1 (MQ=255)
   tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                                  <  1:2084235/67‑1 (MQ=255)
   tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                                  <  1:652622/67‑1 (MQ=255)
   tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                                  <  1:1557900/67‑1 (MQ=255)
   tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                                  <  1:803029/67‑1 (MQ=255)
   tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                                  <  1:1236819/67‑1 (MQ=255)
                 tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGc                                        >  1:317294/1‑63 (MQ=255)
                 tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                  >  1:1912952/1‑69 (MQ=255)
                 tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                  >  1:2412750/1‑69 (MQ=255)
                 tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                  >  1:328482/1‑69 (MQ=255)
                 tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                                  >  1:2774543/1‑69 (MQ=255)
                            ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                       <  1:647061/69‑1 (MQ=255)
                            ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                       <  1:390170/69‑1 (MQ=255)
                            ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                       <  1:3314223/69‑1 (MQ=255)
                            ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                       <  1:915567/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:2895660/69‑1 (MQ=255)
                            ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                       <  1:1625375/69‑1 (MQ=255)
                             tCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt                       <  1:1694716/68‑1 (MQ=255)
                                     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)
                                     ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc              >  1:3384847/1‑69 (MQ=25)
                                     ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc              >  1:3410239/1‑69 (MQ=25)
                                     ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc              >  1:1342487/1‑69 (MQ=25)
                                     ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc              >  1:499818/1‑69 (MQ=25)
                                                gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt   >  1:1978631/1‑69 (MQ=255)
                                                gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt   >  1:1336350/1‑69 (MQ=255)
                                                gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt   >  1:541103/1‑69 (MQ=255)
                                                gcgcCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt   >  1:1921113/1‑69 (MQ=255)
                                                  gccatcgccatTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGAt   >  1:2970276/1‑67 (MQ=255)
                                                  gccatcgccatTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATc  >  1:1775444/1‑68 (MQ=255)
                                                  gccatcgccatTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATc  >  1:3022680/1‑68 (MQ=255)
                                                     |                                                                
AGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCCGCCGCTGCGATC  >  minE/553951‑554068

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: