| 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 | ||
| * | minE | 1,621,158 | 0 | C | A | 72.4% | 29.8 / 14.2 | 29 | A38A (GCG→GCT) | eutQ | conserved 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: ‑ |