Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,955,854 | 0 | T | G | 60.0% | 11.4 / 18.6 | 20 | S90S (TCA→TCC) | ppdC | prepilin‑type N‑terminal cleavage/methylation domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (3/5); new base G (12/0); total (15/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.61e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTG > NZ_CP009273/2955765‑2955930 | ttGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt < 1:338171/90‑1 (MQ=255) ggTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTa > 1:97012/1‑90 (MQ=255) ggTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGGGATCTCGCCCTCTCTGCCCCCGGGGGAAAcca > 2:153478/1‑88 (MQ=255) ggTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGGCATCTCCCCTTCTCTTCCCCCCGGGGGAACTa > 2:185354/1‑90 (MQ=255) aCCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTa < 2:97012/90‑1 (MQ=255) tAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGGGAAACTAGCGGAACGCTGATGc > 2:88383/1‑90 (MQ=255) aGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCCCCTTCTCTGCCCCCGGGGGAAACTAGCGTAACGCTGATggt > 1:248049/1‑88 (MQ=255) gCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGGGAAACTAGCGGAACGCTGATGCTg > 2:138610/1‑90 (MQ=255) gCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCCGGGGAAACTAGCGTAACGCTGATGCTg > 2:168795/1‑90 (MQ=255) ggCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGac < 1:39333/90‑1 (MQ=255) actGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGGGAAACTAGCGTAACGCTGGTGCTGACACATCCCGccc > 1:241452/1‑89 (MQ=255) aCGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGa < 2:264875/90‑1 (MQ=255) ttCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGGGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTc > 2:117166/1‑90 (MQ=255) ttCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGGGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCCGCGATGTc > 2:41473/1‑90 (MQ=255) gCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCAt < 1:135845/90‑1 (MQ=255) ggcATCTCGCCTTCTCTTCCCCCGGGGGAAACTAGCGGAAACCTGATGCTAACACATCCCCCCCCCGGTGTCTGCATTCGGTTGACCTGc > 2:267248/3‑90 (MQ=255) gTCATCTCGCCTTCTCTGCCCCCGGGGGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGc > 1:52519/1‑90 (MQ=255) gTCATCTCGCCTTCTCTGCCCCCCGGGGAAACTAGCGTAACGCTGATGCTGACACATCCCCCCTGCGATGTCTGCATTCGGTTGACCTGc > 2:159984/1‑90 (MQ=255) ctctGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTg > 2:330609/1‑90 (MQ=255) cGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCAt > 2:16929/1‑55 (MQ=255) cGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCAt < 1:16929/55‑1 (MQ=255) | TTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTG > NZ_CP009273/2955765‑2955930 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |