Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 171,738 C→T 58.5% D14D (GAC→GAT frr → ribosome recycling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE171,7380CT58.5% 9.2 / 28.3 41D14D (GAC→GATfrrribosome recycling factor
Reads supporting (aligned to +/- strand):  ref base C (11/6);  new base T (4/20);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 2.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GAGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTGCTTCTCCCAGCCTG  >  minE/171671‑171804
                                                                   |                                                                  
gagaCAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACa                                                                   <  1:210131/69‑1 (MQ=255)
gagaCAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACa                                                                   <  1:1863868/69‑1 (MQ=255)
                                   gATATCAGAAACGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCaa                                 <  1:1197914/68‑1 (MQ=255)
                                   gATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCaa                                 <  1:2506154/68‑1 (MQ=255)
                                   gATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCaa                                 <  1:2273859/68‑1 (MQ=255)
                                   gATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCaa                                 <  1:1317677/68‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:3126476/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:2534295/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:852257/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:717212/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:3078923/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:2323261/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:3150629/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:3184955/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:3232862/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:3382195/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:519035/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:1646443/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:99721/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:1160067/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:143742/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:1465079/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:1491097/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:1720259/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAAt                              <  1:969057/69‑1 (MQ=255)
                                     tatCAGAAAAGATGCTGAAGTACGCATGGATAAATGCGTAGAAGAGTTCAAAACCCAAATCAGCAAAAt                              <  1:2250717/69‑1 (MQ=255)
                                                 tGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCTAATCAGCAAAATACGCACGGGTc                   >  1:8749/1‑68 (MQ=255)
                                                 tGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCg                  >  1:612122/1‑69 (MQ=255)
                                                 tGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCg                  >  1:3043977/1‑69 (MQ=255)
                                                 tGCTGAAGTACGCATGGATAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCg                  >  1:2268958/1‑69 (MQ=255)
                                                           cGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTGCTTCTccc        >  1:3066323/1‑69 (MQ=255)
                                                                 gACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGc                          >  1:3318195/1‑45 (MQ=255)
                                                                 gACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACggg                     >  1:3130582/1‑50 (MQ=255)
                                                                 gACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTGCTTCTCCCAGCctg  >  1:261560/1‑69 (MQ=255)
                                                                 gACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTGCTTCTCCCAGCctg  >  1:3388894/1‑69 (MQ=255)
                                                                 gACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTGCTTCTCCCAGCctg  >  1:3129795/1‑69 (MQ=255)
                                                                 gACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTGCTTCTCCCAGCctg  >  1:2897740/1‑69 (MQ=255)
                                                                 gACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTGCTTCTCCCAGCctg  >  1:1071866/1‑69 (MQ=255)
                                                                 gACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTGCTTCTCCCAGCctg  >  1:2276749/1‑69 (MQ=255)
                                                                 gACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTGCTTCTCCCAGCctg  >  1:2272940/1‑69 (MQ=255)
                                                                 gACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTGCTTCTCCCAGCctg  >  1:1075873/1‑69 (MQ=255)
                                                                   |                                                                  
GAGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTGCTTCTCCCAGCCTG  >  minE/171671‑171804

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: