| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,621,151 | G→A | 75.1% | L41L (CTG→TTG) | eutQ ← | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,621,151 | 0 | G | A | 75.1% | 28.4 / 6.2 | 28 | L41L (CTG→TTG) | eutQ | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (3/4); new base A (3/18); total (6/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.44e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGATCGATTCATCACATTCGGTGATGGTAAAACCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGACATTGCCTG > minE/1621115‑1621217 | ggATCGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGCCACtt > 1:742851/1‑51 (MQ=255)ggATCGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGCCACtt > 1:370004/1‑51 (MQ=255)ggATCGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGCCACtt > 1:2237102/1‑51 (MQ=255) cGATTCATCACATTCGGTGATGTTAAAACCCAACAGATCAGcc < 1:2268125/43‑1 (MQ=25) cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:2285105/43‑1 (MQ=37) cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:767062/43‑1 (MQ=37) cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:595412/43‑1 (MQ=37) cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:526061/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:1678558/43‑1 (MQ=37) cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:120401/43‑1 (MQ=37) cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:1439959/43‑1 (MQ=37) cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:1460292/43‑1 (MQ=37) cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:1570555/43‑1 (MQ=37) cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:1581952/43‑1 (MQ=37) cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:1603741/43‑1 (MQ=37) cGATTCATCACATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:2186255/43‑1 (MQ=37) cGATTCATCACATTCGGTGATGGTAAAACACAACAGATCAGcc < 1:2965043/43‑1 (MQ=25) tcatcaCATTCGGTGATGGTAAAACCCAACAGATCAGcc < 1:468410/39‑1 (MQ=37) aaaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGACAt > 1:3283682/1‑69 (MQ=255) aaaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGACAt > 1:115965/1‑69 (MQ=255) aaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGa < 1:3090407/65‑1 (MQ=255) aaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGa < 1:1554658/65‑1 (MQ=255) aaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGa < 1:3307376/65‑1 (MQ=255) aaaCCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGa < 1:1413214/65‑1 (MQ=255) cAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGACATTGCCTg > 1:1539188/1‑69 (MQ=255) | GGATCGATTCATCACATTCGGTGATGGTAAAACCCAGTAGATCCGCCACTTCACGGGCTTCCGGGGTAATAATGCTGGCGCGCAGAACCACGGACATTGCCTG > minE/1621115‑1621217 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |