Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,184,462 | 0 | C | T | 56.2% | 0.5 / 14.3 | 16 | P158S (CCA→TCA) | yqiI | fimbrial protein |
Reads supporting (aligned to +/- strand): ref base C (0/7); new base T (9/0); total (9/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.74e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATAAAGGGAGCCGGAACCTTCGGGAATGGCTACAAAAAGCCCGCAGTCTGTAGTGGCGGTCCTGAACCGCAATTATCAGTAACTTTTCCAGTACGAGTGCAGCTATATATTAAGCTAGCTAAAAATGCCAATAAAGTAAATAAAAAACTTGTATTACCTGACG > NZ_CP009273/3184375‑3184537 | ataAAGGGAGCCGGAACCTTCGGGAATGGCTACAAAAAGCCCGCAGTCTGTAGTGGCGGGCCCGAAACCCAAATATCAATAAATTTTTCa > 2:332475/1‑90 (MQ=255) aaaGGGAGCCGGAACCTTCGGGAATGGCTACAAAAAGCCCGCAGTCTGTAGTGGCGGTCCTGAACCGCAAATATCAATAAATTTTTCAGt > 1:163602/1‑90 (MQ=255) gAGCCGGAACCTTCGGGAATGGCTACAAAAAGCCCGCAGTCTGTAGTGGCGGTCCTGAACCGCAAATATCAGTAAATTTTTCCGGACGAg > 1:273462/1‑90 (MQ=255) aGCCGGAACCTTCGGGAATGGCTACAAAAAGCCCGCAGTCTGTAGTGGCGGTCCTGAACCGCAATTATCAGTAACTTTTCCAGTACGAGt < 1:177062/90‑1 (MQ=255) aGCCGGAACCTTCGGGAATGGCTACAAAAAGCCCGCAGTCTGTAGTGGCGGGCCCGAACCGCAAATATTAATAAATTTTTCCGTACGAGt > 2:329700/1‑90 (MQ=255) ccGGAACCTTCGGGAATGGCTACAAAAAGCCCGCAGTCTGTAGTGGCGGTCCTGAACCGCAATTATCAGTAACTTTTCCAGTACGAGTGc < 2:274287/90‑1 (MQ=255) ccGGAACCTTCGGGAATGGCTACAAAAAGCCCGCAGTCTGTAGTGGCGGTCCTGAACCGCAATTATCAGTAACTTTTCCAGTACGAGTGc < 2:378450/90‑1 (MQ=255) aCCTTCGGGAATGGCTACAAAAAGCCCGCAGTCTGTAGTGGCGGTCCTGAACCGCAATTATCAGTAACTTTTCCAGTACGAGTGCAGCta < 2:110359/90‑1 (MQ=255) ggCTACAAAAAGCCCGCAGTCTGTAGTGGCGGTCCTGAACCGCCATTATTAGTAACTTTTTCAGTACGAGGGCAGCTATATATTAaacta > 1:380720/1‑90 (MQ=255) ggCTACAAAAAGCCCGCAGTCTGTAGTGGCGGTCCTGAACCCCAAATATTAATAAATTTTTCCGTAAGAGGGGAGCCATATATTaaaata > 2:18442/1‑86 (MQ=255) aaGCCCGCAGTCTGTAGTGGCGGTCCTGAAACCCCATTTTTAATAAATTTTTCCATACGAGGGGGGCCATATATTAAATTAACTaaaaaa > 2:353528/1‑89 (MQ=255) ggTCCTGAACCGCAATTATCAGTAACTTTTTCAGTACGAGGGCAGCTATATATTAAACTAGCTAAAAAAGCCAATAAAATAAATaaaaaa > 1:160990/1‑90 (MQ=255) tCCTGAACCGCAATTATCAGTAACTTTTCCAGTACGAGTGCAGCTATATATTAAGCTAGCTAAAAATGCCAATAAAGTAAATAAAAAACt < 1:18442/90‑1 (MQ=255) tCCTGAACCGCAATTATCAGTAAATTTTTCAGTACGAGGGCAGCTATATATTAAACTAGCTAAAAAAACCAATAAAATaaaaaaaaaaat > 2:157845/1‑88 (MQ=255) aaTTATCAGTAACTTTTCCAGTACGAGTGCAGCTATATATTAAGCTAGCTAAAAATGCCAATAAAGTAAATAAAAAACTTGTATTACCTg < 1:210411/90‑1 (MQ=255) tATCAGTAACTTTTCCAGTACGAGTGCAGCTATATATTAAGCTAGCTAAAAATGCCAATAAAGTAAATAAAAAACTTGTATTACCTGACg < 1:19771/90‑1 (MQ=255) | ATAAAGGGAGCCGGAACCTTCGGGAATGGCTACAAAAAGCCCGCAGTCTGTAGTGGCGGTCCTGAACCGCAATTATCAGTAACTTTTCCAGTACGAGTGCAGCTATATATTAAGCTAGCTAAAAATGCCAATAAAGTAAATAAAAAACTTGTATTACCTGACG > NZ_CP009273/3184375‑3184537 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |