Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,628,879 | C→T | 33.2% | T109T (ACG→ACA) | yhhK ← | conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,628,879 | 0 | C | T | 33.2% | 31.3 / 19.0 | 42 | T109T (ACG→ACA) | yhhK | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (20/8); new base T (7/7); total (27/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.93e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTGACGATTAACACTTCTCCCAGCCGCCCTGTTGTGCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGATT > minE/2628815‑2628941 | ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc > 1:2475894/1‑68 (MQ=255) ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc > 1:2377731/1‑68 (MQ=255) ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc > 1:2279889/1‑68 (MQ=255) ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc > 1:2864037/1‑68 (MQ=255) ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc > 1:803898/1‑68 (MQ=255) ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc > 1:958429/1‑68 (MQ=255) ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc > 1:1161108/1‑68 (MQ=255) gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATcctc < 1:1523523/67‑3 (MQ=255) gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc < 1:1422515/67‑1 (MQ=255) gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc < 1:1685311/67‑1 (MQ=255) gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc < 1:990470/67‑1 (MQ=255) gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc < 1:2993171/67‑1 (MQ=255) gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc < 1:870670/67‑1 (MQ=255) gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc < 1:898170/67‑1 (MQ=255) gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc < 1:3056215/67‑1 (MQ=255) aaaTCCCAGTGCTCGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca < 1:2727212/69‑1 (MQ=255) aaaTCCCAGTGCTCGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca < 1:2690202/69‑1 (MQ=255) aaaTCCCAGTGCTCGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca < 1:1895666/69‑1 (MQ=255) aaaTCCCAGTGCTCGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca < 1:3254474/69‑1 (MQ=255) aaaTCCCAGTGCTAGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca < 1:305723/69‑1 (MQ=255) aaaTCCCAGTGCTAGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca < 1:1030591/69‑1 (MQ=255) aaaTCCCAGTGCTAGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca < 1:1157856/69‑1 (MQ=255) gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa > 1:817649/1‑69 (MQ=255) gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa > 1:2834944/1‑69 (MQ=255) gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa > 1:2152505/1‑69 (MQ=255) gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa > 1:2916705/1‑69 (MQ=255) gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa > 1:1765769/1‑69 (MQ=255) gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa > 1:3174182/1‑69 (MQ=255) cATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAgg > 1:2747724/1‑69 (MQ=255) cATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAgg > 1:1853340/1‑69 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCa > 1:1396426/1‑43 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATc > 1:99189/1‑45 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATc > 1:246366/1‑45 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATc > 1:2457516/1‑45 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATc > 1:2103674/1‑45 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATc > 1:1203590/1‑45 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt > 1:305355/1‑68 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt > 1:679245/1‑68 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt > 1:3001623/1‑68 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt > 1:908715/1‑68 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt > 1:1352616/1‑68 (MQ=255) aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt > 1:1051401/1‑68 (MQ=255) | TTGACGATTAACACTTCTCCCAGCCGCCCTGTTGTGCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGATT > minE/2628815‑2628941 |
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: ‑ |