Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,423,290 | 0 | A | C | 42.3% | 56.9 / 17.9 | 54 | intergenic (+96/+133) | yhdZ/rrfF | putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon |
Reads supporting (aligned to +/- strand): ref base A (7/23); new base C (20/2); total (27/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTG > NC_000913/3423149‑3423429 | tAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCT‑GCCCCCCCCCCAccc > 1:867770/1‑148 (MQ=255) cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCT‑TCCCCCCCCCCCCCCGGcc > 1:11761/1‑148 (MQ=255) cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGt < 2:775733/150‑1 (MQ=255) cGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑TCCCCCCCCCCCcccgcccggcgttcc > 2:161786/1‑144 (MQ=37) gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTCCCCCCCCCCCCCCCGGGCCG‑TGCTTCGCa > 2:951031/1‑149 (MQ=255) aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTGCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑g < 1:333847/151‑1 (MQ=255) ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCC‑TGCCGCCCCCCTCCGTGGCGGTGTTTCGCCAcagt > 2:896252/1‑148 (MQ=25) aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGACTCCT‑GTCCCCCCCCCCCCCCGCCGTTGCtttctacactt > 1:432364/1‑139 (MQ=255) aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCg < 2:960092/150‑1 (MQ=255) gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGTTCAT‑CCCCCCCCCCCCCCCCGCCGTGTCTTTTCAAC‑GCTCAACTcaccccc > 1:796416/1‑145 (MQ=255) gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTcc < 2:482645/151‑1 (MQ=255) aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCcg < 1:440562/151‑1 (MQ=255) aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGTCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCcggcg < 1:947978/150‑1 (MQ=255) ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGAt < 1:823727/150‑1 (MQ=255) cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCC‑TCCCCCCCCCCCCCCGGCCGGTGCTTCGCAAC‑ATTCAAAACCGCTCCccgccggctt > 1:231110/1‑145 (MQ=255) aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCt < 2:360172/151‑1 (MQ=255) aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCt < 2:68758/151‑1 (MQ=255) tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTa < 2:401103/150‑1 (MQ=255) ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCggg < 1:4922/150‑1 (MQ=255) ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCC‑GCCCCCCCCCCCCCCGGCCGTTTCTTTTCAAC‑GTTCAAATCCGCCCCCCGCGGGTTTGTCCCACTCggg > 1:922490/1‑150 (MQ=17) ttAGTACGAAAAGCAAAAAGGCCATCCTGTCGGATGGCCTTTCGCTTGATTTGAGGCCTTGCAGTTTATGGCGGTAGCGCT‑GCCCGCCACCCTCCGGCCCGTTGCGTCGCAAC‑GTTCAAATCCGCGCCCGGCGGATTTGTTCTa < 1:300756/144‑1 (MQ=255) tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCT‑CCCCCCCCCCCCCCCGGCCGTTGCTTCGCAAC‑GCTCAAATCCGCTCCCCGCGGGGTTTTCCtcccccgggggggg > 1:964917/1‑138 (MQ=17) cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑CCCCCCCCCCCCCCCGGCCGTTGCTTTGCAAA‑CGTCAAAAACGCCCCCCGCGGGTTTTTTCTACTCcgggggggg > 2:653670/1‑144 (MQ=17) gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGtt < 2:308054/151‑1 (MQ=255) aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAcc < 1:653603/151‑1 (MQ=255) gCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCg < 2:511073/151‑1 (MQ=255) aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTTCT‑TCCCCCCCCCCCCCCCGCCGCTTCTTCCCCAC‑GCTTCAATCCCCCCCCCGCGGGTTTTTTCCCCTCTGGTGAGTGTTCTccccc > 2:123959/1‑148 (MQ=12) aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCg < 1:744101/149‑1 (MQ=255) aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCC‑CCCCCCCCCCCCCCGGGCCGTTGCTTCGCTAC‑GGTCAAATACGCCCCCGGCGGGGTTTTCCCACTCCGGAGAGGGTTCTCCCac > 2:843920/1‑151 (MQ=14) aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCC‑TCCCCCCACCCCCCCGGGCGCTGCCTTGCAACGGTTCAATCCCCCCCCCGCCGGGTTGTTCCT‑CTCGCGagaggggtttccccc > 1:741968/1‑140 (MQ=9) aaaaGGCCATCCTTTCGGATGGCCTGTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGaca < 1:799081/151‑1 (MQ=255) aaaGGCCTGACTCTCGGATTGCCTTACTCTTGTTTTGATGTTTGACTGTTTATGACGGTAGTCTT‑GTCCGCCGCCCTCCGGGCCGTTGCTTCGCATC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGac < 1:912671/149‑1 (MQ=14) ccttcttttCGTGTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCAGCCCGCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAaacaac < 1:162198/145‑1 (MQ=14) aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCT‑TCCCGCCCCCCTCCCGGGCGTTGCTTTGCAACGGTTCAAAT‑CGCCCCCCGCGGGTTTGTCCTACCCCGGAGAGGGGTCACCCacaaaaaaaa > 2:298524/1‑148 (MQ=255) cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCT‑TCCCGCCCCCCTCCCGGCCGTTGCTTCTCCAC‑GTTTAAATCCTCTCCCCGCGGGTTTTTCCTACTCCGGGGAGGGTTCCCCCacaaacaa > 1:21746/1‑145 (MQ=14) tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCC‑GCCCGCCCCCCCCCCGGCCGTTGTTTCGCCAC‑GTTCAAAACCGCTCCCGGCGGGGTTGTTCTACTCGCGGGAGGGGTCACCCACAAAACACAGATa > 1:370560/1‑150 (MQ=255) aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaa < 2:361996/151‑1 (MQ=255) aTGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGTCCT‑GCCCCCCCCCCCCCCGGCCGGTGGTTTGCAAC‑GATCAAATCCGCCCCCCGCGGGTTTTTTCTTCACCAGCGAGGGTTTACCCACCAACAAAAACTaaaaa > 1:599427/1‑147 (MQ=0) ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAg < 2:262719/151‑1 (MQ=255) cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑TCCCCCCCCCCCCCCGGCCGTTGCTTCGCAAC‑GTTCAAATACCCCCCCCGCGGGTTTTTCCTACTCCGGAGAGTGTTCACCCACAAACAACAGATAAAACAAAAAGcc > 2:841628/1‑151 (MQ=14) cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCCCCCACCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGt < 1:356683/151‑1 (MQ=255) atttgatGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCCCCCCCCCCCGGGGCGGTGCTTCCCAAC‑GTTTAAAACCCCTCCCCGCGGGGTTGTCCCACCCGGGAGAGAGTTTCCCCACAAACCACAGATAAAAAAAAAAGCCCAGTCTTcc > 1:81579/1‑150 (MQ=2) atGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGa < 2:614733/151‑1 (MQ=38) ggCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt < 2:374764/151‑1 (MQ=35) aGTTTATGGCGGGCGTCCT‑GCCCCCCACCCTCCCGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCCCCCGGCGGGTTTGTCCCACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAAGCCCAGTCTTCCGACTGAGCCTTTTGt > 2:80486/1‑151 (MQ=14) aGTTTATGGCGGGCGTCCC‑TCCCGCCCCCCCCCCGGCCGTTGCTTTGCAAC‑GTTCAAATCCGCTCCCGGCGGGGTTTTCCCACCCGGGAGAGGGGTCACCCACCAAAAAAAGATAAAAAAAAAAACCCAGTCCTCCGGCTGGGCCcttttt > 1:333599/1‑149 (MQ=2) gTTTATGTCGGGCGTCCT‑GCCCGCCCCCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGt < 1:747772/150‑1 (MQ=255) tGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGGGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt > 1:932705/1‑151 (MQ=38) tGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt > 1:559542/1‑151 (MQ=38) ggcgggcgTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt < 1:941537/151‑1 (MQ=38) gggcgTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt > 2:549865/1‑151 (MQ=38) ggcgTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATg < 2:197009/151‑1 (MQ=38) ggcgTCCT‑GACCGCCACCCTCCGGGCCGTTACTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTCTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGACCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATg > 1:71242/1‑151 (MQ=255) tCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg > 1:421428/1‑151 (MQ=40) | TAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑GCCCGCCACCCTCCGGGCCGTTGCTTCGCAAC‑GTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTG > NC_000913/3423149‑3423429 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |