Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,423,297 | 0 | G | C | 48.2% | 64.1 / 34.0 | 56 | intergenic (+103/+126) | yhdZ/rrfF | putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon |
Reads supporting (aligned to +/- strand): ref base G (4/25); new base C (27/0); total (31/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.32e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTC > NC_000913/3423153‑3423445 | tCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCcgg > 1:431682/1‑147 (MQ=255) agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTCCCGCCACCCTCCCGGGCGtt > 1:496326/1‑148 (MQ=255) gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGCTCCCTCCCGCCCCCCCCCCGGccccttgt > 2:339764/1‑142 (MQ=255) gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGcc > 1:456993/1‑142 (MQ=255) cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCACCCCCCCGGCCGGTGCt > 2:12164/1‑145 (MQ=37) gCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTc < 2:192734/151‑1 (MQ=255) tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCCCCCCCCCGGCCGGTGCttttcaaaggttaa > 1:528040/1‑138 (MQ=255) ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCGTTGTTTTGCCA‑CGTTCaaa > 2:121236/1‑150 (MQ=255) ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATc < 1:389133/150‑1 (MQ=255) cGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAA‑CGTTCAAATcccccc > 2:133454/1‑147 (MQ=255) ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTGCTGCCCGCCACCCTCCGCGGCCTGGCTTTGCAACCTTTCAAAT‑CTCCCCCCGCGGGTTTGt > 1:449198/1‑151 (MQ=255) ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCACCCCCCCGGCCGGTGCTTCGCAA‑ACTTCAAATCCGCCCCCCGCGGGTTTGt > 1:130040/1‑151 (MQ=18) tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGt < 2:172073/150‑1 (MQ=255) tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCAA‑CGTTCAAATCCCCTCCCCGCGGGttttt > 2:342092/1‑148 (MQ=18) aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCACCCTCCCGGCCGTTGCTTTCCAA‑CGTTCAAAACCGCCCCCCGCGGGGTTGt > 1:446920/1‑149 (MQ=18) ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCCCAA‑CGTTTAAATCCCCCCCCCGCGGGttttttc > 1:632376/1‑145 (MQ=18) tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCCCCCCCCCCCCCCGGCCGTTTCTTCGCCA‑AGTTTAAAAACGCCCCccgcggggtttt > 1:168285/1‑140 (MQ=255) ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCCCCCCCCCGGCCGGTGCTTCGCAA‑CGTTCAAAATCGCCCCCCGCGGGTTTGTTCtcccccgg > 1:300869/1‑143 (MQ=255) tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGga < 1:176765/150‑1 (MQ=255) tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGGTGCTTTGCAA‑CGTTCAAATCCGCCCCCGGCGGGGTTGTCctcctc > 1:426723/1‑146 (MQ=17) gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGTTTCTTCGCAA‑CCTTTCAATCCGCTCCCCGCGGGTTTTTCCTACTCcgg > 2:95826/1‑147 (MQ=17) cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCACCCTCCCGGGCGTTGCTTCGCAA‑CCATCAAATCCCCTCCCGGCGGGTTTGTCCTTCTCCGga > 2:332353/1‑145 (MQ=255) gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCGCCCGCCCCCCCCCCGGCCGGTGCTTCGCAA‑CGGTCAAAACCGCTCCCGGCGGGGTTTTTCTACTCCGGGGAGGGtt > 1:206268/1‑151 (MQ=255) gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGtt < 2:605412/151‑1 (MQ=255) gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAA‑CCGTCCAATACCCCCCCCGCGGGGTTGTTctcctcccgggaggggt > 1:483312/1‑140 (MQ=255) aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGGTGCTTTGCAA‑ACGTCAAATACGCTCCCCGCGGGTTTTTTCTCCTCcgggg > 1:496204/1‑141 (MQ=12) aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAc < 2:94865/151‑1 (MQ=255) aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGGCGGTGCTTCGCAA‑CGTTTAAAACCCCTCCCCGCGGGGTTGTCCTACTCcgggggggggttac > 1:367374/1‑140 (MQ=14) cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGa < 2:196687/151‑1 (MQ=255) aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGGTGCTTCTCAA‑CGTTCAAATCCGCCCCCGGCGGAGTTTTCCTACTCCGGAGAGGGTTCAcc > 2:536218/1‑147 (MQ=14) aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac < 1:327886/151‑1 (MQ=255) gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaacaa < 2:434309/151‑1 (MQ=255) ttCGGATGGCTTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaa < 1:180160/150‑1 (MQ=255) tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaaa < 1:12162/150‑1 (MQ=255) tGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCTGTTCTTTCCCCCCCCCTTCCGGCCGCATCCTCGCAA‑CGGTCAAATTCCCCCCCGGCGGATTTTTCCTACTCCGGGGAATGTTTCCCCGCAAACAACACGTaaaaaaaa > 2:655185/1‑150 (MQ=255) ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaaa < 1:339725/150‑1 (MQ=255) cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGTTGCTTCGCCA‑CGTTTAAATCCCCCCCCCGCGGGTTTTTTCTACTCGGGAGAGGGTTTACCCCCCAACAACAGATAaaaaaaaaa > 2:23731/1‑147 (MQ=9) ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCCGGCCGTTGCTTTGCAA‑CGTTTAAATCCGCCCCCCGCGGGTTTTTTCTCCCCGGGGGGGGGTTCACCCACCAACAACAGAAAAAACAAAAAGccc > 1:9123/1‑150 (MQ=14) ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCTCAA‑AGTTTAAATCCGCCCCCCGCGGGTTTGTCCTACTCGGGAGAGGGTTTACCCACCAACAACAGATAAAAAAAAAAGccc > 1:14404/1‑150 (MQ=255) cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg < 2:528121/150‑1 (MQ=255) cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCTCCCGGCCGGTGCTTCGCCA‑CGTTCAAATCCGCTCCCCGCGGGTTTTTCCTACTCCGGGGAGGGTTTACCGACAAACAAAAGAAAAAACaaaaaagccacgc > 1:542807/1‑141 (MQ=11) atttgatGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCGGGCCGGTGCTTCGCAA‑CGTTCAAAAACCCCCCCCGCGGGGTTTTCCTACTCGGGAGAGTGGTTACCCACAAACAAAAGATAAAAAAAAAAGGCCAGGCtttc > 2:120991/1‑148 (MQ=2) gatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTg < 2:509098/151‑1 (MQ=38) gTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAg < 2:247459/150‑1 (MQ=38) gTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAg < 2:638329/150‑1 (MQ=38) gCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt < 2:637414/151‑1 (MQ=37) gCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCCGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCCCCCGGCGGATTTGTCCTACTCCGGAGAGTGGTCACCGACAAAAAAAAGATAAAACAAAAAGCCCAGTCTTCCGACTGAGCCtttt > 2:1590/1‑151 (MQ=11) aTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt < 2:122096/151‑1 (MQ=38) aTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt < 1:399883/151‑1 (MQ=38) ggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATg < 1:101806/151‑1 (MQ=38) gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt < 2:493003/151‑1 (MQ=40) gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt < 1:128779/151‑1 (MQ=40) ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg < 1:234690/150‑1 (MQ=40) cGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTc > 1:606232/1‑151 (MQ=40) gCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc < 2:41976/151‑1 (MQ=40) cTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAc < 1:310992/149‑1 (MQ=40) ccGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACtctc > 1:215420/1‑151 (MQ=40) | TCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTC > NC_000913/3423153‑3423445 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |