| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,287,710 | 0 | G | A | 38.5% | 44.6 / 13.9 | 39 | intergenic (‑339/‑54) | acs/nrfA | acetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome |
| Reads supporting (aligned to +/- strand): ref base G (16/8); new base A (0/15); total (16/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.52e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCT > NC_000913/4287562‑4287853 | gcTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGa < 2:672784/150‑1 (MQ=255) gTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGaaa < 1:324327/149‑1 (MQ=255) aaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATa > 1:135269/1‑150 (MQ=255) agGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAAc < 2:309052/151‑1 (MQ=255) ggAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCTCTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAAcc < 2:293108/151‑1 (MQ=255) aaGATACTGACTAACTCTAAAGTGGTATTTTACATGAACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAAc > 1:420669/1‑148 (MQ=255) aGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGcc < 1:338220/151‑1 (MQ=255) taaCTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATAAAAAACAACATTCACAG‑GTGGTTTATTTTTTACACAAAGGGGCAAGCAATGTAATGACAGGTGAGGGGGTTTTACGGGCGGAAAAAACAAAAAAAACCGCCATTGc < 1:127635/148‑1 (MQ=255) aCTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGcaac > 1:60582/1‑151 (MQ=255) aaCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTTTGG‑TTATTTGTTACACATAGGGGCGAGCAATTTCATGACAGTGTAGGTGCGTTTACTGTCGTGAAAAACAATAAAAACCGCCATTGcaaca < 2:227394/149‑1 (MQ=255) aaaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGgcgc > 1:103854/1‑150 (MQ=255) tGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTc > 2:116320/1‑151 (MQ=255) ttttACAACCACTTAAAATGGATAAAAGACAACATCCACAGTGTGG‑TTATTTGTTACACAGAGGGGAAAGCAATGTCAGGACCGTGGAGGTGCGGTTACTGTCATGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCgg < 1:331964/148‑1 (MQ=255) cATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGcc > 1:447406/1‑151 (MQ=255) cACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCt > 2:443234/1‑150 (MQ=255) cTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt < 2:294202/150‑1 (MQ=255) ttACAATTAATTAAAGCCAACATTCACAGTGTGG‑TTTTTTTTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGTTTCCTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt < 2:379898/149‑1 (MQ=255) tACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAc > 2:102794/1‑151 (MQ=255) aCAATTGATAAAAAACAACATTCAAAGGGTGT‑TTATTTGTTAAACAAAGGGGCGAAAAATTTAATGACATTGTAGGTGCTTTTTGTGTCGTAAAAACCAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg < 1:485321/148‑1 (MQ=255) tGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGa < 2:243780/150‑1 (MQ=255) ttAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaa > 2:754/1‑150 (MQ=255) aacaTTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac < 2:644970/150‑1 (MQ=255) acaGTGTGG‑TTATTTGTTACAAATAGGGGCGAGAAATGTCATGACAGTGTAGTTGGGTTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCg < 2:634756/151‑1 (MQ=255) tgG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCtt > 2:159364/1‑151 (MQ=255) ggggctATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAAGAATAAAAACCGCCATTGCAACAATGGCGCAACTCGGATGAAGCCCATATGACAAGGATAAAAATAAACCCACGACGTATCt > 1:576509/6‑151 (MQ=255) caaaTAGGGGCGAGAAATGTCATGACGGTGTGGGTGGGGTTACTGGCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg < 2:431119/145‑1 (MQ=255) tcggggCGAGCAATGTCATGACAGTGTGGGTGCGGTTAATGTCGTGATAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAGGCACCTATGACAAGGATAAAAATAAACGCACGCCGTCTCTTCAGCTTATTGATTcc > 1:269691/3‑149 (MQ=255) gtcGAGAAATGTCATGACAGTGAGGGTGGGGTTACTTTGTTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttttgt < 1:568026/149‑3 (MQ=255) gCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttc > 1:533273/1‑151 (MQ=255) cGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct > 1:189073/1‑151 (MQ=255) aatgtaaaaaCAGTGTAGGTGCGTTTACTTTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt < 1:285358/141‑1 (MQ=255) aatgaaaTGACAGTTTGGGGGCGTTTCCTGTGTTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc < 1:517441/145‑1 (MQ=255) aaTGTCATGACGGGGTAGGTGCGTTTCCTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttt < 2:509584/148‑1 (MQ=255) cATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTc > 1:622884/1‑151 (MQ=255) gCGTTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACa < 2:242096/151‑1 (MQ=255) gTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAAc > 2:266183/1‑151 (MQ=255) ttACTGTGGGAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaaa < 2:135105/149‑1 (MQ=255) ttACTGTCGAAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaaa < 2:415989/149‑1 (MQ=255) gTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgct < 2:666183/149‑1 (MQ=255) | GCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGG‑TTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCT > NC_000913/4287562‑4287853 |
| 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: ‑ |