Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,254,414 | 0 | A | C | 57.7% | 5.5 / 31.0 | 26 | V455G (GTG→GGG) | tdcE | 2‑ketobutyrate formate‑lyase/pyruvate formate‑lyase |
Reads supporting (aligned to +/- strand): ref base A (11/0); new base C (0/15); total (11/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.29e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATT > NZ_CP009273/3254332‑3254497 | gTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgcc > 1:72317/1‑90 (MQ=255) cccTTTGTCGTACCCCCCCCCGTCGCCCACCGGCGGTGTTTTTTTCGCCCCGACTTGATTCTTCACCTTCCCGCCCCCCCCGCCGTTAAt < 1:766480/88‑1 (MQ=255) cccTTTGTCGTAACCCAGCACGTCTCCCACCAGGGGGGCTTTTTTCGGCCCGCCCTAAATTTTCATCTTCCCGCCCCCCCCGCCGTTAAt < 1:574514/88‑1 (MQ=255) ttGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCg > 2:802806/1‑90 (MQ=255) acccAGCGCGTCGTCCATCGGGGGTGCTTTTTTCGGCCCGCCCTAATTCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGagcag < 2:535597/88‑1 (MQ=255) acccAGCCCGTCGTCCATCGGGGGTGTTGTTTTCGCCCCGCCCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGagcag < 2:766982/88‑1 (MQ=255) aTCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagcag > 1:535597/1‑90 (MQ=255) aGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGtt > 2:370625/1‑90 (MQ=255) gtcgtcCATCAGCGGTGTTTTTTTCGGCCCGCCCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc < 2:186131/90‑1 (MQ=255) gtcgtcCATCAGCGGTGTTGTTTTCGGCCCGCCCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc < 2:72317/90‑1 (MQ=255) gtcgtcCATCAGCGGTGCTTTTTTCGGCCCGCCCTAAATTTTCAGCTTCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc < 2:100832/90‑1 (MQ=255) gtcgccCACCGGGGGTGTTTTTTCCGCCCCCACTTAATTCTTCAGCTCCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc < 1:662971/90‑1 (MQ=255) tcgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc > 1:186131/1‑90 (MQ=255) tcgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc > 2:39978/1‑90 (MQ=255) tcgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc > 2:312554/1‑90 (MQ=255) ggTGTTTTTTTCGCCCCGCCCTAAATTTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt < 2:88634/90‑1 (MQ=255) ggTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt > 1:182353/1‑90 (MQ=255) tGCTGTTTTCGCCCCGACCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGc < 2:18867/90‑1 (MQ=255) gTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa > 2:36288/1‑90 (MQ=255) ttttCGGCCCGTCCAGATTCTTCCCCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaa < 1:370625/90‑1 (MQ=25) cGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa > 2:778866/1‑90 (MQ=255) aaTCTTCTGCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt < 1:443242/90‑1 (MQ=37) aaTCTTCAGCTCCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt < 2:88471/90‑1 (MQ=37) aGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATc < 2:47326/90‑1 (MQ=255) tctcGTCCACACCGCCGTCAATTACGTAGAGCTGCGTGTTCGCCGGGTTAGCGCGTACCCCAAGGAACTGCATTTGCTTACCAACCActa > 1:343679/1‑88 (MQ=255) tcGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAtt < 2:486043/90‑1 (MQ=255) | GTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATT > NZ_CP009273/3254332‑3254497 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |