Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 899,197 | C→T | 41.2% | P78L (CCT→CTT) | pspC → | transcriptional activator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 899,197 | 0 | C | T | 41.2% | 20.0 / 36.7 | 34 | P78L (CCT→CTT) | pspC | transcriptional activator |
Reads supporting (aligned to +/- strand): ref base C (16/4); new base T (5/9); total (21/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.37e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01 |
TGCTTACATCATTTTGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACCTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACTGGCGGCAAGTGAAACGCGTTTACGCGAG > minE/899129‑899264 | tGCTTACATCATTTTGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACCtt > 1:642532/1‑71 (MQ=255) gCTTACATCATTTTGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACCTTc > 1:914755/1‑71 (MQ=255) aCATCATTTTGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACTTTc > 1:1285337/1‑67 (MQ=255) tcatTTTGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACTTTCCagcagc > 1:1490146/1‑71 (MQ=255) catTTTGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACCTTCCagcagc > 1:1181753/1‑70 (MQ=255) catTTTGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACCTTCCagcagc > 1:346328/1‑70 (MQ=255) catTTTGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGAAGCTACCTTCCagcagc > 1:931608/1‑70 (MQ=255) atTTTGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACTTTCCAGCAGCg < 1:1827590/70‑1 (MQ=255) atTTTGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACCTTCCAGCAGCg < 1:309865/70‑1 (MQ=255) ttttGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACCTTCCAGCAGCGaa > 1:789528/1‑71 (MQ=255) tGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACCTTCCAGCAGCGAATTGCTGGAt > 1:2112669/1‑70 (MQ=255) tGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGATACCTTCCAGCAGCGAATTGCTGGa > 1:591593/1‑69 (MQ=255) cgcTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCACCTACCTTCCAGCAGCGAATTGCTGTATGa < 1:480222/70‑1 (MQ=255) gcTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACTTTCCAGCAGCg < 1:1153800/56‑1 (MQ=255) gcTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACCTTCCAGCAGCGAATTGCTGGATGAAg > 1:1477845/1‑71 (MQ=255) gATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACTTTCCAGCAGGGAATTGCTGGATGAAGTCg < 1:1434519/70‑1 (MQ=255) gATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACTTTCCAGCAGCGAATTGCTGGATGAAGTCg < 1:1628052/70‑1 (MQ=255) aTGCCGGACAACATGGCCTTTGGTGAGCAGCTACTTTCCAGCAGCGAAtt < 1:154486/50‑1 (MQ=255) aTGCCGAACAACATGGCCTTTGGTGAGCAGCTACTTTCCAGCAGCGAAtt < 1:354783/50‑1 (MQ=255) tGCCGGACAACATGGCCTTTGGTGAGCAGCTACTTTCCAGCAGCGAAt < 1:517182/48‑1 (MQ=255) tGCCGGACAACATGGCCTTTGGTGAGCAGCTACTTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGaa < 1:2066681/71‑1 (MQ=255) gCCGGACAACATGGCCTTTGGTGAGCAGCTACCTTCCAGCAGCGaa > 1:1685398/1‑46 (MQ=255) ccGGACAACATGGCCTTTGGTGAGCAGCTACTTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACt > 1:1176133/1‑71 (MQ=255) ccGGACAACATGGCCTTTGGTGAGCAGCTACCTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACt > 1:2163712/1‑71 (MQ=255) aacaTGGCCTTTGGTGAGCAGCTACTTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACTggcgg > 1:166008/1‑70 (MQ=255) aacaTGGCCTTTGGTGAGCAGCTACTTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACTggcgg > 1:163379/1‑70 (MQ=255) aacaTGGCCTTTGGTGAGCAGCTACCTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACTggcggc < 1:1956286/71‑1 (MQ=255) aTGGCCTTTGGTGAGCAGCTACCTTCCAGCAGCGAATTGCTGGATGAAGTCGa > 1:2131689/1‑53 (MQ=255) aTGGCCTTTGGTGAGCAGCTACCTTCCAGCAGCGAATTGCTGGATGAAGTCGa > 1:1722515/1‑53 (MQ=255) tGGTGAGCAGCTACCTTCCAGCAGCGAATTGCTGTATGAAGTCGACCGTGAACTGGCGGCAAGTGaaa > 1:862956/1‑68 (MQ=255) tGGTGAGCAGCTACCTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACTGGCGGCAAGTGaaa > 1:1776804/1‑68 (MQ=255) cagcTACCTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACTGGCGGCAAGTGAAACGCGTTTAcg > 1:355817/1‑71 (MQ=255) cTACTTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACTGGCGGCAAGTGAAACGCGTTTAc < 1:1719054/67‑1 (MQ=255) cTACCTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACTGGCGGCAAGTGAAACGCGTTTACGCga < 1:26903/71‑1 (MQ=255) ccTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACTGGCGGCAAGTGAAACGCGTTTACGCgag < 1:1867517/69‑1 (MQ=255) cTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTg < 1:1756617/36‑1 (MQ=255) | TGCTTACATCATTTTGTCATTTGCGCTTGATCCAATGCCGGACAACATGGCCTTTGGTGAGCAGCTACCTTCCAGCAGCGAATTGCTGGATGAAGTCGACCGTGAACTGGCGGCAAGTGAAACGCGTTTACGCGAG > minE/899129‑899264 |
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: ‑ |