Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,287,719 | 0 | T | A | 25.4% | 127.7 / 21.9 | 55 | intergenic (‑348/‑45) | acs/nrfA | acetyl‑CoA synthetase (AMP‑forming)/cytochrome c552 nitrite reductase |
Reads supporting (aligned to +/- strand): ref base T (39/2); new base A (0/14); total (39/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.76e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGC > NC_000913/4287591‑4287855 | aCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGcc < 2:617383/139‑1 (MQ=255) aCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGcc > 1:617383/1‑139 (MQ=255) aCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCa > 1:383252/1‑140 (MQ=255) ctctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCAt > 1:907367/1‑140 (MQ=255) ctctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCAt > 1:890901/1‑140 (MQ=255) gTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATg > 2:167130/1‑140 (MQ=255) gTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATg > 1:954306/1‑140 (MQ=255) gTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATg > 1:451183/1‑140 (MQ=255) tATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATgg > 1:215975/1‑140 (MQ=255) cATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAAtt > 1:226338/1‑140 (MQ=255) cATGCACTTAAAATTGATAAAAGAAAACATTCAAAGTGTGGTTATTTGTTA‑CAAATGGGGGCGAGAAATGTAATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAAtt < 2:226338/140‑1 (MQ=255) tttACAATAGATAAAAAACAACATCCACAGTGTGGTTATTTGTTA‑CACAGGGGGGCGAGCAATGTCATGACAGTGTGGGTGCGGTTACTGTCGTAAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATg < 2:954306/139‑1 (MQ=255) tttAAAATTGATAAAAGAAAACATCCACAGTGT‑GTTTTTTGTTACCACATGGGGGCGACAAAGTTCATGACGGTGGAGGGGCGTTTACGTTCGGGAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATg < 2:215975/139‑1 (MQ=255) taaaaTTGATAAAAGAAAACATTCACAGTGGGGTTTTTTGTTA‑AAAATAGGGGCGACCAATGTCATAACAGTGTAGGTGCGGTTCCTGTCGTAAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGaa < 2:907367/137‑1 (MQ=255) aCAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAg > 1:1073/1‑140 (MQ=255) aCAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAg < 2:254788/140‑1 (MQ=255) aTTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGccc < 2:383252/140‑1 (MQ=255) gATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTa > 1:510689/1‑140 (MQ=255) ttAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg > 1:437267/1‑140 (MQ=255) gtgtGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac > 2:759861/1‑140 (MQ=255) tgtgTTTTTTTTTTA‑CAAATGGGGGAGACCAATGTCATGACGGT‑TGGGGGGGGTTTTTTTGGGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacg < 2:219979/140‑1 (MQ=255) tgtgGTTTTTTTTTA‑CACAAGGGGGGGAGAAATGTCATGACAGTGTGGGTGGGTTTACTGTCGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacg < 2:361441/140‑1 (MQ=255) ggtgGTTTTTTGTTA‑CACATGGGGGCGAGAAATGTCATGACAGTGTGGGTGGGTTTACTGTCGTAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacg < 2:718905/139‑1 (MQ=255) gtgGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgc > 1:715098/1‑140 (MQ=255) tgGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca > 1:361441/1‑140 (MQ=255) tgGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca > 1:219979/1‑140 (MQ=255) ttttttttA‑CACAGAGGGGCGAGAAATATCAGGATGGTGTGGGGGGTGTTTTTGTCGTAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGcc < 2:510689/138‑1 (MQ=255) tttGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt > 1:813036/1‑140 (MQ=255) ttA‑CACATGGGGGGGAGAAATGTCATGACAGTGTGGGTGCGTTTATTGTGGTAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCt < 2:813036/140‑1 (MQ=255) ttA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATcc > 1:718905/1‑139 (MQ=255) a‑cacATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTc > 2:852434/1‑140 (MQ=255) cacaTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCa > 1:458450/1‑140 (MQ=255) cacaTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCa > 2:166099/1‑140 (MQ=255) aTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCt > 1:144628/1‑140 (MQ=255) aGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTa > 1:669024/1‑140 (MQ=255) ggggCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt > 2:669952/1‑140 (MQ=255) gggCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGcctgtc > 1:182724/1‑135 (MQ=255) gggCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt > 2:541861/1‑140 (MQ=255) cGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAt > 2:741436/1‑140 (MQ=255) cGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAt > 1:136838/1‑140 (MQ=255) aGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTc > 2:911114/1‑140 (MQ=255) aGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTc > 1:547599/1‑140 (MQ=255) aGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTc > 1:781323/1‑140 (MQ=255) gCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTcc > 1:269985/1‑140 (MQ=255) gCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTcc > 1:570326/1‑140 (MQ=255) cATGACAGTGTAGGTGGGGTTACTGGCGTAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttc < 2:136838/140‑1 (MQ=255) cAGTGTGGGGGGGTTTACTGTCGTAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt < 2:523631/140‑1 (MQ=255) cAGTGTGGGGGCGTTTACTGTCGTAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt < 1:142100/140‑1 (MQ=255) aggtgtAGGGGGGGTTACTGTCGAAAAAAAAAAAAAAAC‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt < 2:554069/138‑1 (MQ=255) aaGTGTAGGGGGGTTTCCTGTCGAAAAAAAAAAAAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt < 2:144628/139‑1 (MQ=255) tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc > 1:612581/1‑140 (MQ=255) tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc > 1:800438/1‑140 (MQ=255) gtAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACt > 1:554069/1‑140 (MQ=255) gTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCAc > 1:953971/1‑140 (MQ=255) aCTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCTCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCt > 1:898746/1‑140 (MQ=255) cggaaaaaaaaaaaaaAC‑CCGCCATTGCAACAATGGCTCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACa < 2:898746/136‑1 (MQ=255) aaaaaaaaaaaa‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc < 2:612581/140‑1 (MQ=255) cAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgc > 1:474275/1‑140 (MQ=255) | ACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTA‑CACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAA‑CCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGC > NC_000913/4287591‑4287855 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |