Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,556,207 | 0 | T | C | 29.7% | 60.2 / 14.6 | 37 | S215G (AGT→GGT) | rtcA | RNA 3'‑terminal phosphate cyclase |
Reads supporting (aligned to +/- strand): ref base T (25/1); new base C (0/11); total (25/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.40e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT > NC_000913/3556068‑3556345 | accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑Ct > 1:1040593/1‑140 (MQ=255) gCACTGACGCGCTTTTCACCGACGAAAAAAAAGCGTTGGGTGATTTTTTCTCTTTCGACTTAAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATATTCTGTTCATGC‑GGGAAAAA‑cc‑cccccgc < 2:1543/140‑4 (MQ=255) cTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAg > 1:74485/1‑140 (MQ=255) aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGCGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAgtgt > 2:946016/1‑140 (MQ=255) cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGa > 2:194693/1‑140 (MQ=255) acgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTc > 2:1252873/1‑140 (MQ=255) cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCa > 2:683774/1‑140 (MQ=255) cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCCACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCa > 2:85092/1‑140 (MQ=255) cGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg > 2:210288/1‑140 (MQ=255) ggtattttttcCACTTTCGACTAAAGGCGAGAGGTTTTTCCCCGCCCCCGGGTCGCGCGGCGGTTTATGAATTTTCTGTCCAGGCGGGAAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca < 2:1111473/133‑1 (MQ=255) taattttttCGCTTTCGACTCAAAGCGAGACGTTTTTCCCCGGCCCCGGGCCGCGGGGCAGTTTATGAATTTTCTGTTCAGGAAGGAAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca < 2:926882/135‑1 (MQ=255) ttttttCCTTTTCGACTTAAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGGGGCAGGTTATAAATTTTTTTTTAATGGAGGAAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa < 2:1040593/138‑1 (MQ=255) ttttttCACTTTCGACTTCAAGCGAGACGGTATTCCCCGCCCCCGGGTCGCGCGGCAGGTTATAAATTTTTTGTTCAGGAAGGAAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa < 2:154335/138‑1 (MQ=255) aaaGGGAGAGGGTCTTCCCCGCCCCCTGGTCGCGGGGGAGGTTATAAATTTTCTGTTCA‑GGAGGGAAAAACCC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa < 2:728961/139‑1 (MQ=255) aaaGCGAGACGGTATCCCCCGGCCCCTGGCCGCGGGGCAGTTTATGAATATTTTGTTCAGGCAGGAAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa < 2:1295426/139‑1 (MQ=255) tATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca > 2:1021712/1‑140 (MQ=255) cccGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt > 1:297529/1‑140 (MQ=255) ccGGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTTTTCATGCAGGAAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTc < 2:1241161/140‑1 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 1:1562847/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 1:647640/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 1:194283/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc < 1:1252873/140‑1 (MQ=255) cgcgGCGGGTTATGAATTTTCTGTTCATGCAGGAAAAAC‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg < 1:194693/140‑1 (MQ=255) cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg > 2:1110096/1‑140 (MQ=255) cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg > 1:839071/1‑140 (MQ=255) cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg > 1:725143/1‑140 (MQ=255) gggtttATGAATTTTCTGTTCAGGAAGGGAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 2:74485/136‑1 (MQ=255) gagtttATAAATTTTCTGTTCATGCAGGGAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 2:194283/136‑1 (MQ=255) cAGGTTATGAATTTTCTGTTCATGCAGGAAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 1:210288/140‑1 (MQ=255) tATGAATATTCTGTTCATGCAGGAAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtg < 1:946016/140‑1 (MQ=255) aTGAATATTCTGTTAAAGCAGGGAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtgt < 1:1513510/140‑1 (MQ=255) aaaTTTTCTGTTCAGGAAGGAAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa < 2:1562847/139‑1 (MQ=255) ttCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa > 1:657101/1‑140 (MQ=255) ttCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa > 1:549851/1‑140 (MQ=255) tCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt > 1:433431/1‑140 (MQ=255) tCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt > 1:1588966/1‑140 (MQ=255) gTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc > 1:389032/1‑140 (MQ=255) ttCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAcctgtctc > 1:1480806/1‑133 (MQ=255) ttCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAcctgtctc > 1:1193935/1‑133 (MQ=255) ttCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa > 2:1560643/1‑140 (MQ=255) ttCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa > 2:647620/1‑140 (MQ=255) aatgaGGGAAAAC‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg < 2:839071/136‑1 (MQ=255) aaGGAAAAAC‑CC‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg < 2:1196529/139‑1 (MQ=255) gAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc > 1:1327107/1‑140 (MQ=255) cT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt < 2:1588966/140‑1 (MQ=255) | ACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CT‑CCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT > NC_000913/3556068‑3556345 |
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: ‑ |