Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 38,943 | 0 | T | G | 53.1% | 11.8 / 28.6 | 32 | Y58S (TAC→TCC) | caiB | L‑carnitine CoA‑transferase |
Reads supporting (aligned to +/- strand): ref base T (2/13); new base G (17/0); total (19/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.02e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGG‑CCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCC > NZ_CP009273/38854‑39008 | tttCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGGGGGt > 2:267101/1‑90 (MQ=255) cAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAGTTCGGTTGAAc > 1:1779268/1‑90 (MQ=255) cGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAAt < 1:902767/90‑1 (MQ=255) cGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAAt < 1:303323/90‑1 (MQ=255) cgGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCg < 1:1170643/90‑1 (MQ=255) ttCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGGGGGGTGCGGGGAGTTTGGTTGGACGCGAAAGGGGTCGG‑gccc > 1:1494889/1‑89 (MQ=255) tCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGG‑CCCAg < 2:1714403/90‑1 (MQ=255) ttGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGGGAGTTGCGGGGAGTTTGGGTGAACGCGAATGGGGTCGG‑GCCAGgcggcg > 1:521290/1‑90 (MQ=255) ttGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGG‑CCCAGGCGACg < 1:22524/90‑1 (MQ=255) ttGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGG‑CCCAGGCGACg < 1:1597048/90‑1 (MQ=255) ttGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGGAGGTCGGTTGAAAGCGAATGGGGTCGT‑CCCAGgcggcg > 1:1192187/1‑90 (MQ=255) ttGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAGTTCGGTTGAACGCGAATGGGGTTGG‑CCCCGGCgggg > 2:1705145/1‑87 (MQ=255) gAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGGGAGGTTCGGGGAGTTTGGGTGAAAGCGGAAGGGGGTGG‑CCCCGGGGACggt > 2:1597048/1‑88 (MQ=255) aaTATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGG‑CCCAGGCGACGTtctc < 2:1145494/90‑1 (MQ=255) atatTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGG‑CCCAGGCGACGTTCTCg < 2:535397/90‑1 (MQ=255) atTTAACGACAGCGCGTGCAAATTGCGGCGGGGGGGTTGCGGGGAGGTCGGTTGAAACGGAAAGGGGTGGG‑GCCCGGCGGCGGTCTCgct > 2:32955/1‑88 (MQ=255) ttAACGACAGCGCGTGCAAATTGCGGCGGGGGAGGTGGGGGGAGTTTGGGTGGAAGGGGAAGGGGGCGG‑GCCCGGCGGCGGTCTCGGTcc > 2:1112229/1‑90 (MQ=255) ttAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGG‑CCCAGGCGACGTTCTCGATcc < 1:1082349/90‑1 (MQ=255) ttAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGG‑CCCAGGCGACGTTCTCGATcc < 1:32955/90‑1 (MQ=255) gACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAGTTTGGTTGAAAGCGGAAGGGGTCGG‑GCCCGGGGGCGGTCTCGGTTCCGATa > 2:563073/1‑90 (MQ=255) aCAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAGGTCGGGTGAACGCGAAAGGGGTCGG‑CCCAGGGGGCGGTCTCGGTCCAGATaa > 1:1077876/1‑90 (MQ=255) aGCGCGTGCAAATTGCGGCGGGGGAGGTGGGGGTAGTTTGGTTTGACGCGGAAGGGGTTGG‑GCCCGGGGGCGGTCTCGGTCCCGATaaat > 2:230993/1‑88 (MQ=255) aGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGGAGGTTGGTTTAACGCGAAAGGGGGCGG‑CCCAGGCGACGGTCTCGATCCAGATaaat > 1:513057/1‑88 (MQ=255) aGCGCGTGCAAATTGCGGCGGGAGAGGTGGGGGGAGTTTGGGTGAAAGCGAAAGGGGTCGG‑GCCCGGCGGCGGTCTCCGTCCAGAAAACt > 2:446261/1‑90 (MQ=255) gcgcGTGCAAATTGCGGCGGGAGGGGTGCGGGGAGTTTCGGTGAAAGCGAATGGTGTGGGCCCCAGGC‑ACGCTCTCCATCCAGATAACtt > 1:512724/1‑90 (MQ=255) gcgcGTGCAAATTGCGGCGGGAGAGTTGCGGGGAGTTCGGTTGAACGCGAATGGGGGCGG‑CCCAGGCGGCGTTCTCGATCCAGATAACtt > 1:1276085/1‑90 (MQ=255) gcgcGTGCAAATTGCGGCGGGAGAGGTGGGGGGAGGTTGGTTGAAAGCGGATGGGGTGGG‑GCCCGGGGGCGGTCTCGATCCAGATAAAtt > 2:1479271/1‑90 (MQ=255) gcgcGTGCAAATTGCGGCGGGAGAGGTGCGGGGGGGTTGGGTGAAAGCGGAAGGGGTCGG‑GCCCGGCGACGTTCTCGATCCCGATAAAtt > 2:430765/1‑90 (MQ=255) gcgTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGG‑CCCAGGCGACGTTCTCGATCCAGATAACTTcc < 1:1695219/90‑1 (MQ=255) tGCAAATTGCGGCGGGAGAGGTGCGGGGAGTTCGGGTGAACGCGAAAGGGGTCGG‑GCCCGGCGGCGGTCTCGATCCCGATAAATTCcgcc > 2:1793228/1‑89 (MQ=255) aaaTTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGG‑CCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGccc < 1:430765/90‑1 (MQ=255) agTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTc < 1:392549/35‑1 (MQ=255) agTTGCGGGTAGTTCGGGTGAACGCGAATGGTGTc > 2:392549/1‑35 (MQ=255) | TTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGG‑CCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCC > NZ_CP009273/38854‑39008 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |