Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 40,515 | 0 | G | T | 57.1% | 23.5 / 39.9 | 28 | Q473K (CAA→AAA) | caiT | L‑carnitine/gamma‑butyrobetaine antiporter |
Reads supporting (aligned to +/- strand): ref base G (2/10); new base T (16/0); total (18/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.03e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGT > NZ_CP009273/40434‑40600 | gTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑GTTGGATT‑gggt > 1:69800/1‑88 (MQ=255) gcgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATT‑GGGTTTAgg > 2:727131/1‑90 (MQ=255) tttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTTATT‑GGGTTTAGggcgcc > 2:763962/1‑90 (MQ=255) tttAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATT‑GGGTTTAGgccccc > 2:576672/1‑90 (MQ=255) gagaGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCcagc < 1:472138/90‑1 (MQ=255) agagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑ATTTTATC‑GGGTTCAGGCCGCCGCGGGccacca > 2:364907/1‑90 (MQ=255) gagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTTATC‑GGGTTCAGGCCGCCCCGCGGcaggcg > 2:266868/1‑87 (MQ=255) gagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATC‑GGGTTTAGGCCGCCCGGCGCccgcag > 2:472138/1‑90 (MQ=255) gagCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCcagcag < 1:83786/90‑1 (MQ=255) gCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTTATCGGGTTTCA‑GCCGCCCCGCGCCCCCAGaa > 1:331876/1‑90 (MQ=255) gCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAAAA‑GGGTTTAGGCCGCCCCGGGCCCGCAGaa > 2:301514/1‑90 (MQ=255) gCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑GTTTAATC‑GGGTTTTGGCCGCCCAGGGGCCGCAGaa > 2:409706/1‑90 (MQ=255) aCCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCCAGCAGAACAATa < 1:284218/90‑1 (MQ=255) aCCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCCAGCAGAACAATa < 2:624556/90‑1 (MQ=255) aaTGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAAAC‑GGGTTTAGGCCCCCGAGCGCCCGCAGAACAAAAACGAt > 2:733458/1‑90 (MQ=255) gTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATC‑GGGTTCAGGCCGCCCGGCGCCAGCAGAACAAAACCGATaaa > 1:321461/1‑89 (MQ=255) tGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATC‑GGGTTCAGGCCGCCGAGCGCCCGCAGAACAAAAACCATAAAGc > 2:683034/1‑90 (MQ=255) tGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGc < 1:409706/90‑1 (MQ=255) tGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGc < 2:178829/90‑1 (MQ=255) gACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTTGAATC‑GGTTT‑AGGCCGCCCAGGGCCCGAAGAAAAAAAACGATAATGcc > 2:389708/1‑90 (MQ=255) gACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATC‑GGTTTTAGGCCGCCCGGGGCCCCCAGAAAAAAAACCATAATGcc > 2:621890/1‑90 (MQ=255) gACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATC‑GGGTTTAGGCCCCCGAGCGCCCGCAGAAAAAAAACCATAAAGcc > 2:169852/1‑90 (MQ=255) gACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTTAATA‑GGGTTTAGGGCGGCGGGCGGCCGCAGAAAAAAAACCATaaaggc > 2:732134/1‑86 (MQ=255) cATCCTCCGGCGATAATGGCGG‑TTTTAATT‑GGTTTTAGGCCGCCGGGGGCCCGCAGAAAAAAAACGAAAATGCCAACCAAAAATTACCaa > 2:193099/1‑90 (MQ=255) tcctccGGCGATAATGGCGG‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAAcc < 1:576672/90‑1 (MQ=255) tcctccGGCGATAATGGCGG‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAAcc < 1:763962/90‑1 (MQ=255) cGGCGATAATGGCGG‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATAc < 1:301514/90‑1 (MQ=255) g‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAggt < 1:352720/90‑1 (MQ=255) | GTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGG‑TTTGAATC‑GGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGT > NZ_CP009273/40434‑40600 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |