| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,952,826 | T→C | 19.1% | E218G (GAA→GGA) | rho ← | transcription termination factor Rho |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,952,826 | 0 | T | C | 19.1% | 36.7 / 5.9 | 21 | E218G (GAA→GGA) | rho | transcription termination factor Rho |
| Reads supporting (aligned to +/- strand): ref base T (6/11); new base C (2/2); total (8/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.18e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01 | |||||||||||
GGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTGGCACTCAGGGTTGTTACGGGTGATGTTGGCCGCAATGTTCTGCAGCATGATCGTCTTACCGGCTTTCGGCGGGGCAACGAT > NC_002947/5952692‑5952953 | ggATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTc < 2:28060/136‑1 (MQ=255) tAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGt < 2:365465/136‑1 (MQ=255) aCCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCa > 2:281961/1‑136 (MQ=255) tCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCtcgg > 1:273476/1‑135 (MQ=255) gcgcTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATcagc < 2:421633/136‑1 (MQ=255) cTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAgg > 2:410554/1‑136 (MQ=255) ttGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGa < 2:249808/136‑1 (MQ=255) cTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCTCTCTTCGATCAGCAGGACGATc < 2:309126/136‑1 (MQ=255) tctcGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAg < 2:390366/136‑1 (MQ=255) cGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTCCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTg > 2:326224/1‑136 (MQ=255) cGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTCCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTg > 2:328223/1‑136 (MQ=255) cGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTCCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTGGCACTCAGGgtt < 1:328223/136‑1 (MQ=255) cGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTCCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTGGCACTCAGGgtt < 1:326224/136‑1 (MQ=255) tGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTGGCACTCa > 1:420400/1‑123 (MQ=255) tGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTGGCACTCa < 2:420400/123‑1 (MQ=255) tcgtcgAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTGGCACTCAGGGTTGTTACGGGTGATGTTGGCCGCAATGTTCt < 1:410554/136‑1 (MQ=255) aCTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTGGCACTCAGGGTTGTTACGGGTg < 1:387731/96‑1 (MQ=255) aCTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTGGCACTCAGGGTTGTTACGGGTg > 2:387731/1‑96 (MQ=255) ccACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTGGCACTCAGGGTTGTTACGGGTGATGTTGGCCGCAATGTTCTGCAGCATGATCGTCTTACCGGCTTTcg < 1:281961/136‑1 (MQ=255) gcacgGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTGGCACTCAGGGTTGTTACGGGTGATGTTGGCCGCAATGTTCTGCAGCATGATCGTCTTACCGGCTTTcggcgg < 1:79698/136‑1 (MQ=255) gcTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCACGTGGCACACAGGGTTCTTACGGGTGATGTTGGCCGCAATGTTCTGCAGAATGATCGTCTTACCGGCTTTCGGCGGGGCAACGAt > 1:340891/1‑136 (MQ=255) | GGATAACCACATCTTTCTTGTGCTCGACCAGGCGCTTGGCCTTCTCGATCACCATCTCGGCAACCTGCACGTGGCGGGTTGGCGGCTCGTCGAAGGTGGAAGCAACCACTTCGCCACGCACGGTGCGCTGCATTTCGGTCACTTCTTCCGGGCGCTCGTCGATCAGCAGGACGATCAGGTGGCACTCAGGGTTGTTACGGGTGATGTTGGCCGCAATGTTCTGCAGCATGATCGTCTTACCGGCTTTCGGCGGGGCAACGAT > NC_002947/5952692‑5952953 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |