Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,246,967 A→G 100% L130L (CTA→CTG znuC → zinc transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,246,9670AG96.7% 61.5 / ‑5.1 28L130L (CTA→CTGznuCzinc transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (16/11);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGT  >  minE/1246913‑1247005
                                                      |                                      
tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:647382/68‑1 (MQ=255)
tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:1684465/68‑1 (MQ=255)
tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:1946544/68‑1 (MQ=255)
tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:199684/68‑1 (MQ=255)
tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:400016/68‑1 (MQ=255)
tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:2062162/68‑1 (MQ=255)
tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:2066186/68‑1 (MQ=255)
tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:738896/68‑1 (MQ=255)
tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:3049573/68‑1 (MQ=255)
tcTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:2940353/68‑1 (MQ=255)
 cTGATTAACGCACCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCAtt                          <  1:676046/68‑1 (MQ=255)
            aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:1277028/1‑66 (MQ=255)
            aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:3243932/1‑66 (MQ=255)
            aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:3055903/1‑66 (MQ=255)
            aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:2502318/1‑66 (MQ=255)
            aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:894/1‑66 (MQ=255)
            aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:2300454/1‑66 (MQ=255)
            aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:2181258/1‑66 (MQ=255)
            aCCGATGCAAAAGCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCgac               >  1:1571005/1‑66 (MQ=255)
                  gCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCa           <  1:2303009/66‑1 (MQ=255)
                        gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAgcg                                  >  1:2008672/1‑37 (MQ=37)
                        gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:676705/1‑69 (MQ=255)
                        gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:995268/1‑69 (MQ=255)
                        gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:316605/1‑69 (MQ=255)
                        gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:3075790/1‑69 (MQ=255)
                        gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:2439770/1‑69 (MQ=255)
                        gCTCTCGGGTGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:1428572/1‑69 (MQ=255)
                        gCTCTCGGGAGGCGAAACGCAGCGTGTACTGTTAGCGCGAGCATTGTTAAATCGACCGCAATTATTAGt  >  1:700543/1‑69 (MQ=255)
                                                      |                                      
CCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGT  >  minE/1246913‑1247005

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: