Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,449,251 | 0 | A | C | 50.9% | ‑6.6 / 26.4 | 56 | S18R (AGC→CGC) | PP_5556 | bacteriophage lysis protein |
Reads supporting (aligned to +/- strand): ref base A (19/8); new base C (24/4); total (43/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.05e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGG > NC_002947/3449136‑3449383 | aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCG‑CCCGGGcc > 2:15725/1‑135 (MQ=255) aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCCCGCTGGCGGCG‑CCCGGGcc > 2:680005/1‑135 (MQ=255) gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCGGCG‑CCCCGGCCGCCAGGctg > 1:114518/1‑136 (MQ=255) gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCG‑CCCGGGCAGCCAGGgtg > 1:670836/1‑136 (MQ=255) gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCACCGCTGGCGGCG‑CCCGGGCAGCCgggggg > 1:203691/1‑134 (MQ=255) gCAACCGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCGCAGCGCAGGCCGCG‑CCCGGCCAGccacggat > 1:150799/1‑130 (MQ=255) gtgtGGAGCACCGCCAGGCACTCAACACGCGCGCGCGGGTGGTGGTGCTGTGGGCGATTGGGCCGCCGGCGTGGTCGCCGGGGCGATGGGGGGGGCCGCCCGGGATGGGGGCG‑CCCGGGCAGCCTGGGTGTGGCAg < 2:442084/133‑1 (MQ=255) cgggaggaGCCCAAGGCCCTCAACAAGCGCGCGCTCGAGGTGCTGCTGTGCGCGCGGGGGCCACCCGCCGGGGCGCGGTGGCGCTGCGGCGGGCCGCCAGCGCGGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGc < 2:1040136/134‑1 (MQ=255) ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCGGCG‑CCCGCGCAGCCAGGGTGTGGCAGGCCa > 1:461477/1‑136 (MQ=255) gtggcccGCACGGCCCTATACACGCGCGCGCTCGGGGGGCGGCTGTGCGCGCGGGGGCCACCCGCCGGGCCGCCGGGGCGGTGCGGCGGGCCGCCCGCGCGGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCaa < 2:712214/130‑1 (MQ=255) aTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTTGTCGCCGTTGCGCTGCTTCTGGCCGCCCGCGCGGGCGGCG‑CCCCGGCCGCCGGGGGGTGGCAGGCCAAcccc > 1:673688/1‑136 (MQ=255) aTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCG‑CCCGGGCAGCCAGGGTGTGGCAGGCCAAcccc > 1:368557/1‑136 (MQ=255) cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCG‑CCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCgata > 2:302378/1‑136 (MQ=255) cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCG‑CCCGGGCCGCCGGGGTGGGGCCGGCCAACGCCgata > 1:207639/1‑136 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCG‑CCCGGGCCGCCGGGGTGTGGGAGgccgccgccgaaaa > 1:384516/1‑131 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCG‑CCCCGGCAGCCTGGGTGTGGCAGGCAAACCCCTATaa > 2:8123/1‑136 (MQ=255) gCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTAGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgc < 2:97009/136‑1 (MQ=255) tctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCG‑CCCCGGCCGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAg > 2:258268/1‑136 (MQ=255) aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCG‑CCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTc > 2:794372/1‑136 (MQ=255) cAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCG‑CCCGGGCAGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCg > 1:1080182/1‑136 (MQ=255) cgcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCAGGTCCCCCTTCCCCAGCAGCAGACCCCCACCACAATCCGCC‑CCCCCGCAGCCAGGGTGTTGCAGGCCACCCCCCTAAAGGCCCACCGCGCCGATCac > 2:139199/1‑135 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAgg > 1:259042/1‑136 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCG‑CCCGGGCCGCCGGGGTGTGGCCGGCCCACGCCTATAAGGCGCCGCTCGCTGGTCAgg > 1:437474/1‑136 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCAATAAGGCGCAGCTCGCTGATCAgg > 2:1037384/1‑136 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCG‑CCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAgg > 1:430169/1‑136 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCCCGGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACCCCGATAAGGCGCAGCTCGCTGATCAgg > 2:36045/1‑136 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCG‑CCCCGGCAGCCGGGGTGTGGCAGGCCAACGCCAATAAGGCGCAGCTCGCTGATCAgc > 2:22166/1‑135 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCGGGCGGCG‑CCCCGGCAGCCGGGGTGTGGGAGGCCCACGCCAATAAGGCGCCGCTCGCTGATAAgg > 2:78718/1‑136 (MQ=255) gcTAGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCG‑CCCGGGCAGCCTGGCTGTGGCAGGCCACCCCCGAAAAGGAGCAGCGCGCGGATCAgc > 1:60104/1‑135 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCG‑CCCGGGCCGCCGGGGTGTGGCAGGCCAACCCCGATAAGGCGCAGCTCGCTGATCAGGcc > 1:459085/1‑136 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCCCGCAGGCGGCG‑CCCGGGCCCCCTGGGTGTGGCAGGGCAACGCCTATAAGGCGCCCGTCGCTGATCAGccc > 2:404086/1‑136 (MQ=255) gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCAGATCAGGCCgc > 2:995510/1‑135 (MQ=255) cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCG‑CCCGGGCCGCCGCGGTGTGGGAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCCATGagatacg > 1:548939/1‑130 (MQ=255) aGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCG‑CCCGGGCCGCCTGGGGGTGGCAGGCCCCCGCCGATAAGGCGCAGCTCGCTGATCAGGCCGCTGACTTCGGc > 1:397365/1‑136 (MQ=255) gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCG‑CCCGGGCCGCCGGGGTGTGGGAGGCCAACGCCGATAATGCGCAGCGCGCGGATCAGGCCGATGACTTCggggagga > 2:902071/1‑130 (MQ=255) gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCG‑CCCGGGCACCCGGGGTGTGGCCGGCCCACCCCGATAAGGCGCCGCTCGCGGATCAGGCCGAGGGCGTCGGCCAGCa > 2:960625/1‑136 (MQ=255) gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCCCAGGCGGCG‑CCCCGGCATCCAGTGTGTGGCAGGCCAACGCCTAAAAGGCGCAGCTCGCTGATAAGGCCGCAGCCTTCGGCCAGCa > 2:855405/1‑136 (MQ=255) gcTTGGGCCCCCCGCGGGGGCGCGGGGGCGCGGGGGCGGGCCGCCCGCGGTGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt < 2:563235/136‑1 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt > 1:22139/1‑136 (MQ=255) ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCG‑CCCGGGCCGCCGGGGGGTGGCCGGCCCCCCCCTATAAGGCGCAGCTCGCTGATCAGGCCGCAGACGTCGGCCAGCAGTTGGCg > 2:990896/1‑136 (MQ=255) cGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCG‑CCCGGGCCGCCGGGGTGTGGGAGGCCAACCCCGATAAGGCGCAGCTCGCTGATCAGGCCGCTGACGTCGGCAAGCAGTTGGCGGGGAAgg > 1:240347/1‑136 (MQ=255) tCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCG‑CCCGGGCCGCCGGGGTGTGGGAGGCCCACGCCGATAAGGCGCCGCGCGCGGATCAGGCCGCTGACTTCGGCCAGCAGGTGGGGGGGAAGGGTcgcgc > 2:663013/1‑136 (MQ=255) gcggggTGGCGGTGCTGGGGGCCGCGAGGGCGGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc < 2:627784/131‑1 (MQ=255) gcGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc < 1:1013991/135‑1 (MQ=255) gcGCCGTGGCGCGGCTGCTGGCCGCCAGCGCTGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc < 1:1002966/135‑1 (MQ=255) gTGGCGCGGCTGCCGGCCGCCAGCGCGGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACg < 1:351045/136‑1 (MQ=255) ggCGCTGCTGCTGGCCGCCCGCGCGGGCGGCG‑CCCGGGCCGCCGGGGGGTGGGGGGCCCACCCCGATAAGGCGCCGCTCGCGGGGGAGGGCGCTGGCGTCGGGAAGCAGTTGGGGGGGGAGGCTCGCGCGCACggg > 2:917696/1‑135 (MQ=255) ggCGCTGCTGCTGGCCGCCCGCGCGGGCGGCG‑CCCCGGCCGCCGGGGGGTGGCAGGCCAACCCCGATAAAGCCCCCCTCGCGGATAAGGCCCCAGACTTCGGGACACAGGTGGCGGAGAAGGATCGCGCGCCCggg > 2:919423/1‑135 (MQ=255) ggCGCTGCTGCTGGCCGCCAGCCCTGGCGGCGCCCCGGGCCG‑CGGGGGGTGGGCGGGCAACGCCAAAAAGGCGCCGCTCGGGGGTGAGGCCGCGGCCTAGGGGAAGCAGTTGGCGGAGCAGGAGAGCGCGCCCgcg > 2:452185/1‑134 (MQ=255) gcggcgctgctgCCCGCCCGCGGGGCCGGCG‑CCCGGGCAGCCTGGGTGTTGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTc < 1:754005/131‑1 (MQ=255) cgcTGCTGCTGGCCGCCCGCGCTGGCGGCG‑CCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGGTCAGGCCGATGACTTCGGCCAGCAGTTGGGGGTGAAGGGTCGCGCGCGCggggt > 2:907716/1‑133 (MQ=255) ctgGCCGCCCGCGCTGGCGGCG‑CCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCGCGCTGATCAGGCCGCTGACTTCGGCAAGCAGTTGGCGGTGAAGGATCGCGCGCCCGGTGTGGAGCggg > 2:1104532/1‑136 (MQ=255) gggCCGCCTGCGCGGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCTGga < 2:63522/135‑1 (MQ=255) ggCCGCCAGCGCTGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCGGATCAGGCCGATGAGTTCGGCAAGCAGTTGGCGGAGGAGGATCGCGCGCACGGTCTGGAGCGGgag > 1:1073162/1‑136 (MQ=255) ccgcccgcGCGGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa < 1:663013/136‑1 (MQ=255) ccgcccgcGCGGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa < 1:504696/136‑1 (MQ=255) ccgccAGCGCGGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa < 1:1112921/136‑1 (MQ=255) cccAGCGGGGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAgg < 2:642787/135‑1 (MQ=255) | ATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCG‑CCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGG > NC_002947/3449136‑3449383 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |