Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,254,414 | 0 | A | C | 52.6% | 13.7 / 23.5 | 19 | V455G (GTG→GGG) | tdcE | 2‑ketobutyrate formate‑lyase/pyruvate formate‑lyase |
Reads supporting (aligned to +/- strand): ref base A (9/0); new base C (0/10); total (9/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATT > NZ_CP009273/3254328‑3254497 | ggCTGTCCTTCTTTTTGTCGTAACCCACCCCGTCTTCCATCAGGGGTGTTTTTTTCGCCCCGCCCTAAATCTTCAGCTTCTC‑TCcccccc < 2:22302/90‑1 (MQ=255) gTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgcc > 1:22302/1‑90 (MQ=255) ccttttGTCGAAACCCCGCCCGTCCCCCACCGGCGGTGCTTTTTCCGGCCCGCCCTAAATCTTCAGCTCCTCGTCCCCCCCGCCGTTAAt < 2:357417/87‑1 (MQ=255) tGTCGTAACCCCGCCCGTCGCCCACCGGGGGTGCTTTTTTCGGCCCGACCTGAATCTTCAGTTTCTCCCCCCCCCCGCCGTTAATTGCGt < 2:520244/90‑1 (MQ=255) aTCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTTAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagcag > 2:553546/1‑90 (MQ=255) gCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcgc > 1:224097/1‑90 (MQ=255) ggTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt > 2:4466/1‑90 (MQ=255) ggTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt < 1:16025/90‑1 (MQ=255) tGCTTTTTCCGCCCCGCCCGGAGTCTTCCGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGc < 2:50411/90‑1 (MQ=255) gCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCa > 1:144302/1‑90 (MQ=255) gTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa > 1:459873/1‑90 (MQ=255) ttttCGGCCCGCCCTGATTCTTCACCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaa < 2:450065/90‑1 (MQ=38) gcccTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt < 1:102919/88‑1 (MQ=255) aCCTAATTCTCCATCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt < 1:297941/90‑1 (MQ=255) cTGATTCTTCAGCTCCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTg < 2:459873/90‑1 (MQ=37) tCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaa > 2:172793/1‑90 (MQ=255) tCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaa > 2:262362/1‑90 (MQ=255) tCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaa > 2:302908/1‑90 (MQ=255) aGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATc < 1:597351/90‑1 (MQ=255) tcGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAtt < 1:550989/90‑1 (MQ=255) | GGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATT > NZ_CP009273/3254328‑3254497 |
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: ‑ |