Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,287,719 | 0 | T | A | 37.0% | 58.1 / 30.7 | 46 | intergenic (‑348/‑45) | acs/nrfA | acetyl‑CoA synthetase (AMP‑forming)/cytochrome c552 nitrite reductase |
Reads supporting (aligned to +/- strand): ref base T (25/4); new base A (0/17); total (25/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.42e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGC > NC_000913/4287584‑4287855 | cTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATAAAAGACAACATTCACGGTGTGGTTTTTTGTTACACATGGGGGCGAGCAATGTCATGACGGTGTGGGTGCGGTTA‑CTGTCGTAAAAAACAaaaaaa < 2:95725/140‑1 (MQ=255) tGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATaaaaa > 2:260760/1‑140 (MQ=255) aCTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAAcc < 2:513645/140‑1 (MQ=255) cTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCg > 2:532877/1‑140 (MQ=255) cTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCg > 1:73280/1‑140 (MQ=255) aaCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGcc > 1:136872/1‑140 (MQ=255) gTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATAGCAACAATg > 1:415651/1‑140 (MQ=255) ttttttACATGCATTTACAATTGATAAAAAACAACATCCACTGTGGGTTTTTTTGTTACACATAGGGGCGAGCAATGTCATGCCAGTGAAGGTGGGTTTC‑CTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATgg < 2:98035/138‑1 (MQ=255) tATTTTACATCCACTTAAAATTGATTAAAAACAACATCCAAAGGGTGTTTATTTGTAACTAATAGGGGCGAGAAATGTAATAACAGTGAGGGTGGGGTTA‑TTGTCGTAAAAAACAAAAAAAGCCGCCATTGCAACAATgg < 2:138074/140‑1 (MQ=255) ttttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGgcg < 2:46038/140‑1 (MQ=255) cATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAAtt > 1:112132/1‑140 (MQ=255) cATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAAtt > 1:665251/1‑140 (MQ=255) cATGCACTTAAAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAAtt < 1:390576/140‑1 (MQ=255) aCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGAt > 2:316053/1‑140 (MQ=255) tACAATTTATAAAAGACAATATTCACAGTGGGGTTATTTGTTAAAAATAGGGGCGAGCAATGTAATGACAGTGTAGGTG‑GGTTATCTGCCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGaa < 2:349217/140‑1 (MQ=255) tACAATTGATTAAAGAAAACATTCACAGTGTGTTTTTTTGTTACACATGGGGGCGAGCAATGTCATGACAGTGTAGGGGCGGTTA‑CTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGaa < 2:490467/140‑1 (MQ=255) tACAATTGATAAAAGACAACATTCACAGTGTGTTTTTTTGTTACAGATGGGGGCGAGCAATGTCATGACAGTGTGGGTGCGTTTA‑TTGTCGTAAAAAACAAAAAAAACCGCCATAGCAACAATGGCGCAATTCGGATGaa < 2:415651/140‑1 (MQ=255) aCAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAg > 1:51792/1‑140 (MQ=255) aCAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAg > 1:634953/1‑140 (MQ=255) aCAATTGATAAAAGACAACATTCACAGTGTGTTTATTTTTTACAAATAGGGGCGAGCAAGGTCATGCCAGTGTAGGTGCGGTTA‑CTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAg < 2:51792/140‑1 (MQ=255) aTTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt > 2:309430/1‑140 (MQ=255) ttAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg > 1:490467/1‑140 (MQ=255) aaGACAACATACACAGTGTGGTTATTTGTTACACATGGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACa < 2:610070/140‑1 (MQ=255) gtgGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgc > 2:310272/1‑140 (MQ=255) ggtttATTTGTAACACAAGGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca < 2:687975/137‑1 (MQ=255) ggTTTTTTGTTCCACATGGGGGGGACCAATGTCATGACAGTGTGGGGGGGTTTA‑TGCTGGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcac < 2:266503/140‑1 (MQ=255) tttGTTACAAATGGGGGGGAGCATTTTCATGACGGTGTGGGTGGGGTTA‑CTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt < 1:316053/140‑1 (MQ=255) tGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTAt > 1:448658/1‑140 (MQ=255) aacacAGGGGGGGGAGAAATGTAATAAAGGTGGGGGGGGGTTTA‑TTGTCGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCtt < 2:193472/139‑1 (MQ=255) acacATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTc > 1:553699/1‑140 (MQ=255) acaaATGGGGGCGAAAAATGTCATGACAGTGTAGGTGCGTTTC‑CTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTc < 2:665251/140‑1 (MQ=255) aaatgggggCGAGCATTTTCATGACAGTGTGGGTGCGGTTA‑TTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAg < 1:260760/135‑1 (MQ=255) gggCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt > 2:325373/1‑140 (MQ=255) cGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGCTAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAt > 1:58957/1‑140 (MQ=255) tGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttt > 1:193072/1‑140 (MQ=255) aggtgtAGGTGGGTTTA‑TTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt < 2:741912/138‑1 (MQ=255) aggtgtAGGGGGGGTTA‑TTGTGGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt < 2:48903/138‑1 (MQ=255) aaGTGTGGGGGCGGTTA‑CTGTCGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt < 2:701443/139‑1 (MQ=255) tgtAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc > 1:331324/1‑140 (MQ=255) tgtAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc > 2:434021/1‑140 (MQ=255) ggggTTC‑CTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt < 2:125978/138‑1 (MQ=255) cGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTc > 1:48903/1‑140 (MQ=255) gTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCAc > 2:453661/1‑140 (MQ=255) tA‑CTGTCGTGAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGc < 1:434021/140‑1 (MQ=255) aaCAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgct > 1:387247/1‑140 (MQ=255) cAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgc > 2:325197/1‑140 (MQ=255) cAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgc > 1:467705/1‑140 (MQ=255) | CTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTA‑CTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGC > NC_000913/4287584‑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: ‑ |