Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 955,151 | 0 | G | A | 21.9% | 70.2 / 17.2 | 41 | intergenic (+105/‑208) | cyoE/alaC | protoheme IX farnesyltransferase/aminotransferase |
Reads supporting (aligned to +/- strand): ref base G (18/14); new base A (5/4); total (23/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-02 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTT > NC_002947/955006‑955288 | gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTc > 3:279671/1‑149 (MQ=255) gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTc < 3:331869/149‑1 (MQ=255) cAGAGCGCTGCGGAGGTGGTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCt < 4:279671/149‑1 (MQ=255) gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGt < 3:66721/149‑1 (MQ=255) cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAG‑GAACCTTGTAGGAGC‑gg < 3:206053/149‑1 (MQ=38) gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAG‑GAACCTTGTAGGAGC‑gg > 6:106018/1‑147 (MQ=37) gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAG‑GAACCTTGTAGGAGC‑gg < 5:106018/147‑1 (MQ=37) tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTGTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGTAGGAGC‑GGGTTTa < 8:172660/148‑1 (MQ=255) acgcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTT‑CTGTGCTGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAG‑GAACCTTGTAGGAGc < 6:67899/128‑1 (MQ=21) acgcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTT‑CTGTGCTGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAG‑GAACCTTGTAGGAGc > 5:67899/1‑128 (MQ=21) cTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGTAGGAGG‑GGGTTTACCCGCGAAGAGGgcagca > 5:256819/1‑147 (MQ=255) cTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAAACCCCTACTAAAA‑AGTTACTGTAGGAAC‑GGGTTTTACCGCGAAAAGGCCAGcaa > 7:322822/1‑147 (MQ=255) ctgcCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGaaaaata > 5:328101/1‑148 (MQ=255) gcCTGTTTCACGAAAACCCCGGCCACTGTTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGTAGGAGC‑GGGTTTAccc > 1:212845/1‑122 (MQ=255) gcCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGTAGGAGC‑GGGTTTAccc < 2:212845/121‑1 (MQ=255) gTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGTAgg > 4:238572/1‑104 (MQ=255) gTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGTAgg < 3:238572/104‑1 (MQ=255) cACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAAACCACGACAAAAA‑GGTTACGGTATGAGCGGGGTTTCCCCGCAAAACGGCCACCCAACAAAAATAAAAAAAcagt > 2:119325/1‑146 (MQ=255) aCGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTTAAACCCCTACAAAAA‑CGTTGCTGTAGGAGC‑GGGTTTAACCCCGAAGAGGCCAGCCCACAAAAAAAAAAATACCTGAt > 8:33945/1‑149 (MQ=255) aaaCCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTTAAACCCCAACTACAAGGGTTGGGGTGGGGGG‑GGGTTT‑CACGGGAATGGGGCCGGACAGGAAAATAAAAATACCGGAattt > 8:22775/1‑148 (MQ=255) ccccGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTGCTACAA‑GGTTGCTGTTGGAGC‑GGGTTTTCCCGCGTAGAGGGCAGGACAGAAAAATAAAAATAACTGATTTTaaa > 3:320859/1‑148 (MQ=37) tgCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGTTAAACCCGCTCCTACAA‑GGTTCCTGTAGTAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGa < 5:336775/149‑1 (MQ=255) tgCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGTAGGAGG‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCt < 2:131251/147‑1 (MQ=255) ccggggtTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCAACAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGaa < 3:238479/146‑1 (MQ=255) ggggTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTATTCGCGGGGTAAACCC‑CGCCGAAAA‑AGGTGAGGTTGGGGC‑GGGGGTAACCGCGCACAGGGCCGCCCAGAAAAAAAAAAAAAAAAGTTTTTTAAATAAAAAATCTGtaaaa > 2:318287/1‑148 (MQ=255) ggggTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCATACAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa < 7:33945/148‑1 (MQ=255) gCCGGCCTCTTCGC‑GGGTAAACCCGCTCCAACAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTGAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACa > 8:285139/1‑149 (MQ=255) ttCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATAACTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAAAAGGAAATTATCCTTTACAGGCAtaggg > 8:287727/1‑145 (MQ=255) gc‑gGGTAAACCCGCTCCTACAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt < 6:328101/149‑1 (MQ=255) gggTAAACCCGCTCCTACAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCg > 4:263184/1‑149 (MQ=255) aCCCGCTCCTACAG‑GAACCTTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTa > 3:65325/1‑149 (MQ=37) cccGCTCCTACAG‑GAACCTTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTa > 6:169798/1‑148 (MQ=37) cccGCTCCTACAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt > 3:230040/1‑149 (MQ=255) cccGCTCCTACAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt > 8:284471/1‑149 (MQ=255) ccGATCATACAC‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc < 1:119325/149‑1 (MQ=255) gCTCCTACAG‑GAACCTTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt < 4:65325/149‑1 (MQ=38) gCTCCTACAG‑GAACCTTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct > 5:344520/1‑148 (MQ=38) aCAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg > 4:281308/1‑134 (MQ=255) aCAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg < 3:281308/134‑1 (MQ=255) cAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc < 3:179926/91‑1 (MQ=255) cAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc > 4:179926/1‑91 (MQ=255) | GCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTG‑TGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAA‑GGTTCCTGTAGGAGC‑GGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTT > NC_002947/955006‑955288 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |