| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 997,447 | 0 | T | G | 55.6% | 4.2 / 11.5 | 19 | V62G (GTA→GGA) | elfG | fimbrial protein |
| Reads supporting (aligned to +/- strand): ref base T (5/3); new base G (10/0); total (16/3) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.86e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGGGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGA > NZ_CP009273/997358‑997516 | cGGGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGt < 2:178658/90‑1 (MQ=255) ggCAGTGTCATACTACTCAAGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAg > 2:283060/1‑90 (MQ=255) aTACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGTGTTTAAGaaa > 2:143668/1‑88 (MQ=255) actCAGGGGGATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTTAAGGAAGAGGCAAATACGGCAGGTGAATGTCAAGACaaaaat > 2:394957/1‑87 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGt > 2:508769/1‑74 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGt < 1:508769/74‑1 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAAttt > 1:256030/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGGGTTAAAGAAAAAttt > 2:317166/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAttt > 2:23580/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAatt > 2:605038/1‑90 (MQ=255) ggggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGAc > 1:115361/1‑79 (MQ=255) ggggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGAc < 2:115361/79‑1 (MQ=255) ggggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCa > 1:34317/1‑90 (MQ=255) aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGGGGa > 2:527563/1‑90 (MQ=255) cGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcg < 2:632413/90‑1 (MQ=255) tataTTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAAcgaa > 1:42834/1‑90 (MQ=255) gTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGAAAAATTTTATCAGTGGAACGAATCGAATGat > 1:370305/1‑90 (MQ=255) tGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCg > 1:325462/1‑90 (MQ=255) gaggaAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGGGGAACGAATCGAATGATTATTATTTTTCCTGTGATTGCGAt > 1:546757/1‑90 (MQ=255) aagaagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat < 2:248328/59‑1 (MQ=255) aagaagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat > 1:248328/1‑59 (MQ=255) agaagCAAATACGGCAGGGGGAGTTAAAGACAAATTTTTTCAGTGGAAAGAATCGAATGATTATTTTGTTTCCTGTGATTGCGATAAAGa > 2:552334/1‑90 (MQ=255) | CGGGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGA > NZ_CP009273/997358‑997516 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |