Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,287,719 | 0 | T | A | 43.6% | 57.8 / 39.8 | 55 | intergenic (‑348/‑45) | acs/nrfA | acetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome |
Reads supporting (aligned to +/- strand): ref base T (25/6); new base A (0/24); total (25/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.42e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01 |
AGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCT > NC_000913/4287575‑4287868 | agGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑c > 2:891732/1‑151 (MQ=255) aGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGcc > 2:488518/1‑151 (MQ=255) tacctctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACGTTGTAGGTGCGGTTACGGTCGTAAAAAACAAAAAAAA‑CCGCCATTg < 2:657732/142‑1 (MQ=255) tAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGcaaca > 1:547687/1‑150 (MQ=255) aaCTCTAAAGTGGTATTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTTACACATGGGGGCGAGCAATGTAATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGcaaca < 2:861967/149‑1 (MQ=255) aaCTCAAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAAt < 2:610462/151‑1 (MQ=255) aaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACACTAAAAA‑CCGCCATTGCAACAATGGCGCa > 1:310494/1‑150 (MQ=255) aGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAAt > 2:344526/1‑151 (MQ=255) ggTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCg > 1:898338/1‑151 (MQ=255) aTTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGAt > 1:921975/1‑151 (MQ=255) tACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGa > 2:657261/1‑149 (MQ=255) tGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAATATCAATAAAAA‑CCGCCATTGCAACAATGGCGCCATTCGGAt > 2:81417/1‑143 (MQ=255) cACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCt > 2:405085/1‑150 (MQ=255) aCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt > 2:81124/1‑151 (MQ=255) aCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt > 1:416449/1‑151 (MQ=255) tACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTAACACATAGGGGCGAGTAATGTAAGGACAGTGTAGGTGCGGTTATTGTCGTAAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg < 1:116861/149‑1 (MQ=255) aaTTGATGAAAGCCAACATTCACAGTGTGATTATTTTTAACACATAGGGGCGAGAAAGTTCATGCCAGTGTAGGTGCGGTTCCTTTCGGGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAg < 2:248982/151‑1 (MQ=255) tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAAAGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa > 2:582170/1‑150 (MQ=255) acaacaTTCACAGTGTGGTTATTGGTTACACATAGGGGCGAGCAATGCCATGACAGTGTAGGTGCGTTTACTGTCGTGAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATaa < 2:603997/150‑1 (MQ=255) aacaTTCGCAGTGTGATTATTTGTTCCAAAAGGGGACGATCAACGTAATGACTGTGTATGACCGGTTATTGGCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacg < 1:217112/151‑1 (MQ=255) tCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcacg < 2:692588/150‑1 (MQ=255) acaGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCg > 1:531357/1‑151 (MQ=255) caGTGTGTTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGTTTACTGTCGTAAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCg < 2:802694/150‑1 (MQ=255) ttttttGTTACAAAGAGGGGGGGGAAATTTCAAGACATTGTAGGTCCGGTTACTGTCGGGAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt < 2:581838/140‑1 (MQ=255) tttGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt > 1:117751/1‑151 (MQ=255) ttACAAATGGGGGCGAGCAATGTAATGACGGGGTGGGTGCGGTTACTGTCGGAAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt < 1:244624/149‑1 (MQ=255) ttACAAATGGGGGCGAGCAAGGTAATGACAGTGTAGGTGGGGTTCTTGTGGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt < 1:922411/149‑1 (MQ=255) acaagggggggAGAAATGTCATGACGGTGTGGGGGCGGTTACTGTCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGa < 2:358193/143‑1 (MQ=255) caTGGGGGCGAGAAATGTCATGACAGTGTGGGGGCGTTTACTGTCGTGAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt < 1:146854/149‑1 (MQ=255) atgggggCGAACAATGTCATGACAGTGTAGGTGGGTTTACTGGCGGAAAAAACAATAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt < 1:637252/141‑1 (MQ=255) aTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCt > 2:47434/1‑151 (MQ=255) ggggggAGAATTTTCATGCCAGTGTGGGTGGGGTTCTGGTGGGAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt < 2:75940/145‑1 (MQ=255) gggaagaaaTGCCATGACGGGTTGGGGGCGTTTACTGTGGGAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt < 2:523561/144‑1 (MQ=255) cAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt > 1:10117/1‑151 (MQ=255) aaTTTCATGACAGTGGAGGTGGTTTTCCTTTCGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttc < 2:547262/145‑1 (MQ=255) aaTTTCATACCAGTGGAGGGGCGTTTACTGTGGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttt < 2:414793/149‑1 (MQ=255) aaTGTCATGACGGTGTGGGGGGGTTTCCTGTCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt < 1:728837/150‑1 (MQ=255) aTGACGGTGTAGGTGGGTTTATTGTCGTAAAAAACAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttt < 1:925508/142‑1 (MQ=255) tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt > 1:612753/1‑150 (MQ=255) tGACAGGGTGGGTGCGTTTCCTGTCGTGAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCa < 1:618468/145‑1 (MQ=255) gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCCCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCC‑TTTTTCTTTTTCACTTCTGtt > 1:758981/1‑151 (MQ=255) aCAGTGTAGGTGGGTTTCCTGCGGTAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGt < 1:725623/150‑1 (MQ=255) cAGTGTAGGTGCGGTTACTGTCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTc < 1:213688/151‑1 (MQ=255) gtgtGGGTGGGTTTACTGTGGAAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt < 1:224574/148‑1 (MQ=255) ggtagggTGCGGTAAGTGTCGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAAGGGCGCAATTCGGATGAAGCCCCTAGGAAAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt < 2:41817/143‑1 (MQ=255) gtgggTGCGGTTACTGTCGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTg < 2:513653/141‑1 (MQ=255) tggggggggTTGCGGTAGTCAAAAAAAAAAAAAG‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAAAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTg < 1:224561/136‑1 (MQ=255) ggtttactgtgtaaaaaaaaaaaaaaaa‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa < 2:178750/137‑1 (MQ=255) ggTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa > 1:798708/1‑150 (MQ=255) ttctggccgaaaaaaaaaaaaaaaa‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAAc < 2:517299/139‑1 (MQ=255) actgtggaaaaaaaaaaaaaaaa‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACgg < 2:179021/141‑1 (MQ=255) cGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTc > 1:671175/1‑151 (MQ=255) tGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTccc > 1:175345/1‑151 (MQ=255) gAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTccc > 2:351337/1‑150 (MQ=255) aaaaaaaaaaaaaa‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgct < 1:145400/143‑1 (MQ=255) aaaaaaaaTAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGc < 2:189514/151‑1 (MQ=255) aaaaCAAAAAAAATCCGCCATTGCAACAATGGCGCATTTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGc < 2:276959/151‑1 (MQ=255) aaTAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAAcc > 1:531298/1‑151 (MQ=255) aTAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCt > 2:351023/1‑151 (MQ=255) | AGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCT > NC_000913/4287575‑4287868 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |