Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 630,548 | 0 | A | T | 53.3% | 3.1 / 17.5 | 15 | D95E (GAT→GAA) | ybdM | ParB/RepB/Spo0J family partition protein |
Reads supporting (aligned to +/- strand): ref base A (1/6); new base T (8/0); total (9/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.40e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCAC > NZ_CP009273/630464‑630631 | aCCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGGAAAACAACCAAAAATTGCATGGCGTTTTTTTCCGt > 2:120881/1‑90 (MQ=255) atgatgatgAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAAtgggtt < 1:363324/90‑1 (MQ=255) tgatgaGCCTTTACCAATTTCATGGCGGTGAAAACCATCCCCGATTTTCATGGGATTTTTTTTCTTTTGTGTGAGCACAAatggttgggt > 2:576786/1‑90 (MQ=255) tgatgaGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCCTGGCGTTTTTTTCCGTATGTGTGACCACAAaggggtgggg > 1:757806/1‑89 (MQ=255) tgatgaGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCCTGGCATTTTTTTTCGTATGTGTGACCACAAaggggtgggg > 1:200293/1‑89 (MQ=255) tgatgaGCCTTTACCAATTTCATGGCGGCGAAAAAAATCCCCGATTTTCATGGCATTTTTTTTCGCCTGTGTGACCACAAtaggttgggg > 2:363324/1‑89 (MQ=255) gatgaGCCTTTACCAATTTCATGGCGGTGAAAAACAACCACGATTTTCATTGCATTTTTTTTCCTATGTGTGACCACAAagggttgggga > 2:276180/1‑88 (MQ=255) tttACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATc < 1:336806/90‑1 (MQ=255) tttACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATc < 1:576786/90‑1 (MQ=255) tttACCAATTTCATGGCGGGGAAAAACAACCACGAGTTTCCTGGGATTTTTTTCCTTATGGGTGACCACAAAGGGGTGGGTAAAAACATc > 2:136398/1‑90 (MQ=255) ttACCAATTTCATGGCGGTGAAAAACATCCACGGTTTTCCTGGGATTTTTTTTCGTATGTGTGACCACAATGGGGTGGGTAAAACCATCt > 2:669925/1‑90 (MQ=255) ggCGGTGGAAAACATCCCCGATTTTCATGGCATTTTTATCCCTATGGGTGGCCACAATGGGGTGGGGAAAACCCTCTATTTTAATGGGtt > 2:339468/1‑90 (MQ=255) cATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGtt < 1:669925/90‑1 (MQ=255) cGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGtttttttt < 2:757806/90‑1 (MQ=255) aTTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCAc < 1:481548/90‑1 (MQ=255) | ACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCAC > NZ_CP009273/630464‑630631 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |