Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,556,205 | 0 | A | C | 41.9% | 54.6 / 19.7 | 43 | S215R (AGT→AGG) ‡ | rtcA | RNA 3'‑terminal phosphate cyclase |
Reads supporting (aligned to +/- strand): ref base A (24/1); new base C (0/18); total (24/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.12e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCG > NC_000913/3556071‑3556344 | tCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTccc > 2:20829/1‑140 (MQ=255) ccggccttGACGCGCTTTTCACCGACGAAAAAAAAGCGTTCGGTGATTTTTTCACTTTCGCCTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCA‑GGAAAAAcccccc < 2:514597/133‑1 (MQ=255) gCACTGCCGCGTTTTTCACCGACGACAAAAAAGCGTTGGGTGATTTTTTCCCTTTCGACTTCAAGGGAGCCGGTATTCCCCGCCCCCGGGCCGCGGGGCAGTTTATAAATATTCTGTT‑CATGAGGGAAAAAACCCCCCGc < 2:114835/140‑1 (MQ=255) aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAgtgt > 1:1263190/1‑140 (MQ=255) aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAgtgt > 1:646704/1‑140 (MQ=255) cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCATGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGa > 2:120979/1‑140 (MQ=255) cACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCa > 1:1730289/1‑140 (MQ=255) gacgacAAAAAGGCGTTCGTTGATTTTTTCTCTTCCGACTTCAAGGGAGACGTTATCCCCCGCCCCCTGGTCGCGGGGCAGTTTATAAATATTTTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCt < 2:1314092/140‑1 (MQ=255) aaaaaaaGCGTTCGGTGATTTTTTCTCTTTCGACTTAAGGGGAGAGGGTATTCCCCGCCCCTTGGTCGCGCGGCGGTTTATGAATTTTTTTTT‑CAGGCAGGAAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCa < 2:80641/139‑1 (MQ=255) ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg > 1:863855/1‑140 (MQ=255) ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg > 1:1515443/1‑140 (MQ=255) tCGTTGTTATTTTCTTTTTCGACTTCAAGAGAGCCGGTGTCACCCGGCCCCCCGGCGCGGGGAAGTTTAGTAATTTTTTGCT‑GATGCGGGGAAAATCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg < 2:764534/140‑1 (MQ=255) ggTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTT‑CATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca < 1:531798/140‑1 (MQ=255) tattttttCTTTTTCGATTTCAAGGGAGCCGTTTTTCCCCGCCCCCGGGTCGCGCGGCAGGTTATAAATTTTCTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac < 2:1587717/136‑1 (MQ=255) ttttttCTCTTTCGACTAAAAGCG‑GGCGGTTTTTCCCCGCCCCCGGGCCGGGGGGCGGTTTTAAAATTTTTGTTAAATGGGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa < 2:423306/138‑1 (MQ=255) tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa > 1:1693716/1‑140 (MQ=255) ttttCGACTTCAAGCGAGACGGTATTCCCCGCCCCCTGGCCGCGCGGCAGTTTATGAATATTCTGTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg < 2:528881/139‑1 (MQ=255) ttCGACTTCAAGGGAGACGGTATTCCCCGCCCCCTGGTCGCGGGGCAGTTTATAAATTTTCTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa < 2:1468693/140‑1 (MQ=255) ttCGACTTAAAGCGGGACGTTTTCCCCGGGCCCCTGGTCGCGGGGCAGTTTATAAATTTTCTGTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa < 2:1391384/140‑1 (MQ=255) aCTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc > 1:1065179/1‑140 (MQ=255) cAAGCGAGCGGTTATCCCCCGCCCCCGGGCCGCGCGGCAGTTAATAATTTTTCTGTC‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa < 2:863855/140‑1 (MQ=255) aaaGCGAGACGGTATTCCCCGGCCCCTGGTCGCGGGGCTGTTTATGAATTTTTTTTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa < 2:1543575/139‑1 (MQ=255) agggagaCGGTATCCCCCGGCCCCGGGTCGCGCGGAAGTTTATAATTTTTTTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg < 2:381976/137‑1 (MQ=255) gagcgggTATCCCCCGGCCCCGGGCCGGGGGGAGGTTTATAAATTTTCTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt < 2:1730289/135‑1 (MQ=255) gaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt > 2:1307005/1‑140 (MQ=255) ggTTTTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca < 1:1594844/140‑1 (MQ=255) tACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATg > 1:1423621/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGGGATTTCACGCTCAGCAACATGGCTCGGCCCCCGACCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 2:1083691/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 2:213816/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 2:1412072/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCCCCACGCATCTGCACAATGTTCCcgc > 1:1715265/1‑140 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGAGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc > 1:1360271/1‑140 (MQ=255) ccccGGGCCGCGGGGCAGGTTATAAATTTTCTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc < 2:793549/140‑1 (MQ=255) ggTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg > 1:149148/1‑140 (MQ=255) cgcgGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg > 1:1080043/1‑140 (MQ=255) gCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGc > 1:614657/1‑140 (MQ=255) gaGGTTATGAATTTTCTGTT‑CAGGGGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 2:973954/139‑1 (MQ=255) aagtttATAAATTTTCTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 1:1695622/136‑1 (MQ=255) aaGGTTATAAATATTCTTTT‑CATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa < 1:1144456/139‑1 (MQ=255) aGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa > 1:580202/1‑140 (MQ=255) ttCTGTT‑CATGGAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa < 1:1748869/140‑1 (MQ=255) tCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACACTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt > 2:1696652/1‑140 (MQ=255) ttGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg < 1:1402342/139‑1 (MQ=255) ttt‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc < 2:1423621/139‑1 (MQ=255) tt‑CATGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa < 2:1515443/140‑1 (MQ=255) tt‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa > 2:1187468/1‑140 (MQ=255) tGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg > 2:78232/1‑140 (MQ=255) ggaaGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg < 2:549597/137‑1 (MQ=255) ggCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg < 2:1080043/139‑1 (MQ=255) gCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCGACGGGc > 1:863508/1‑140 (MQ=255) gAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc > 2:574192/1‑140 (MQ=255) aCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCg > 2:785628/1‑140 (MQ=255) | TCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCG > NC_000913/3556071‑3556344 |
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: ‑ |