Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,560,039 | 0 | T | G | 37.0% | 27.9 / 13.0 | 27 | G394G (GGT→GGG) | intZ | CPZ‑55 prophage; putative phage integrase IntZ |
Reads supporting (aligned to +/- strand): ref base T (8/9); new base G (10/0); total (18/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.83e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGC > NC_000913/2559974‑2560102 | tACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGGaa > 1:3137192/1‑68 (MQ=255) ccTGAAGCAACGGGGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACg < 1:2498642/68‑1 (MQ=255) ccTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCCGGGAACg > 1:2690572/1‑68 (MQ=255) ccTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACg > 1:3229295/1‑68 (MQ=255) ccTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGGAACg > 1:612198/1‑68 (MQ=255) aGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAAt < 1:91210/68‑1 (MQ=255) gCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATc > 1:187280/1‑68 (MQ=255) cAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGGAAAGGCaagac > 1:1674093/1‑65 (MQ=255) aCGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCctc < 1:292863/68‑1 (MQ=255) ggCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTa > 1:1563385/1‑68 (MQ=255) ggCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGAAACGGCAATCTTCTa > 1:1310889/1‑68 (MQ=255) aaTGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCgg < 1:529773/68‑1 (MQ=255) gatgCANGATTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGgaac > 1:2032051/1‑67 (MQ=255) gatgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaat > 1:1482409/1‑68 (MQ=255) tgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaatga < 1:3268787/68‑1 (MQ=255) tgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaatga < 1:283434/68‑1 (MQ=255) gCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCCATCCTCTAGAACCCgaatgaa > 1:2472709/1‑68 (MQ=255) aGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGaata < 1:365674/68‑1 (MQ=255) gATTACGCCCGCGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGaataat < 1:1601209/68‑1 (MQ=255) aTTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAAGGAATAAtt > 1:3222922/1‑68 (MQ=255) aTTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCCGAATGAATAAtt > 1:2363229/1‑68 (MQ=255) aTTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAAtt > 1:2320369/1‑68 (MQ=255) aTTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATACTGTCGAACCGGCAGGaataaat > 1:223580/1‑66 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCg > 1:3278795/1‑68 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCg > 1:3107588/1‑68 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCg > 1:841214/1‑68 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTTGAAACGGAATGAAAAATTAAACg > 1:1692218/1‑68 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCg > 1:439325/1‑68 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCg > 1:2833263/1‑68 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCg > 1:1335103/1‑68 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAAACGGAATGAATAATTAAGCg > 1:396124/1‑68 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAAACGGAATGAATAATTAAGCg > 1:780061/1‑68 (MQ=255) cGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCttt > 1:3332338/1‑68 (MQ=255) gCTATTGATTCAATATTGGCAGGCAACGGCAATCCCCTAGAACCGGAATGAATAATTAACCGGCtttt > 1:65261/1‑68 (MQ=255) ttGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGc > 1:262908/1‑68 (MQ=255) ttGATTCAATATTGGCAGGGAACGGCAATCCCCTAGAACCCGAATGAATAAATAAGCCGCTTTTTTGc > 1:1254149/1‑68 (MQ=255) tGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCa > 1:44643/1‑68 (MQ=255) ttCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAg > 1:375756/1‑68 (MQ=255) ttCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCCTTTTTGCAAAg > 1:2416886/1‑68 (MQ=255) atatTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCa < 1:1188055/68‑1 (MQ=255) tatTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCCAAGCGCaa > 1:927426/1‑68 (MQ=255) atTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAAt < 1:394576/68‑1 (MQ=255) cAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTCAAAGCGCAATTTTGc > 1:1245835/1‑68 (MQ=255) | TACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGC > NC_000913/2559974‑2560102 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |