Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,475,282 | 0 | T | G | 55.2% | 29.6 / 35.6 | 29 | T163P (ACC→CCC) | pepA | leucyl aminopeptidase |
Reads supporting (aligned to +/- strand): ref base T (3/10); new base G (16/0); total (19/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.43e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCGAGATCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTC > NZ_CP009273/4475196‑4475364 | gCCGAGATCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCgg < 1:170887/90‑1 (MQ=255) gaTCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGt < 1:87727/90‑1 (MQ=255) tgctgcTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAacac < 2:757393/90‑1 (MQ=255) tgctgcTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAacac < 2:303937/90‑1 (MQ=255) tgctgcTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAacac < 2:161196/90‑1 (MQ=255) ctgctTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCa < 2:619006/90‑1 (MQ=255) cggcggCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCGGGCTGGGAACCCTCTTACGccggg > 1:335097/1‑86 (MQ=255) cggcggCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTTAAAACCATCTTTCGCAGCg > 2:727193/1‑90 (MQ=255) ggcggcAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGGTGAACACCATCTTTCGCCGCgg > 1:351905/1‑90 (MQ=255) cggcAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGgac < 1:96645/90‑1 (MQ=255) tGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCCCGGTGGACACCATCTTACGCCGCGGGCGACGCGGGTCGCTCttt > 1:439664/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGTCGGGACGTTGAAAACCATCTTTCGCAGCGGACGACGCGGGTCGCTCTTTTTTGgtttt > 2:700252/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAacac < 1:410923/47‑1 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAacac > 2:410923/1‑47 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGGGGGACGGTGAAAACCACCTTACCCCGCGGGCGACGCGGGTTGCTCTTGTTTGttttt > 1:943529/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGGGGCGCGGTGAACACCCTCTTTCCCAGCGGGCGACGCGGGTCGCTCTTGTTTGttttt > 1:521693/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGGACGTTGAAAACCCTCTTTCCCCCCGGGCGACGCGGGTCGCTCTTTTTTGttttt > 2:670100/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCGCGTTGAACACCCTCTTACGCAGCGGGCGACGCGGGTCGCTCTTTTTTGttttt > 1:406378/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCCCCGTGAACACCACCTTACCCAGCGGGCGGCGCGGGTCGCTCTTTTTTGttttt > 1:555442/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGGAAACCATCTTACGCAGCGGGCGGCGCGGGTTGCTCTTGTTTGttttt > 2:678138/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGGAAACCACCTTTCCCAGCGGGCGGCCCGGGTTGCTCTTGTTTGttttt > 2:431454/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGAACACCCCCTTACCCAGCGGGCGACGCGGGTCGGCCTTGTTCGTTTTc > 1:454409/1‑90 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGAACACCATCTTTCGCAGCGGACGGCGCGGGTCGCTCTTGTTCGTTTTc > 2:638439/1‑90 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTc > 2:146460/1‑90 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGGTGAAAACCCTCTTACGCAGCGGGCGACGCGGGTTGCTCTTTTTCGttttt > 2:874055/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGGTGAAAACCATCTTACGCAGCGGGCGACGCGGGTCGCTCTTGTTTGttttt > 2:862430/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGCTGGGGCCGCTGGAAATCAAATTTCGCCGCGGACGCCGCGGGGCGGTCTTGTTCCTGTTc > 2:870352/1‑90 (MQ=255) cgcTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGc < 2:619897/90‑1 (MQ=255) cggcggGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGtc < 1:874055/90‑1 (MQ=255) | GCCGAGATCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTC > NZ_CP009273/4475196‑4475364 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |