Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 941,713 | 0 | T | G | 43.3% | 89.4 / 18.6 | 60 | V252G (GTG→GGG) | dmsA | dimethyl sulfoxide reductase, anaerobic, subunit A |
Reads supporting (aligned to +/- strand): ref base T (5/29); new base G (26/0); total (31/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.48e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC > NC_000913/941571‑941860 | tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGGGCTGTTTGGCAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCtttc > 1:209442/1‑148 (MQ=255) cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTa < 2:421999/151‑1 (MQ=255) ctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTtattat > 2:309202/1‑151 (MQ=255) tccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAatat > 2:799825/1‑150 (MQ=255) ccGCGCAATTTTCGGAAGGTTTGAACTATACCGACGGCGGGTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTTGTAATAACCCTGGCGAAACGCGAATGAGTTGCGGTGGGGTGACTTACTATCt < 1:453584/151‑1 (MQ=255) gcgcAAATTGCGGAAGGTTTGAAGTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTg < 1:1030742/151‑1 (MQ=255) cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGa < 2:324209/151‑1 (MQ=255) gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGaa < 2:646408/151‑1 (MQ=255) cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc < 2:845956/151‑1 (MQ=255) cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc < 2:536403/151‑1 (MQ=255) cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTaa > 2:559174/1‑150 (MQ=255) gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTAATTTCTTTAAAAAGCAc > 1:721896/1‑148 (MQ=255) ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAGGCACGcc > 2:571763/1‑148 (MQ=255) ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa < 1:1007812/150‑1 (MQ=255) ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGa < 1:945466/151‑1 (MQ=255) ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAGGCACGCCAGAAATCTAAAGCACGCCTg > 1:322662/1‑150 (MQ=255) ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAAGCACGCCCGAAATCTAATGCCCGCATGatc > 2:408338/1‑150 (MQ=255) gcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTGAACAGGGACGCCCGAAATCTAATGCCCGCACGatc > 1:464246/1‑149 (MQ=255) gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcat < 1:625315/151‑1 (MQ=255) ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc < 2:484013/151‑1 (MQ=255) gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAAGCACGCCAGAAAACTAATGCCCGCATGATCATCATCg > 2:246865/1‑151 (MQ=255) aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAt < 2:89966/151‑1 (MQ=255) aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAt < 1:215078/151‑1 (MQ=255) aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCATTTTATCTTTAAAAAGCACCCCCGAAATCTAATGCCCCCAAGatc > 1:945904/1‑142 (MQ=255) cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATTTTTTACAGGCACCCCCGAAAACTAAAGCCCCCCTGATTATCAACggatcgc > 1:752693/1‑143 (MQ=255) cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTAATATTTTGAAAAAGCACGCCAGGAATCCAATGCCCGCCTGAGCATCATCTGTCCCCGCtaatac > 1:615119/1‑147 (MQ=255) ggATAGCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCGTGATCATCATCGTTCCGCGCTATACCGACAcc < 1:67179/151‑1 (MQ=255) atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATCCCGACACCgg < 1:596228/151‑1 (MQ=255) atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCgg < 2:790467/151‑1 (MQ=255) atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAAAGGCACGCCAGAAAACTAATGCCCGCATGGTCCTCATCGATCCGCGCTATACCGACACCgg > 1:697647/1‑151 (MQ=255) atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGt < 2:662974/150‑1 (MQ=255) aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCgg > 2:40863/1‑151 (MQ=255) aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATTCGCGCCATACCGACACCGGTGCCggg > 2:176228/1‑151 (MQ=255) tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAACGACAACGGGGCCggg > 2:73015/1‑145 (MQ=255) aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGaa < 2:800616/150‑1 (MQ=255) aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAg < 1:84788/151‑1 (MQ=255) cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTACTATCTTTAAACGGCAAGCCAGAAAATTAATGGACGCAAGGTCCTCCACGATCCGCACTACACTGACAACCGGGCCgg > 1:819210/1‑140 (MQ=255) tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCCCCTGATCATCCTCGATCCGCGCTATACCCACACCCGGGCCGGGCGCGCAGATg > 2:509668/1‑150 (MQ=255) ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAg < 2:249253/151‑1 (MQ=255) gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTGGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTTCCGGGCGCGAAGATGAGt < 1:165181/151‑1 (MQ=255) gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGACACGCCCGAAATCTAATGCCCGCATGATCATTATTGATCCGCGCCATACCGACACCCGGGCCGGGCGCGAAGAagag > 1:427240/1‑150 (MQ=255) aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTg < 1:662452/150‑1 (MQ=255) tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCCGAAAACTAATACCCCCATGATCATCATCGATCCGCGCCAAACCCACACCGGGGCCGGGCGCGAAGATGAgggg > 2:269688/1‑149 (MQ=255) gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTGTACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt < 1:1007608/151‑1 (MQ=255) gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATccc < 1:479097/150‑1 (MQ=255) tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCATCCTCGATCCCCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTccc > 2:843097/1‑149 (MQ=255) tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt < 2:467921/150‑1 (MQ=255) gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTATAATCTTGAACAAGCACGCCAGAAAACTAATGGCCGCATGATAATCATCGATCCGCGCTATACCGACACCCGGGCCGGGCGCGGAGAGGAGTGGGGCCCTAt > 1:5573/1‑148 (MQ=255) taaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTTCTATCCTGAAAAAGCACGCCAGAAATATAATGCCCGCCTGCACATCATTACTCCGCGCCATACCGACACCGGGGCCCGGCGCGCAGATGAGAGGGTCCCTCTTTGTcc > 1:810581/1‑150 (MQ=255) cccTGGCGAAACGCGAATGAGTGGCGGTTGGGGGGCTTACTTTCTTGAACAGGCACCCCCGAAAACTAATACCCCCCTGATCATCCTCGATCCCCGCCATACCGGCACCGGGGCCGGGGGCGAAGATGAGTGGATccccctt > 2:757510/1‑138 (MQ=255) cccTGGCGAAACGCGAATGAGTGGCGGTGGTGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATGGATCCGCGTTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg < 1:788700/150‑1 (MQ=255) cccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAAAAGCACCCCCGAAATCTAAAGCCCGCCTGATCATCATCGATCCCCGCTATACCGACACCCGGGCCGGGGGCGAAGATGAGTGGATCCCCATTTGTCCggg > 2:380804/1‑150 (MQ=255) ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCCCGCCCGAAAACTAATGTCCGCATGATCATCATCGATTCGCGCTATACCGCAACCCGGGGCGGGCGCGAAGATGAGTGGGGCCCTATTCGTCCggg > 1:1016742/1‑149 (MQ=255) ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg > 2:202732/1‑149 (MQ=255) cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTCGTATCTTGAAACAGCACGCCAGAAATCTAATGCCCGCATGGTCATCATCGATCCGCGCTATCCCCCCCCCCGCGCCGGGCGCGAAAATGAGTGGATCCCCACTCGGCCggg > 1:437577/1‑148 (MQ=255) aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGGTGc < 1:680246/150‑1 (MQ=255) aaTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg < 1:451778/151‑1 (MQ=255) tGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGGCGCACTGGTTaa > 1:302930/1‑150 (MQ=255) ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg < 1:945103/151‑1 (MQ=255) gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc > 2:201098/1‑151 (MQ=255) | TCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC > NC_000913/941571‑941860 |
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: ‑ |