| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 3,965,557 | (T)8→9 | coding (74/1266 nt) | rhlB ← | ATP‑dependent RNA helicase RhlB |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,965,549 | 1 | . | T | 85.0% | 58.9 / 4.4 | 20 | coding (82/1266 nt) | rhlB | ATP‑dependent RNA helicase RhlB |
| Reads supporting (aligned to +/- strand): ref base . (2/1); new base T (7/10); total (9/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.66e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 | |||||||||||
CGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTT > NC_000913/3965472‑3965624 | cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑ttt < 2:777069/81‑1 (MQ=255)cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑ttt < 1:777068/81‑1 (MQ=255)cgcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑tt < 1:777070/80‑1 (MQ=255) gcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑tttt > 1:777118/1‑81 (MQ=255) gcCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑tttt < 1:777071/81‑1 (MQ=255) ccTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑ttttt > 2:777109/1‑81 (MQ=255) gCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑tttttttt > 1:777105/1‑81 (MQ=255) cccGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑ttttttttt > 1:777113/1‑80 (MQ=255) cccGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑ttttttttt < 2:777073/81‑2 (MQ=255) cccGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑ttttttt < 2:777072/79‑1 (MQ=255) ccGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑tttttttttc < 1:777074/81‑3 (MQ=255) aCGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCTT‑‑TTTTTTTTCAAGGGc > 1:777112/1‑81 (MQ=255) cACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTc > 2:777119/1‑80 (MQ=255) ccGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAAcc > 1:777091/1‑81 (MQ=255) gCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGa < 1:777075/81‑1 (MQ=255) gTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATg > 1:777111/1‑81 (MQ=255) gTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATg > 1:777126/1‑81 (MQ=255) gCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAgg > 2:777122/1‑81 (MQ=255) ggAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGc < 1:777076/81‑1 (MQ=255) ggCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGt < 1:777077/81‑1 (MQ=255) gCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTc < 1:777078/81‑1 (MQ=255) gCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTc < 2:777079/81‑1 (MQ=255) aaTGGGCGTACAGTTATGAAACC‑C‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctg < 2:777080/81‑1 (MQ=255) gggCGTCCAGTTATGAACCC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt < 1:777081/81‑1 (MQ=255) cGTACAGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg > 1:777125/1‑81 (MQ=255) gTACAGTTATGAAACC‑CTTTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttctg < 1:777082/81‑1 (MQ=255) aGTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaaa < 2:777083/81‑1 (MQ=255) gTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAt < 2:777084/81‑1 (MQ=255) gTTATGAAACC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAt > 1:777124/1‑81 (MQ=255) aaaCC‑C‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtttt > 2:777129/1‑81 (MQ=255) aaCC‑CT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtttt < 2:777085/81‑1 (MQ=255) | CGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACC‑C‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTT > NC_000913/3965472‑3965624 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |