Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 333,785 | 0 | G | T | 52.6% | 20.7 / 36.0 | 19 | V433L (GTG→TTG) | fdrA | acyl‑CoA synthetase FdrA |
Reads supporting (aligned to +/- strand): ref base G (0/9); new base T (10/0); total (10/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCT > NZ_CP009273/333701‑333862 | ccGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTc < 2:34682/90‑1 (MQ=255) gACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGGGGTGGGTTTTTTGCTCGGGTTTTGa > 1:634811/1‑90 (MQ=255) gACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGGGGTCGTGGTGGGTTTTTTGGTCGGGTTTggg > 2:13807/1‑89 (MQ=255) gACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGGGGTCGGGGTGGGTTTTTTGGTCGGGTTTTGa > 2:49080/1‑90 (MQ=255) aCCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGAt < 2:343347/90‑1 (MQ=255) aCCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGAt < 1:386632/90‑1 (MQ=255) ccccACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATc < 1:13807/90‑1 (MQ=255) ccccACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATc < 2:133844/90‑1 (MQ=255) accaACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTTTGCTCGGGTTTGGATCGCAt > 1:125392/1‑90 (MQ=255) cATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGGTGGGTTTTTTGGTCGGGTTTGGATTGCATGAAGaatc > 1:458963/1‑87 (MQ=255) aGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTc < 1:595820/90‑1 (MQ=255) aaGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGaagcga < 1:246358/90‑1 (MQ=255) ggCGGTGGTCGTGGTGGGTTTTTTGTTCGGGTTTGGGTCGCAAGAAGAACCGGGGGGGTCCCCCCACGAGACGATCAAAAaaacgtaaag > 2:532004/1‑88 (MQ=255) ggCGGTGGTCGTGATGGGTTTTTTGCTCGGGTTTGGATCGCATGAAGAACCGGGCGGCTCCCCCACCGAGACGAACAAAAaagcgaaaac > 2:154338/1‑88 (MQ=255) ggCGGTGGTCGGGGTGGGTTTTTTTCTCGGGTTTGGATCGCATGAAGATCCGGGCGGGGCCCCCACCGAGACGATCAAAAaaacgaaaac > 2:685620/1‑88 (MQ=255) ggCGGTGATCGTGGTGGATTTTTTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCCCCATCGAGACGATCAAAAaagcgaaaac > 1:320667/1‑88 (MQ=255) ggCGGTGATCGTGATGGATTTTTTGCTCGGGTTTTGGTCGCATGAAGAACCGGTCGGCTCCCCCCTCGAGACGAGCAAAAaagcgaaaga > 2:480636/1‑89 (MQ=255) cGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATcgccgc < 2:266231/90‑1 (MQ=255) gTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCt < 2:634811/90‑1 (MQ=255) | CCGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCT > NZ_CP009273/333701‑333862 |
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: ‑ |