Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,318,440 | 0 | T | G | 62.5% | 1.9 / 11.0 | 16 | V29G (GTG→GGG) | yciO | L‑threonylcarbamoyladenylate synthase |
Reads supporting (aligned to +/- strand): ref base T (2/4); new base G (10/0); total (12/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.24e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCA > NZ_CP009273/1318354‑1318527 | tATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGAt > 1:44533/1‑90 (MQ=255) atTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAAc < 2:243818/83‑1 (MQ=255) atTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAAc > 1:243818/1‑83 (MQ=255) aCCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCg > 2:556542/1‑90 (MQ=255) aCCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGAAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCg > 2:2282/1‑90 (MQ=255) acaGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTg > 1:454220/1‑90 (MQ=255) acaGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTg > 2:483261/1‑90 (MQ=255) caGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGt > 2:602852/1‑90 (MQ=255) aGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTgga > 2:467560/1‑88 (MQ=255) aCGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAAt > 1:785331/1‑90 (MQ=255) aCGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAAt > 2:495775/1‑90 (MQ=255) gTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTg < 1:660785/90‑1 (MQ=255) tGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTTTGCGCTCGGCTGGAAAATTGAAGATAAAAACGCGATGGAGCGTAttt > 1:24759/1‑90 (MQ=255) gTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTc > 1:674499/1‑90 (MQ=255) gggTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGcc < 2:569711/90‑1 (MQ=255) ggTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCa < 2:785331/90‑1 (MQ=255) | TATGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCA > NZ_CP009273/1318354‑1318527 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |