| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 3,965,557 | (T)8→7 | coding (74/1266 nt) | rhlB ← | ATP‑dependent RNA helicase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,965,550 | 0 | T | . | 97.0% | 129.3 / ‑2.7 | 33 | coding (81/1266 nt) | rhlB | ATP‑dependent RNA helicase |
| Reads supporting (aligned to +/- strand): ref base T (0/0); major base . (13/19); minor base C (1/0); total (14/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.24e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCG‑TTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAA > NC_000913/3965413‑3965683 | gagagaAGATAATGAAACGTTGACGTAAGAAACGCCATCG‑TTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtt < 2:417675/139‑1 (MQ=255) aagaTAATGAAACGTTGACGTAAGAAACGCCATCG‑TTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttttttt > 2:690270/1‑139 (MQ=255) agaTAATGAAACGTTGACGTAAGAAACGCCATCG‑TTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCtttttttc > 1:1010935/1‑138 (MQ=255) aTGAAACGTTGACGTAAG‑AACGCCATCGTTTTTCCGGGAACGGGTTTGCCCTTGCCCGGC‑AAGTCCCCACCCGCCAGCCTCAACGGGAGGGCCAGTGCCTTAATGGGCGTACAGTTATGAAACCCCTTTTTTTAAAggc > 2:33538/1‑138 (MQ=255) gaaacctttGACG‑AAGAAACGCCATCG‑TTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCt < 1:33538/132‑1 (MQ=255) gTTGACGTAAGAAACGCCATCG‑TTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCtt < 1:151499/135‑1 (MQ=255) gTTGACGTAAGAAACGCCATCG‑TTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCtt > 2:151499/1‑135 (MQ=255) gTTGACGTAAGAAACGCCATCG‑TTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTAc > 1:208879/1‑139 (MQ=255) gTTGACGTAAGAAACGCCATCG‑TTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTAc > 2:276921/1‑139 (MQ=255) ttGACGTAAGAAACGCCATCG‑TTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACa < 2:896525/139‑1 (MQ=255) tCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGAGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCg > 1:764946/1‑139 (MQ=255) cGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa < 2:1010935/139‑1 (MQ=255) ccGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACtt < 1:276921/139‑1 (MQ=255) ccGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACtt < 2:208879/139‑1 (MQ=255) tGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttc > 2:902184/1‑139 (MQ=255) cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg > 2:1010191/1‑139 (MQ=255) tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGTAAACTTCTGTTCTGTTAAATGTGTTTTGCt < 2:182989/139‑1 (MQ=255) tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCt < 2:535598/139‑1 (MQ=255) cACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTc > 2:43371/1‑139 (MQ=255) gACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATa < 2:546079/139‑1 (MQ=255) ccAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTgg < 1:273805/139‑1 (MQ=255) cAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGt < 1:690270/139‑1 (MQ=255) tCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGttt < 1:493851/119‑1 (MQ=255) tCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGttt > 2:493851/1‑119 (MQ=255) ggAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAg < 1:902184/139‑1 (MQ=255) aaGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt < 2:20562/96‑1 (MQ=255) aaGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt > 1:20562/1‑96 (MQ=255) gAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTAc > 1:91510/1‑139 (MQ=255) gaatggCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCgg < 1:53180/49‑1 (MQ=255) gaatggCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCgg > 2:53180/5‑53 (MQ=255) cAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATc < 1:43371/139‑1 (MQ=255) aaaCCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGa > 2:455238/1‑46 (MQ=255) aaaCCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGa < 1:455238/46‑1 (MQ=39) aaaCCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGc < 1:1055759/82‑1 (MQ=255) aaaCCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGc > 2:1055759/1‑82 (MQ=255) aaaCCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGaaa < 2:91510/139‑1 (MQ=255) | GAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCG‑TTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCATAGTGTGGTAGAATATCAGCTTACTATTGCTTTACGAAAGCGTATCCGGTGAAA > NC_000913/3965413‑3965683 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |