Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 96,629 C→T 100% P237P (CCG→CCA yacF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE96,6290CT82.4% 67.4 / 6.4 45P237P (CCG→CCAyacFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (7/1);  new base T (37/0);  total (44/1)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTCACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCCGGTACCTGTCCGTTTTCAGTGTCCAGCGGCATAAAACGAATGGCAAAA  >  minE/96563‑96677
                                                                  |                                                
aTTCACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCTgg                                                >  1:2383596/1‑69 (MQ=255)
aTTCACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCTgg                                                >  1:96343/1‑69 (MQ=255)
aTTCACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCTgg                                                >  1:3050772/1‑69 (MQ=255)
aTTCACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCTgg                                                >  1:2709396/1‑69 (MQ=255)
aTTCACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCTgg                                                >  1:2556140/1‑69 (MQ=255)
aTTCACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCTgg                                                >  1:1990024/1‑69 (MQ=255)
aTTCACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCTgg                                                >  1:1834211/1‑69 (MQ=255)
aTTCACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCTgg                                                >  1:1197928/1‑69 (MQ=255)
   cACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCCGGTAc                                             >  1:661023/1‑69 (MQ=255)
   cACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCCGGTAc                                             <  1:2135804/69‑1 (MQ=255)
                tCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCCGGTACCTGTCCGTTTTc                                 >  1:1761715/1‑68 (MQ=255)
                tCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCCGGTACCTGTCCGTTTTCa                                >  1:70551/1‑69 (MQ=255)
                tCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCCGGTACCTGTCCGTTTTCa                                >  1:1823630/1‑69 (MQ=255)
                tCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCCGGTACCTGTCCGTTTTCa                                >  1:595234/1‑69 (MQ=255)
                tCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCCGGTACCTGTCCGTTTTCa                                >  1:1180688/1‑69 (MQ=255)
                tCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCCGGTACCTGTCCGTTTTCa                                >  1:1063366/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGGaaa   >  1:593225/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:409522/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:1495104/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:2811762/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:3017058/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:3352220/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:2566303/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:1060800/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:2477692/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:1575641/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:1627587/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:1290548/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:1278684/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:2114372/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaa   >  1:1994847/1‑69 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaa    >  1:1001274/1‑68 (MQ=255)
                                             cAGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaa    >  1:405150/1‑68 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTTGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaaa  >  1:2445112/1‑69 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaaa  >  1:3321179/1‑69 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaaa  >  1:3061779/1‑69 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaaa  >  1:3422043/1‑69 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaaa  >  1:2064336/1‑69 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaaa  >  1:2214083/1‑69 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaaa  >  1:1890583/1‑69 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaaa  >  1:2717120/1‑69 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaaa  >  1:1692995/1‑69 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaaa  >  1:2230702/1‑69 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGCGGCATAAAACGAATGGCaaaa  >  1:25706/1‑69 (MQ=255)
                                              aGTTCGAAATCCAGACGTTCTGGTACCTGTCCGTTTTCACTGTCCAGAGGCATAAAACGAATGGCaaaa  >  1:613682/1‑69 (MQ=255)
                                                                  |                                                
ATTCACCGTAATAGTTTCTGACATTTTTACTCCTTAGCAACAGGCCAGTTCGAAATCCAGACGTTCCGGTACCTGTCCGTTTTCAGTGTCCAGCGGCATAAAACGAATGGCAAAA  >  minE/96563‑96677

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: