Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,287,719 | 0 | T | A | 20.5% | 80.8 / 11.0 | 39 | intergenic (‑348/‑45) | acs/nrfA | acetyl‑CoA synthetase (AMP‑forming)/cytochrome c552 nitrite reductase |
Reads supporting (aligned to +/- strand): ref base T (24/7); new base A (0/8); total (24/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.05e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAAC > NC_000913/4287589‑4287849 | tAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAATAAAAA‑CCGc < 2:3783/140‑1 (MQ=255) aaCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGcc > 1:949945/1‑140 (MQ=255) ctctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCAt > 2:613402/1‑140 (MQ=255) aaaaGTGTTTTTTTACATGCACTTAAAATTGATAAAAGACAACATTCAAAGTGTGGTTTTTTGTT‑ACACATGGGGGCGAGCAATGTCATGACAGTGTAGGGGCGTTTACTGTGGTAAAAAACAAAAAAAA‑CCGCCATTGc < 2:949945/139‑1 (MQ=255) tgttATTTTACATGCATTAACAATTAATAAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACAAAAAGGGGCGAGCAATGTCATGACAGTGTAGGTGGGTTTACTGTCGTAAAAAACAATAAAAA‑CCGCCATTGcaacaa < 2:853164/137‑1 (MQ=255) tATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATgg > 2:32593/1‑140 (MQ=255) tATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATgg > 2:631932/1‑140 (MQ=255) ttttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGgcg > 2:257252/1‑140 (MQ=255) ttttACATCCACTTACAATTGATTAAAGACAACATTCACAGTGTGTTTTTTTTTA‑ACACATGGGGGGGAGCAATGTCAAGAAAGTGTGGGTGGGGTTACTGTCGGAAAAAACAAAAAAAA‑CCGCCATTGCAACAATGgcg < 2:767354/140‑1 (MQ=255) tGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCg > 2:611121/1‑140 (MQ=255) aCAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAg < 1:495951/140‑1 (MQ=255) aaaGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAc < 1:551266/140‑1 (MQ=255) aaaGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAc < 1:257252/140‑1 (MQ=255) aaGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACa < 2:913779/140‑1 (MQ=255) aaGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACa < 2:99058/140‑1 (MQ=255) caTTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaa > 1:179317/1‑140 (MQ=255) gtgtGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac > 1:549850/1‑140 (MQ=255) tgGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca > 1:23910/1‑140 (MQ=255) tattGTTAACACAGGGGGGGGAGAATTGTAAGAACAGGGTGGGGGGGTTTTCTGTCG‑GAAAAACAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGcc < 2:549850/138‑1 (MQ=255) tt‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCt > 2:526305/1‑140 (MQ=255) t‑ACACATGGGGGAGAGCAATGTCATGACAGGGTGGGTGGGTTTACTGTCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCtt < 2:81993/140‑1 (MQ=255) acacATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTc > 1:490326/1‑140 (MQ=255) aGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTc > 1:21629/1‑140 (MQ=255) aGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTc > 1:535451/1‑140 (MQ=255) tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctt > 1:630267/1‑140 (MQ=255) tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctt > 1:333861/1‑140 (MQ=255) gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttt > 2:385783/1‑140 (MQ=255) cggtgtGGGTGGGTTTACTGTCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt < 1:374278/138‑1 (MQ=255) aggtgtGGGTGGGGTTACTGTCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt < 1:572849/138‑1 (MQ=255) tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc > 1:310392/1‑140 (MQ=255) tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc > 2:325886/1‑140 (MQ=255) gtAGGTGGGTTTATTGTGATAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACt < 2:535451/140‑1 (MQ=255) gtAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACt > 1:22097/1‑140 (MQ=255) gTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTg > 2:12944/1‑140 (MQ=255) gCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt > 2:764845/1‑140 (MQ=255) cGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTc > 1:572564/1‑140 (MQ=255) tttACTGTCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCAc < 1:32593/139‑1 (MQ=255) gTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCAc > 2:101343/1‑140 (MQ=255) gaaaaaaaaaaaaaaAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa < 2:21629/137‑1 (MQ=255) gAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAAc > 1:290623/1‑140 (MQ=255) | TAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAAC > NC_000913/4287589‑4287849 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |