Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,806,632 | 0 | C | T | 41.4% | 16.1 / 12.5 | 29 | S181L (TCG→TTG) | coaBC | bifunctional phosphopantothenoylcysteine decarboxylase/phosphopantothenate‑‑cysteine ligase CoaBC |
Reads supporting (aligned to +/- strand): ref base C (3/14); new base T (12/0); total (15/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.05e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGT > NZ_CP009273/3806546‑3806708 | cAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGcc < 2:279389/90‑1 (MQ=255) gTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCg > 1:68516/1‑90 (MQ=255) tGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCAGTTTTTGccccg > 1:245577/1‑88 (MQ=255) aGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAAc < 2:145168/90‑1 (MQ=255) gTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACg > 1:32445/1‑90 (MQ=255) tCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGa < 2:32445/90‑1 (MQ=255) ttGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGaa < 1:244798/90‑1 (MQ=255) ttGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGaa < 2:351736/90‑1 (MQ=255) ttGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGaa < 2:245880/90‑1 (MQ=255) gATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCAACGACCTTAAAAATCTg > 1:144183/1‑90 (MQ=255) gATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTg < 1:141674/90‑1 (MQ=255) tatCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCATTTTTTGCCCCTCAACCACCTGAAAAATCTGaa > 2:14419/1‑90 (MQ=255) tCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAAAATCTGAACATTAt > 1:101525/1‑90 (MQ=255) cTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCCTTTTTTTCCCGTCCACCACCCGAAACAACTGAACAATATGa > 2:391046/1‑90 (MQ=255) gggCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGCATTTTTTGCCCGCCCACGACCCGAAAAATCTGGAAATTATGAtt > 2:69294/1‑90 (MQ=255) aaTGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTAccgc > 1:103868/1‑90 (MQ=255) aaTGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGCATTTTTTTCCCCTCCACGACCTGAAAAAACTTAACATTATGATTAacgc > 2:343870/1‑90 (MQ=255) gATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCg < 1:343870/90‑1 (MQ=255) gTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGAtt > 1:218495/1‑73 (MQ=255) gTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGAtt < 2:218495/73‑1 (MQ=255) cATTGTGGATATGGCGGGAGCGCCTTTTTTTCCCCCCAACGACCCGAAAAAACTTAAAATTATGATTACCCCCCGCCCCGCGCGGGaaac > 2:196839/1‑88 (MQ=255) gtgGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTc < 1:24656/90‑1 (MQ=255) ggATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGa > 1:45835/1‑90 (MQ=255) tGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCgg < 2:68516/90‑1 (MQ=255) tGGCGGTAGCGAATTTTTCGCCCGTCAACGTCCTGAACCATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCgg < 2:386486/90‑1 (MQ=255) cGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGc < 2:112850/90‑1 (MQ=255) gTAGCGCATTTTTTGCCCGTCAACGACCTGAAAAATCTTAAAATTATGATTACCCCCCGCCCCACGCGGGAACCGCCCGATCCGGTGCGt > 2:385451/1‑90 (MQ=255) gcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCg > 1:388471/1‑58 (MQ=255) gcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCg < 2:388471/58‑1 (MQ=255) | CAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGT > NZ_CP009273/3806546‑3806708 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |