Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,423,297 | 0 | G | C | 43.4% | 83.3 / 29.5 | 53 | 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 (10/20); new base C (23/0); total (33/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.22e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCG > NC_000913/3423150‑3423446 | aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCggg < 1:686315/150‑1 (MQ=255) cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCCCCCGGCCg > 1:714769/1‑149 (MQ=255) gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGCTTCTC‑CCCCCCCCCCTCCCGGcc > 2:108976/1‑146 (MQ=255) tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCACCCTCCCGGCCGGTGCtt > 1:388594/1‑148 (MQ=255) aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCa < 1:275120/151‑1 (MQ=255) aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTG‑CCCCCCCCCCCCCCCGCCGCTGCtt > 2:286358/1‑144 (MQ=255) ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCTC‑CCCCCCCCCCCCCCGGGCGGTGGTTctccacagt > 2:53223/1‑143 (MQ=25) ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTG‑CCCGCCCCCCTCCCGGCCGGTGCTTCGCCAccattaaa > 2:688465/1‑143 (MQ=255) tttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAAt < 1:9080/150‑1 (MQ=255) aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCCCCCGGCCGTTGCTTCGCAACG‑TTTAAAAcccc > 1:366678/1‑149 (MQ=255) aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGCCG‑TGCTTTTTAACG‑TTCAAAtac > 2:177662/1‑147 (MQ=255) aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCcg < 1:290970/151‑1 (MQ=255) ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGTCCG‑CCCCCCCCCCCCCCGGCCCCTGCTTCCCCACC‑TTCAAAAtccctcccc > 2:260383/1‑149 (MQ=255) gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCACCCCCCCGGCCGTTGCTTCGCAAAG‑TTCAAATCCCCAccc > 1:215759/1‑148 (MQ=21) tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCT‑CCCCCCCCCCCCCCGGCCGGTGCTTTGCAACG‑GTCAAATCCGCcccc > 2:895254/1‑144 (MQ=21) cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCG‑CCCCCCCCCCCCCCGGCCCGTGCTTTGCAACA‑GTCAAAATCGCCCCccgccggtt > 2:811707/1‑145 (MQ=255) gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTAc < 1:785590/150‑1 (MQ=255) aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTc < 2:463048/151‑1 (MQ=255) aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCg < 1:552472/151‑1 (MQ=255) gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCggg < 1:385786/151‑1 (MQ=255) aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCG‑CCCGCCCCCCCCCCGGCCGTTGCTTTGCAACGCTTCAAA‑CCGCCCCCCGCGGGTTTGTTCTACTcccggg > 1:690313/1‑145 (MQ=17) gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTG‑CCCCCCCCCCCCCCGGCCGGTGCTTCGCAACG‑TTTAAATTCCCCCCCCGCGGGGTTTTCCTCCTCcgggg > 2:429747/1‑147 (MQ=255) tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgt < 2:838772/150‑1 (MQ=255) gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCG‑CCCCCCACCCCCCGGGCCGTGTCGTCGCCACG‑TTCAAAATCGCTCCCCGCGGGTTTTTTCCACTccgcggaggatt > 2:139822/1‑140 (MQ=255) aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCC‑CCCCCCCCCCCCCCGGCCGCTGCTTCTCCAAG‑GTTAAATCCCCCCCCCGCGGAGTTGTTCCACTCCGGGGAgagt > 1:657964/1‑148 (MQ=12) aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCCG‑CCCCCCACCCACCCGGCCGGTGCTTCGCAACG‑GTCCAATCCGCTCCCCGCGGGGTTGTGCTACTCCGGAGAGGGttt > 1:548654/1‑148 (MQ=17) aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGac < 2:694080/151‑1 (MQ=255) aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCG‑CCCCCCCCCCCCCCGGCCGGTGCTTCGCAACC‑GTCACATCCGCCCCCCGCGGGGTTGTTCTTCTCCGGGGAGGGTTTACCgcc > 2:52728/1‑149 (MQ=9) aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac < 2:150841/151‑1 (MQ=255) ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca < 1:821266/151‑1 (MQ=255) gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaacaa < 2:370690/151‑1 (MQ=255) gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCG‑CCCCCCCCCCCCCCGGCCCGTGCTTCGCAACG‑TTTAAAACCCCCCCCCGCGGGGTTTTTCTACTCCGGAGAGGGGTTACCGacaaaca > 1:184802/1‑150 (MQ=9) cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTT‑CCCCCCCCCCCCCCGGCCGTTTCTTCGCAACG‑TTCAAATCCCCCCCCCGCGGGTTTTTCCTTCTCCGGAGAGGGTTTACCCaccaaaaa > 1:475418/1‑149 (MQ=14) aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCACCCCCCCGGCCGGTTCTTCGCAACG‑TTCAAATCCGCCCCCCGCGGGTTTTTCCTACTCCGGAGAGGGGTTACCCacaaaaaa > 1:589438/1‑148 (MQ=14) tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCCCCCCCCGGGCCGTTGCTTTCCAACG‑TTCAAATCCGCCCCCCGCGGGTTTTTCCTACCCGGGAGAGTGTTTACCCACAAaccaca > 1:445181/1‑146 (MQ=14) ggAGGGCCTCTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAc < 1:571585/150‑1 (MQ=255) cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCCG‑CCCCCCCCCCCCCCGGCCGCTGTTTCGCAAAG‑TTTAAAACCGCCCCCCGCGGGTTTTTCCCTATCGGGGGGGAGGTCTCCGCCAAACAACAAAAAAAACAAAAAGCCAAggtt > 1:508793/1‑148 (MQ=255) tgatttgatGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCt < 2:114086/149‑1 (MQ=40) atttgatGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg < 2:205311/151‑1 (MQ=38) gatGTCTGGCAGTTCATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACt < 1:829481/150‑1 (MQ=38) atGCCTGGCAGTTTATGGCGGGCGGCCTG‑CCCCCCCCCCCCCCGGCCGGTGCTTTTCAAAC‑TTTAAATCCCCCCCCCGCGGGGTTTTTCTTCTCCGGGGAGGGGTTACCCACAAACAAAAAATAAAAAAAAAAGGCCCGGCTTcccgc > 2:338722/1‑145 (MQ=0) tCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCCCCCCCCCGGCCGTTTCTTTGCAACA‑GTCAAAACCGCCCCCCGCGGGTTTTTTCTACTCCGGGGAGGGGTCTCCCACAAAAAACAGAAAAAAAAAAAAGCCCAGGCTTTCGGCTGGGcc > 1:131438/1‑151 (MQ=2) aGTTTATGGCGGGCGTCCTG‑CCCGCCCCCCCCCCGGCCGGTGCTTCGCAACG‑TTCAAATCCGCTCCCCGCGGGTTTGTCCTCCTCCGGGGAGGGGTCACCGACAAAAAAAAAATAAAAAAAAAAGGCCAGGCTTCCCACTTGGCCcttttt > 2:999981/1‑149 (MQ=2) tATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt < 2:43270/150‑1 (MQ=38) ggcgggcgTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt > 2:429807/1‑151 (MQ=38) ggcgTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATg > 1:335800/1‑151 (MQ=38) gTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt > 2:679118/1‑151 (MQ=40) g‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc < 1:453485/149‑1 (MQ=40) g‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc > 1:425888/1‑149 (MQ=40) g‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa < 2:98371/150‑1 (MQ=40) ccGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGtt < 2:421880/151‑1 (MQ=40) cACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt > 2:102610/1‑151 (MQ=40) cACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt > 2:730549/1‑151 (MQ=40) cGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCg > 2:668911/1‑151 (MQ=40) | AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCG > NC_000913/3423150‑3423446 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |