Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,498,447 | 0 | C | T | 56.0% | 2.9 / 28.3 | 25 | S311S (TCG→TCA) | yhfS | aminotransferase class I/II‑fold pyridoxal phosphate‑dependent enzyme |
Reads supporting (aligned to +/- strand): ref base C (3/8); new base T (14/0); total (17/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.53e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCA‑CCCA‑CCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAG > NZ_CP009273/3498358‑3498500 | ttAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTc < 1:514385/90‑1 (MQ=255) aaTCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTTATATTTTGACTCTGc > 1:541581/1‑90 (MQ=255) aTCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGca < 2:466477/90‑1 (MQ=255) cagaTGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGca‑ccca‑c < 1:422131/87‑1 (MQ=255) acaATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGca‑ccca‑c < 1:585546/90‑1 (MQ=255) aTGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCA‑CCCA‑CCgg > 1:499662/1‑90 (MQ=255) gTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCA‑CCCA‑CCGGGt < 1:455683/90‑1 (MQ=255) gTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCA‑CCCA‑CCGGGt < 2:175661/90‑1 (MQ=255) gTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCA‑CCCA‑CCGGGt < 2:378789/90‑1 (MQ=255) tCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCA‑CCCA‑CCGGGTaa < 1:249265/90‑1 (MQ=255) cTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTTTTTTTTGACTCTGCA‑CCCC‑CCCGGGaag > 2:585546/1‑90 (MQ=255) cTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTTTATTTTGACTCTGCC‑CCCC‑CCcggggag > 2:190960/1‑86 (MQ=255) cTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTTTATTTTGACTCCGCC‑CCCC‑CCGGGTgag > 2:248608/1‑90 (MQ=255) tGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTGACTCTGCA‑CCCC‑CCCGGGaagg > 2:492380/1‑90 (MQ=255) gATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTTTATTTTGACTCTTCACCCCC‑CCGGGAaagg > 2:297729/1‑90 (MQ=255) tgtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTGACTCTGCA‑CCCC‑CCCGggaagggaag > 1:173063/1‑90 (MQ=255) tgtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTGACTCTGCA‑CCCACCCGggaaaggaaa > 1:699107/1‑90 (MQ=255) tgtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGGTTTTTATTTTTACCCCTCA‑CCCACCCGggaaaggaaa > 2:357056/1‑90 (MQ=255) tgtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTTATATTTTGTCTCTTCA‑CCCC‑CCCGggaagggaag > 1:599501/1‑90 (MQ=255) tgaggATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTGACTCTGCC‑CCCA‑CCCgggaagggaaa > 1:626281/4‑89 (MQ=255) gtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCA‑CCCA‑CCGGGTaaggcaagg > 1:59842/1‑90 (MQ=255) tgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTTTATTTTTTCTCTTCC‑CCCC‑CCGGGGaaggcacagg > 2:697298/1‑86 (MQ=255) tgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTTTATTTTTACTCTTCA‑CCCC‑CCCGGGaaaggaaggg > 2:81384/1‑89 (MQ=255) tgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTTGACTCTGCA‑CCCA‑CCGGGTaaagcaaggc > 1:128079/1‑90 (MQ=255) ggCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTCATATTTCGACTCTGCC‑CCCA‑CCGGGGAAGGGAAGGGACCGCGCtttt > 1:364358/1‑90 (MQ=255) ggAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCA‑CCCA‑CCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAg < 2:59842/90‑1 (MQ=255) | TTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCA‑CCCA‑CCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAG > NZ_CP009273/3498358‑3498500 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |