Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,318,440 | 0 | T | G | 57.1% | 6.4 / 13.7 | 14 | V29G (GTG→GGG) | yciO | L‑threonylcarbamoyladenylate synthase |
Reads supporting (aligned to +/- strand): ref base T (1/5); new base G (8/0); total (9/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.00e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCA > NZ_CP009273/1318356‑1318527 | tGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTg > 1:137147/1‑90 (MQ=255) aGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGtt > 1:861126/1‑90 (MQ=255) tttttATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAAc < 1:696159/90‑1 (MQ=255) aTCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAAc > 1:95646/1‑79 (MQ=255) aTCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAAc < 2:95646/79‑1 (MQ=255) gATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTTTCCAACTGATTCCGGGTATgtg > 2:448402/1‑88 (MQ=255) caGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGt > 1:152384/1‑90 (MQ=255) caGCAACGTCTGATCAACCAGGCGGTGGAGATCGGGCGTAAAGGCGGGGGGATTGTTTATCCAACCGTTTCCGGCTATCCGCTCGGCTgg > 1:836410/1‑89 (MQ=255) aCGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAAt > 1:408962/1‑90 (MQ=255) agaTCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGc < 1:847246/90‑1 (MQ=255) tCGTGCGTAAAGGCGGGGGGATTGTTTATCCAACTGGTTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGGGCgga > 1:616832/1‑88 (MQ=255) tGCGTAAAGGCGGGGGGATTGTTTATCCAACTGATTCCGGCTATGGGCTTGGCTGTAAAATTGAAGGTAAAAACGGGATGGGGCGAAttt > 2:726033/1‑90 (MQ=255) tAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCg < 2:523517/90‑1 (MQ=255) ggTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCa < 2:137147/90‑1 (MQ=255) | TGAGCCAGTTTTTTTATATTCATCCTGATAACCCACAGCAACGTCTGATCAACCAGGCGGTGGAGATCGTGCGTAAAGGCGGGGTGATTGTTTATCCAACTGATTCCGGCTATGCGCTCGGCTGTAAAATTGAAGATAAAAACGCGATGGAGCGTATTTGTCGTATTCGCCA > NZ_CP009273/1318356‑1318527 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |