Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,621,158 C→A 72.4% A38A (GCG→GCT eutQ ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,621,1580CA72.4% 29.8 / 14.2 29A38A (GCG→GCTeutQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (3/18);  total (7/22)
Fisher's exact test for biased strand distribution p-value = 6.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCGATTCATCACATTCGGTGATGGTAAAACCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGACATTGCCTGTTCG  >  minE/1621115‑1621221
                                           |                                                               
ggATCGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGCCACtt                                                          >  1:370004/1‑51 (MQ=255)
ggATCGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGCCACtt                                                          >  1:742851/1‑51 (MQ=255)
ggATCGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGCCACtt                                                          >  1:2237102/1‑51 (MQ=255)
    cGATTCATCACATTCGGTGATGTTAAAACCCAACAGATCAGcc                                                              <  1:2268125/43‑1 (MQ=25)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:2186255/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:3215529/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:3197070/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:265259/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:526061/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:2285105/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:1678558/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:1603741/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:1581952/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:1570555/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:595412/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:1460292/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:1439959/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:767062/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:120401/43‑1 (MQ=37)
    cGATTCATCACATTCGGTGATGGTAAAACACAACAGATCAGcc                                                              <  1:2965043/43‑1 (MQ=25)
        tcatcaCATTCGGTGATGGTAAAACCCAACAGATCAGcc                                                              <  1:468410/39‑1 (MQ=37)
                            aaaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGACAt            >  1:115965/1‑69 (MQ=255)
                            aaaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGACAt            >  1:3283682/1‑69 (MQ=255)
                             aaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGa               <  1:3307376/65‑1 (MQ=255)
                             aaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGa               <  1:3090407/65‑1 (MQ=255)
                             aaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGa               <  1:1554658/65‑1 (MQ=255)
                             aaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGa               <  1:1413214/65‑1 (MQ=255)
                                  cAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGACATTGCCTg      >  1:1539188/1‑69 (MQ=255)
                                      aGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGACATTGCCTGTTcg  >  1:840107/1‑69 (MQ=255)
                                          ccGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGAc              >  1:943683/1‑53 (MQ=255)
                                          ccGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGAc              >  1:2326351/1‑53 (MQ=255)
                                           |                                                               
GGATCGATTCATCACATTCGGTGATGGTAAAACCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGACATTGCCTGTTCG  >  minE/1621115‑1621221

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: