| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,620,869 | 0 | A | G | 47.4% | 2.3 / 22.7 | 19 | M560V (ATG→GTG) | rhsB | rhs element protein RhsB |
| Reads supporting (aligned to +/- strand): ref base A (5/5); new base G (7/2); total (12/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.50e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAGCCGTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGA > NC_000913/3620790‑3621167 | gAGCCGTTACGGTCAGTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAATTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCAGGTGACCTGAACGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcgc > 1:408388/1‑299 (MQ=11) gTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGATTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAg < 2:4940/298‑1 (MQ=14) gTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTa > 1:452808/1‑300 (MQ=31) cAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGAAACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGa > 2:281133/1‑296 (MQ=14) aattGCTTAGCTTCACCGACTGCTCCGTTTAGGTAACCCGGTGTGACCATTACCGCTTCGGGCAAACGACGGTGGTGCACCGCGAGGAAGGGTTGAGTCAGTGCCGCGCATGCGATAGCCGTGGACAGGTAATTGCC‑TGTAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCTTCATTGCCCCGGACTGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATAc < 2:55885/297‑1 (MQ=11) aGTTGCTGAGCTTCACCGACTGTTCCTGTTATGTAACCCGTTATGACCATGACCGTTTCGGGCAGATGACGGCGGTGCATCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAACCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAATACGGTATGAATACAACGTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCATAAACTGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATAc < 2:55897/299‑1 (MQ=11) aGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATAc < 2:707099/299‑1 (MQ=31) gTTGCTGAGCTTCACTGACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAACACGa > 2:442793/1‑300 (MQ=255) gctgaactccgcggACTGTTCTTGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACTGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCAGGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGCGGAAAGGCCGTCCGTACCACGCAGGGCTGGCTGACGCGCAGTATGGAATACAAtgc < 2:356476/287‑1 (MQ=11) actgACTGTTCCGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACgg < 2:635090/297‑1 (MQ=9) cGGTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTAGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCCGTATGGAATACGATGCTTCCGGACGGGTCATCCGcat > 2:616475/1‑298 (MQ=11) ggTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGcc < 1:281133/298‑1 (MQ=31) gTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTCATTGCCGTG‑AAAGACACGCAGGGCCATGACACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGNGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGCTGCTGCCGGACGGGTCATCCGCCTGa > 1:518236/1‑300 (MQ=11) gTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGNGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGa > 1:154168/1‑300 (MQ=9) gTTATGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGAAGGCAGCAGAAAAGGGACACAGTACGATGCGTGGNGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGACTGa > 1:517080/1‑300 (MQ=11) aTGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCa > 2:414802/1‑300 (MQ=31) aTGTAACCCGCTATGACCATGACCGTTTTGGTCAGGTGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCa > 1:59882/1‑300 (MQ=9) cGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCAGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCAGTCCGTCCCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAAc > 2:343302/1‑300 (MQ=14) gATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGa > 1:608116/1‑300 (MQ=31) | GAGCCGTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTG‑AAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGA > NC_000913/3620790‑3621167 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |