| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,560,039 | 0 | T | G | 35.1% | 64.1 / 12.2 | 37 | G394G (GGT→GGG) | intZ | CPZ‑55 prophage; putative phage integrase IntZ |
| Reads supporting (aligned to +/- strand): ref base T (11/13); new base G (13/0); total (24/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.69e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATT > NC_000913/2559982‑2560098 | gCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTTACGGCa > 1:2964920/1‑65 (MQ=255)gCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCa > 1:3243664/1‑65 (MQ=255)gCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGGAAAGGCa > 1:167333/1‑65 (MQ=255) ggCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCct > 1:15132/1‑65 (MQ=255) cGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCctct < 1:228416/65‑1 (MQ=255) aaTGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGaag > 1:3417073/1‑64 (MQ=255) aaTGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGCAACGGCAATCCCCTAGAAc > 1:1506635/1‑65 (MQ=255) gttgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAGCGGCAATCCTCTAGAACCgg < 1:2973572/63‑1 (MQ=255) gatgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCgg < 1:1706391/65‑1 (MQ=255) gatgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCgg > 1:1336185/1‑65 (MQ=255) tgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaa > 1:725463/1‑65 (MQ=255) gCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaat > 1:2326178/1‑65 (MQ=255) gCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGgaat > 1:867870/1‑65 (MQ=255) aGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaatga < 1:460143/65‑1 (MQ=255) aGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaatga < 1:578540/65‑1 (MQ=255) aGGATTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGgaatga > 1:2970388/1‑65 (MQ=255) gATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGgaatgaat > 1:2778384/1‑65 (MQ=255) aTTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGaata > 1:2895255/1‑65 (MQ=255) aTTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTCGTACCGGAACGaata > 1:2501621/1‑65 (MQ=255) aTTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGaata > 1:645501/1‑65 (MQ=255) aTTACGCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCCCTAGAACCGGAAAGgata > 1:970191/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTGACGGCAATCCTCTAGAAACGGAATGAataaataa > 1:2458607/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTAACGGGAATCCTCTTGAACCGGAATGAAAAAATaa > 1:1118569/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAGTCCTCTAGAACCGGAAGGAATAATTaa > 1:171293/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTaa > 1:1198970/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTaa > 1:1374890/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTaa > 1:420075/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTaa > 1:3434016/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTaa > 1:3137876/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTaa > 1:2954495/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTaa > 1:322671/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTaa > 1:1039747/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTaa > 1:782582/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTaa > 1:3459780/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTaa > 1:1544504/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAACGGCAATCCTCTAGAAACGGAATGAATAATTaa > 1:2102009/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAAAGGCAATTCTCTAGAACCCGGATGAATAATTaa > 1:1190149/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGGAAAGGCAAGACTCTAGAAGGGGGAAGAATAATTaa > 1:1197171/1‑65 (MQ=255) gCCGACGCTATTGATTCAATATTGGCAGGCTACGGCAACCCTCTTGAACCGCAGGGAATAATTaa > 1:2332264/1‑65 (MQ=255) cGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGc < 1:3430928/65‑1 (MQ=255) ttGATTCAATATTGGCAGGGAACGGCAACCCTCTCGAAACCGAATGAATAAATAAGCGGCttttt > 1:1449320/1‑65 (MQ=255) tGATTCAATATTGGCAGGGAACGGCGATCCTCTAGAACCGGAATGAATAAGTAAGCGGGgttttt > 1:1491372/1‑65 (MQ=255) ttCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCa > 1:1845480/1‑65 (MQ=255) ttCAATATTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCAGCTTTTTTGCa > 1:148212/1‑65 (MQ=255) cAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCaaa < 1:1628429/65‑1 (MQ=255) atatTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgc < 1:1575727/65‑1 (MQ=255) atatTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgc < 1:581120/65‑1 (MQ=255) atatTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgc > 1:1751487/1‑65 (MQ=255) atatTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgc < 1:561324/65‑1 (MQ=255) atatTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgc < 1:54578/65‑1 (MQ=255) atatTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgc < 1:38635/65‑1 (MQ=255) atatTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgc < 1:3097219/65‑1 (MQ=255) atatTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgc < 1:2824267/65‑1 (MQ=255) tatTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgcg > 1:945416/1‑65 (MQ=255) atTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgcgc > 1:3399878/1‑65 (MQ=255) atTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAgcgc > 1:2602890/1‑65 (MQ=255) atTGGCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAACCGGCCTTTTTTCAAAgcgc > 1:553581/1‑65 (MQ=255) ttGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCa < 1:2418861/65‑1 (MQ=255) gtcAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAAt > 1:828758/3‑65 (MQ=255) ggCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAAt < 1:2098459/65‑1 (MQ=255) ggCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAAt < 1:1484488/65‑1 (MQ=255) ggCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTCAAAGCGCAAt > 1:1921694/1‑65 (MQ=255) gacAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAAt < 1:2848205/63‑1 (MQ=255) gCAGGGAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAAtt > 1:2668417/1‑65 (MQ=255) | GCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATT > NC_000913/2559982‑2560098 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |