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% | 8.8 / 21.0 | 18 | V62G (GTA→GGA) | elfG | fimbrial protein |
Reads supporting (aligned to +/- strand): ref base T (3/5); new base G (10/0); total (13/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.54e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGGGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTC > NZ_CP009273/997358‑997523 | cGGGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGt < 1:829897/90‑1 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAAttt > 2:105493/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTTAAGAAAAAatt > 1:488478/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAatt > 1:798221/1‑90 (MQ=255) aGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATc < 1:105493/90‑1 (MQ=255) ggggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAACAAAATTTATCa > 1:171916/1‑90 (MQ=255) aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGa > 2:127070/1‑90 (MQ=255) aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGa > 2:751623/1‑90 (MQ=255) aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGa > 2:540320/1‑90 (MQ=255) tataTTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGAc > 1:141585/1‑69 (MQ=255) tataTTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGAc < 2:141585/69‑1 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGGGTTAAAGACAAAATTTTTCAGTGGGAagaatcgaat > 1:416952/1‑90 (MQ=255) cGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGAAAAATTTTATCAGTGGAAAGAATCGAATGa > 1:47682/1‑90 (MQ=255) ccTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTg < 2:488478/90‑1 (MQ=255) aggaAGAAGCAAATACGGCGGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATa < 2:514987/90‑1 (MQ=255) aggaAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATa < 1:127070/90‑1 (MQ=255) aagaagCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAg > 1:260531/1‑90 (MQ=255) agaagCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTTCTGTGATTGCGATAAAGa > 1:163433/1‑90 (MQ=255) aaTACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAAAGAATCGAATGATTATTATGTTTTCTGTGATTGCGATAAAGACAAAGTc > 1:529939/1‑90 (MQ=255) | CGGGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTC > NZ_CP009273/997358‑997523 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |