Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 553,997 T→C 40.9% L153L (CTA→CTG ybiY ← predicted pyruvate formate lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE553,9970TC40.9% 12.8 / 29.1 44L153L (CTA→CTGybiYpredicted pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base T (9/17);  new base C (8/10);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCC  >  minE/553929‑554056
                                                                    |                                                           
tatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGt                                                             <  1:720823/69‑1 (MQ=255)
tatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGt                                                             <  1:2978445/69‑1 (MQ=255)
tatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGt                                                             <  1:2665644/69‑1 (MQ=255)
tatTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGt                                                             <  1:2179227/69‑1 (MQ=255)
      cACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgcgc                                                        <  1:944754/68‑1 (MQ=255)
      cACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgcgc                                                        <  1:796589/68‑1 (MQ=255)
      cACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgcgc                                                        <  1:724942/68‑1 (MQ=255)
      cACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgcgc                                                        <  1:2456371/68‑1 (MQ=255)
      cACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgcgc                                                        <  1:2339821/68‑1 (MQ=255)
       aCGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgcgc                                                        <  1:2182003/67‑1 (MQ=255)
       aCGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgcgc                                                        <  1:682752/67‑1 (MQ=255)
       aCGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgcgc                                                        <  1:311202/67‑1 (MQ=255)
       aCGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgcgc                                                        <  1:246454/67‑1 (MQ=255)
       aCGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAgcgc                                                        <  1:3099543/67‑1 (MQ=255)
                      aGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCgg                                         >  1:1570467/1‑67 (MQ=255)
                      aGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGc                                        >  1:1040649/1‑68 (MQ=255)
                      aGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGc                                        >  1:692743/1‑68 (MQ=255)
                      aGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGc                                        >  1:517335/1‑68 (MQ=255)
                         tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                      <  1:2084235/67‑1 (MQ=255)
                         tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                      <  1:652622/67‑1 (MQ=255)
                         tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                      <  1:2799231/67‑1 (MQ=255)
                         tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                      <  1:1236819/67‑1 (MQ=255)
                         tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                      <  1:116816/67‑1 (MQ=255)
                         tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                      <  1:1557900/67‑1 (MQ=255)
                         tttCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTg                                      <  1:803029/67‑1 (MQ=255)
                                       tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGc                            >  1:317294/1‑63 (MQ=255)
                                       tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                      >  1:328482/1‑69 (MQ=255)
                                       tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                      >  1:2412750/1‑69 (MQ=255)
                                       tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                      >  1:2774543/1‑69 (MQ=255)
                                       tGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTcacca                      >  1:1912952/1‑69 (MQ=255)
                                                  ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt           <  1:3078550/69‑1 (MQ=255)
                                                  ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt           <  1:3054992/69‑1 (MQ=255)
                                                  ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt           <  1:1625375/69‑1 (MQ=255)
                                                  ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt           <  1:3118498/69‑1 (MQ=255)
                                                  ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt           <  1:3314223/69‑1 (MQ=255)
                                                  ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt           <  1:915567/69‑1 (MQ=255)
                                                  ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt           <  1:2895660/69‑1 (MQ=255)
                                                  ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt           <  1:390170/69‑1 (MQ=255)
                                                  ctCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt           <  1:647061/69‑1 (MQ=255)
                                                   tCGTGGCTGGCTTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGt           <  1:1694716/68‑1 (MQ=255)
                                                           ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc  >  1:3384847/1‑69 (MQ=25)
                                                           ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc  >  1:1342487/1‑69 (MQ=25)
                                                           ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc  >  1:499818/1‑69 (MQ=25)
                                                           ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc  >  1:3410239/1‑69 (MQ=25)
                                                           ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc  >  1:2664145/1‑69 (MQ=25)
                                                           ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc  >  1:2590993/1‑69 (MQ=25)
                                                           ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc  >  1:289767/1‑69 (MQ=25)
                                                           ggctTGCAGCAGCGCCGTCGCCATTTCCGGCTGCATAAACGGTTCACCGCCAGAAAGCGTTAAgccgcc  >  1:2630402/1‑69 (MQ=25)
                                                                    |                                                           
TATTTCCACGGCACATGCAGACAGGTTTCTACCGCAGTATGAATGCCTGCCTCGTGGCTGGCTTGCAGTAGCGCCATCGCCATTTCCGGCTGCATAAAGGGCTCACCACCCGAAAGCGTTAAACCGCC  >  minE/553929‑554056

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: