Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,318,440 | 0 | T | G | 52.9% | 9.3 / 16.0 | 17 | V29G (GTG→GGG) | yciO | L‑threonylcarbamoyladenylate synthase |
Reads supporting (aligned to +/- strand): ref base T (3/5); new base G (9/0); total (12/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.05e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTT > NZ_CP009273/1318357‑1318514 | gAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGt > 1:617405/1‑90 (MQ=255) aGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATc > 1:397616/1‑90 (MQ=255) atatTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGAtt < 1:221229/90‑1 (MQ=255) atTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGGGGAGATCGTGCGTAAAGGCGGGGGGANTTTTTATCCAACTTATTcc > 2:96981/1‑90 (MQ=255) aCGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCg < 1:399511/80‑1 (MQ=255) aCGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCg > 2:399511/1‑80 (MQ=255) aCGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAAt > 2:479652/1‑90 (MQ=255) aCGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGGAAAAt > 2:344595/1‑90 (MQ=255) tCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGa > 1:464318/1‑90 (MQ=255) aaCCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGGAGATaaa > 2:256246/1‑90 (MQ=255) ggAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGa < 1:549987/90‑1 (MQ=255) ggAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGa > 1:63279/1‑90 (MQ=255) agaTCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGc < 1:353769/90‑1 (MQ=255) aTCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGt > 2:642389/1‑90 (MQ=255) gTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTAtt < 2:580237/90‑1 (MQ=255) tGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTAttt > 1:459488/1‑90 (MQ=255) tGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTAttt > 1:271213/1‑90 (MQ=255) | GAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTT > NZ_CP009273/1318357‑1318514 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |