Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,556,205 | 0 | A | C | 36.7% | 36.6 / 10.1 | 30 | S215R (AGT→AGG) ‡ | rtcA | RNA 3'‑terminal phosphate cyclase |
Reads supporting (aligned to +/- strand): ref base A (18/1); new base C (0/11); total (18/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.20e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC > NC_000913/3556071‑3556343 | ccGGCCCTGACGCGCTTTTCACCGACGAAAAAAAAGCGTTCGGTAATTTTTTCCCTTCCGACTTAAAGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAAC‑CTccc < 2:1102559/139‑1 (MQ=255) cACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGcc > 1:470421/1‑140 (MQ=255) gACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtg > 1:1211959/1‑140 (MQ=255) ccccGACAAAAAAAAAGCGTTGGGTAATTTTTTCCCTTTCGACTTCAAGGGAGACGTTATCCCCCGCCCCCTGGCCGCGCGGCAGTTTATAAATTTTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCa < 2:1279919/138‑1 (MQ=255) cACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCa > 1:173582/1‑140 (MQ=255) ccGACGAAAAAAAGGCGTTCGGTGATTTTTTCCCTTTCGCCTTCAAGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACg < 2:444624/140‑1 (MQ=255) ggTTGGGTGTTTTTTTCCCTTTCGACTGAAGGAGACAAGGTATTCCCCGCCCCCTGGTCGCGGGGCATTTTATGAATTTTCTTTTGAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg < 2:611528/139‑1 (MQ=255) cGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg > 1:449668/1‑140 (MQ=255) attttttCATTTTCGACTTCAGCAGGGGCGGTATCACCCGCCCCCTGGCCCCGTGCCAGGTTATGAATATTTTTTTTATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc < 2:1093586/137‑1 (MQ=255) ttttttCTTTTTCGACTACAAGCGAGGCGTTTTTCCCCGGCCCCTGGTCGCGGGGCGGTTTATGAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa < 2:423128/138‑1 (MQ=255) ttttttCCCTTTCGACTTAAGGCGAGAGGGTTTTCCCCGCCCCCGGGTCGCGGGGCAGTTTATAAATTTTCTTTTCAGGCAGGAAAACC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa < 2:672495/138‑1 (MQ=255) cGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGAGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACt > 1:151456/1‑140 (MQ=255) cAAGCGAGACGGTATTCCCCGGCCCCGGGCCGCGAGGCGGGTTATAAATTTTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa < 2:151456/140‑1 (MQ=255) aaaGCGAGACGGTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa < 2:201311/139‑1 (MQ=255) gagcCGGTATTCCCCGGCCCCGGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt < 2:449668/140‑1 (MQ=255) gagaGGATATTCCCCGGCCCCTGGCCGCGCGGCAGGTTATGAATATTCTGTTCATGAAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt < 2:1211959/140‑1 (MQ=255) cccGGGTGGCGGGGCAGGTTAAAAATATTCTGTTCTTGCAGGAAAAGC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgcg < 1:303198/140‑1 (MQ=255) tGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc > 1:1274473/1‑140 (MQ=255) ggTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg > 1:164614/1‑140 (MQ=255) cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg > 2:721072/1‑140 (MQ=255) gagtttATAAATTTTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 2:653837/136‑1 (MQ=255) gaGGTTATGAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 2:164614/139‑1 (MQ=255) cAGTTTATAAATTTTTTTTTCAGGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 2:591545/140‑1 (MQ=255) aaGGTAATAAATTTTCTGTAAACGAAGGAAAAAC‑CACCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 1:913224/139‑1 (MQ=255) taattatTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGACGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtgat < 2:408121/136‑3 (MQ=255) atTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg > 2:443366/1‑140 (MQ=255) ttttGTTCAGGCAGGAAAACC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa < 2:303392/137‑1 (MQ=255) ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa > 1:71270/1‑140 (MQ=255) ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa > 1:654110/1‑140 (MQ=255) tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt > 1:1150302/1‑140 (MQ=255) tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt > 2:633160/1‑140 (MQ=255) tGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc > 1:478424/1‑140 (MQ=255) ttCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa > 2:323662/1‑122 (MQ=255) ttCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 1:323662/122‑1 (MQ=255) ttCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa < 2:25510/140‑1 (MQ=255) ttCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa < 2:71270/140‑1 (MQ=255) tgggaGGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg < 2:641735/136‑1 (MQ=255) aTGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg > 1:25510/1‑140 (MQ=255) tGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg > 2:897023/1‑140 (MQ=255) gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc > 1:303392/1‑140 (MQ=255) acgaaaaac‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCga < 2:1150302/136‑1 (MQ=255) aa‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc < 2:654110/140‑1 (MQ=255) | TCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC > NC_000913/3556071‑3556343 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |