Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,287,719 | 0 | T | A | 45.5% | 36.0 / 33.1 | 44 | intergenic (‑348/‑45) | acs/nrfA | acetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome |
Reads supporting (aligned to +/- strand): ref base T (18/6); new base A (0/20); total (18/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.57e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGC > NC_000913/4287570‑4287861 | aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAtc < 2:123535/151‑2 (MQ=255) aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAt > 2:713605/1‑150 (MQ=255) aaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATaa < 1:291448/151‑1 (MQ=255) aaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATaa > 1:879867/1‑150 (MQ=255) gagGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATaaaa > 1:888291/1‑150 (MQ=255) aacttaaaggtgtTATTTTACATGCACTTACAATTGATAAAAGACAACATTAACGGTGTGGTTATTTTTTACACAAAGGGGCGAAAAAGGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAAt < 1:161620/142‑1 (MQ=255) aaCTATAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAAt > 1:128953/1‑151 (MQ=255) ctctAAAGTGTTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACGGTGTGGTTATTTGTTACACATGGGGGCGAGCAATGTCATGACAGTGTGGGGGCGT‑TTACTG‑TCGTAAAAAACAAAAAAAA‑CCGCCATTGCAACAATg < 2:394747/150‑1 (MQ=255) ctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGgcg > 1:439139/1‑151 (MQ=255) gTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAAtt > 1:803802/1‑151 (MQ=255) gTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAAtt > 1:543911/1‑151 (MQ=255) gttATTTTACATGCATTTACAATTGATTAAAGAAAACATTAACAGTGTGGTTTTTTGTTCCCCATAGGGGCAAGCAATGTAATGACAGGGTGGGTGCGG‑TTACTG‑TCGTAAAAAAAAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTc < 1:135019/148‑1 (MQ=255) aTTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGAt > 1:478580/1‑151 (MQ=255) tttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATg > 1:646436/1‑150 (MQ=255) tttACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTTTTTTTTACACATGGGGGCGAGCAATGTCATGACAGTGGAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATg < 1:543863/150‑1 (MQ=255) tgaccTTAAAATTGATTAAAGAAAACATGCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTTTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGccc < 1:101047/146‑1 (MQ=255) gCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGT‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGcccc < 2:882840/150‑1 (MQ=255) aCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt > 2:120177/1‑151 (MQ=255) aCAATTGATAAAAGACAACATTCACAGTGTGGTTATTTTTTACACATAGGGGAGAGCAATGTCATGACAGTGTAGGTGCGT‑TTACTG‑TCGTGAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAc < 2:80705/150‑1 (MQ=255) aTTAATAAAAGACAACATTCACAGGGTGGTTATTTGTTACACATAGGGGCGAGAATTGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAgg < 2:202167/151‑1 (MQ=255) aTTAAAGACAACATTCACAGTGTGGTTATTTTTTACAAAAAGGGGCGAGAAATGTCATGACAGTGAAGGTGGGG‑TTACTG‑TCGTAAAAAAAAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaa < 1:506690/151‑1 (MQ=255) aaaacaacaTTCAAAGTGTGGTTATTTGTTACACATGGGGGCGACCATTGTCATGACAGTGTAGGTGCGT‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAAt < 2:491560/148‑1 (MQ=255) acaacaTTCACAGTGTGGTTATTTTTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATaaa < 2:846760/151‑1 (MQ=255) aacaTTCACAGTGTGTTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac < 2:886179/150‑1 (MQ=255) gttttttttttACAAATGGGGGAGAGCATTTTAATGACAGT‑TAGGGGGGG‑TTTATGTTCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATCCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTa < 1:225834/141‑1 (MQ=255) ttttttGTAACAAAAGGGGGGGAGAAATGTCATGAAAGTGGAGGTGGGG‑TTACTG‑TCGTAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCtt < 2:639802/145‑1 (MQ=255) ttttttGTAACAAAAGGGGGCGAGAAATGTCAGGACGGTGTAGGGGGGT‑TTACTG‑TCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCa < 1:146825/147‑1 (MQ=255) tttttACAAATGGGGGCGAACAATGTCATGACGGTGTGGGTGGGT‑TTACTG‑TCGTAAAAAACAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt < 2:414692/147‑1 (MQ=255) ttACACATGGGGGCGAAAAAGGTCATGACGGTGTAGGTGCGT‑TTACTG‑TCGGAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt < 1:110386/149‑1 (MQ=255) cacaggggggCGAGAAATGTAATGAGGGGGGGGGGGGGT‑TTACGG‑TCGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt < 2:872830/145‑1 (MQ=255) aTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCt > 1:256326/1‑151 (MQ=255) tAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtt > 2:4096/1‑151 (MQ=255) gggggaagaattGTCATG‑CAGTGGTGGGGCGGTTTACTG‑T‑GTAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttt < 1:1467/140‑1 (MQ=255) ggggCGAAAAATTTAAGGACGGTGTGGGGGGGT‑TTCCTG‑TGGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttt < 1:128085/150‑1 (MQ=255) ggggCGAAAAATGTCATGACAGTGTGGGTGCGT‑TTACTG‑TCGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt < 1:346519/145‑1 (MQ=255) ggggCGAAAAATGTAATGACGGTGTGGGGGCGG‑TAACTG‑TGGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttt < 1:197988/151‑1 (MQ=255) aaTGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc > 1:19295/1‑151 (MQ=255) aTGACAGTGT‑GGTGGGGTTTACTG‑CGGAAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt < 1:869496/151‑1 (MQ=255) gtAGGTGCGT‑TTCCGG‑TGGAAAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACg < 2:91265/150‑1 (MQ=255) gtttcgtgtaataaaaaaCAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAAc < 2:519317/135‑1 (MQ=255) aCTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc > 1:294098/1‑151 (MQ=255) tggaaaaaaaaaaaaaaaa‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgct < 1:513511/143‑1 (MQ=255) cGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTc > 1:771293/1‑151 (MQ=255) gTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTcc > 2:819441/1‑151 (MQ=255) gAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCg > 1:765515/1‑151 (MQ=255) gAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCg > 1:693807/1‑151 (MQ=255) aaaaaCAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGc > 2:545210/1‑151 (MQ=255) | AAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTG‑TCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGC > NC_000913/4287570‑4287861 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |