Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,806,632 | 0 | C | T | 57.1% | 8.0 / 26.3 | 28 | S181L (TCG→TTG) | coaBC | bifunctional phosphopantothenoylcysteine decarboxylase/phosphopantothenate‑‑cysteine ligase CoaBC |
Reads supporting (aligned to +/- strand): ref base C (1/11); new base T (16/0); total (17/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.59e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATA > NZ_CP009273/3806553‑3806712 | aGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAAc < 1:558026/90‑1 (MQ=255) aGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAAc < 1:766019/90‑1 (MQ=255) cAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGCCTTTTTTGCCCCGCAACGAc > 1:7616/1‑90 (MQ=255) ggCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGGGGGGGTTTTTTTCCCCTCCACGACCt > 2:463357/1‑90 (MQ=255) ttGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGaa < 1:20626/90‑1 (MQ=255) ttGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGaa < 2:263658/90‑1 (MQ=255) gtggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCCTTTTTTTCCCCTCCACCACCTGaaaa > 2:20626/1‑89 (MQ=255) gtggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCCTTTTTTGCCCGTCAACCACCCGaaaa > 1:558412/1‑89 (MQ=255) gATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGGGGGCATTTTTTTCCCCTCCACCACCCGAAAAAACTg > 2:419830/1‑90 (MQ=255) gATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGGGGGCATTTTTTGCCCCCCCACGACCCGAAACATCTg > 2:568444/1‑90 (MQ=255) gATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGCGGTTTTTTCCCGCCAACGACCCGAAACATATg > 2:691349/1‑90 (MQ=255) gATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGCCTTTTTTGCCCCCCCACGACCCGAAACATCtt > 2:169339/1‑89 (MQ=255) ggTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGGGGGCCTTTTTTTCCCCCCCACCACCCGAAACATCTTAAAAtt > 2:511413/1‑90 (MQ=255) cGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTTCCCCCCCAAGGCCCTAAACAACTGAAAAATATGATTATCGCCGGccc > 2:374645/1‑90 (MQ=255) gATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCg < 1:787902/90‑1 (MQ=255) ttAACCATTGTGGATATGGCGGGAGCGCCTTTTTTTCCCCCCCACCACCCCAAAAAACTTAACATTATTATTTACGCCCGCCcgccgccg > 2:726203/1‑88 (MQ=255) ccATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAAc < 1:70508/90‑1 (MQ=255) ttGTGGATATGGCGGTAGCGCGTTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCCCCGGCCCCACGCGGGAACCGc > 1:358260/1‑90 (MQ=255) ttGTGGATATGGCGGTAGCGCATTTTTTTCCCGTCAACGACCTGAAACATCTTAACATTATGATTACCCCCCGCCCCACGCGTGaaaccc > 1:778634/1‑86 (MQ=255) ttGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAAAATTATGATTTCCGCCCGCCCCACGCGGGAACCGc > 2:169577/1‑90 (MQ=255) ttGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCCACGACCCGAAAAAACTGAAAATTATGATTACCCCCCGCCCCACGCGGGAACCGc > 2:130608/1‑90 (MQ=255) gtgGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAAAATCTGAACATTATGATTAACGCCGGCCCCACGCGGGAACCGccc > 2:25645/1‑88 (MQ=255) gtgGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTc < 1:130608/90‑1 (MQ=255) ggATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGa < 1:511413/90‑1 (MQ=255) tatGGCGGTAGCGCATTTTTTTCCCCCCAACGACCCGAAAAATCTGAAAAATATTATTTCCCCCCGCCCGACGCGGGAACCCCCCGAAcc > 2:60120/1‑90 (MQ=255) tGGCGGTAGCGCATTTTTCGCCCGTCAACGCCCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCgg < 2:326723/90‑1 (MQ=255) tGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCgg < 1:169577/90‑1 (MQ=255) cgcATTTTTCGCCCTTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTtata < 2:778634/90‑1 (MQ=255) | AGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATA > NZ_CP009273/3806553‑3806712 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |