Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | NC_000913 | 3,620,870 | 0 | T→G | 100.0% | 0.1 | 5 | M560R (ATG→AGG) | rhsB | Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (0/3): ref base (2/0): total (2/3) | ||||||||||
Fisher's exact test strand distribution p-value = 1.00e-01 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
TCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCAC > NC_000913/3620816‑3621155 | tCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAACCGGGACACCGTACGACGCGTGGGGAAAGGCCGTCCGTACCCCGCAGGGCGGGCCGACGCCCCGTATGGCATACggttgcgaccggc > 1:72264/1‑288 (MQ=11) tCACCGACTCGTCCGGTTAATGAACGCGTTATGACCTTGACCGCTGCGAGAAGAGGAAGGCGGGGCACCGGGAGAAAGGGCTGAGTCGGTACCGCGCATACGACAACCGTGGAGAGTGAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGAc < 2:6076/300‑1 (MQ=11) ttATGGACAGGACCGGTTCGGGCAGAGGGCGGCGGGGCACCGCGAGGAAGGGCTGAGTCAGTACCGGGCATACGACAGCCGTGGACAATTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACgg < 1:183380/300‑1 (MQ=2) gaccatgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGCCCGCCGTCATTGCCCCGGACGGCAGCCGAAACAGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGGACCACGCAGGACGGGCAGACGGGGAAACTGGAATAAGATAATGCAGGCAGGATCACCCGACAGACCACTGAAAAACGCAGc > 2:375672/1‑300 (MQ=11) acGCTGCGGGCAGAGGACGGAGGGGGACGGCGAGGAAGGGCGGAGTCAGGACCGCGCATATCACAGCCGTGGACAGTTAATTGCCGGGAAAGACACGCAGGTCCAGGACACGCGGTATGAAGACAACATGGCCGGGGACCTGACGGGCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcaccac < 1:163286/299‑1 (MQ=11) | TCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCAC > NC_000913/3620816‑3621155 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
Reads not counted as support for junction |
read_name Not counted due to insufficient overlap past the breakpoint. |
read_name Not counted due to not crossing MOB target site duplication. |