Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,770,559 | G→T | 10.3% | T2918K (ACG→AAG) | pvdD ← | non‑ribosomal peptide synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,770,559 | 0 | G | T | 10.3% | 107.9 / 4.0 | 39 | T2918K (ACG→AAG) | pvdD | non‑ribosomal peptide synthetase |
Reads supporting (aligned to +/- strand): ref base G (17/18); new base T (2/2); total (19/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01 |
CCAGGTAGAGTTCTCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTGGGTTGCCG > NC_002947/4770419‑4770699 | ccAGGTAGAGTTCTCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTaa > 1:755705/1‑150 (MQ=255) cAGGTAGAGTTCCCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAAc < 1:272117/150‑1 (MQ=255) tAGAGTTCCCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACg < 1:457455/150‑1 (MQ=255) agagTTCTCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACgg < 1:286362/150‑1 (MQ=255) gTTCTCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGttt > 1:612103/1‑150 (MQ=255) acCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATa > 1:152526/1‑150 (MQ=255) gcggcAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCACCGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCAt < 1:292023/150‑1 (MQ=255) gcggcAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCAt < 2:363874/150‑1 (MQ=255) tGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTgatgat < 2:152526/150‑1 (MQ=255) aGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTtcg < 2:431371/120‑1 (MQ=255) aGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTtcg > 1:431371/1‑120 (MQ=255) tCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGc > 1:1121786/1‑150 (MQ=255) aaCGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGTTGAGCCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTg > 1:618545/1‑150 (MQ=255) ccATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGt < 1:1038640/150‑1 (MQ=255) cATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTa > 1:1002770/1‑150 (MQ=255) cATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTa < 2:174604/150‑1 (MQ=255) aTAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATg > 1:1120433/1‑150 (MQ=255) aGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATgcg > 1:843202/1‑150 (MQ=255) aCGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTc < 1:92198/150‑1 (MQ=255) gCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCAt > 1:1119134/1‑150 (MQ=255) ccgCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTGGAAAGAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATa < 1:248584/148‑1 (MQ=255) ccTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTg > 1:452383/1‑150 (MQ=255) aaCGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGccc > 1:1109709/1‑89 (MQ=255) aaCGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGccc < 2:1109709/89‑1 (MQ=255) aaCGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGc > 2:463633/1‑150 (MQ=255) gcgcATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCAc < 2:755705/150‑1 (MQ=255) gcgcATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCAc < 2:1119134/150‑1 (MQ=255) gACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCAAGATCGCCGGGGCCACCACGATTTCGGTCGGCCCGTAGCCATTGATGATGTACTGCGGCTGCGAGTAGGCATGCCCCTGTGCATAGCTGGCCTGCGGCACCGCATCCCCGCCCaa > 1:615876/1‑150 (MQ=255) aCCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAg < 2:1002770/150‑1 (MQ=255) gcgcCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCTTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCaa > 1:1032279/1‑150 (MQ=255) gcgcCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCTTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCaa > 1:1033753/1‑150 (MQ=255) gcAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGAt < 2:1120433/150‑1 (MQ=255) ggCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGc < 1:751165/84‑1 (MQ=255) ggCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGc > 2:751165/1‑84 (MQ=255) ggCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGATgcgc < 1:463633/150‑1 (MQ=255) gcgGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCATCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTggg > 2:1044544/1‑150 (MQ=255) tttCCAGATCTTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTGGGTTGcc < 2:1033753/150‑1 (MQ=255) tttCCAGATCTTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTGGGTTGcc < 2:1032279/150‑1 (MQ=255) ttCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCGCCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTGGGTTGCCg > 2:987584/1‑150 (MQ=255) | CCAGGTAGAGTTCTCCCGCCACACCGGGCGGCAGTGGGTTCAGGTCAACGTCCATGACATAGGCGCTACGCTGCCCTACCAACGTACCTATTGGCGCATAGACCGCGCCGCAAGGCTGCTGTGGCGCGGCTTTCCAGATCGTCGGGGTAACCACGGTTTCGGTCGGCCCATAGCCATTGATGATGTACTGCGGCTGCAAGTAGGCATGCGCCTGTTCATAGCTGGCCTGCGGCACCGCATCACCGCCGAAGCAATAGATGCGCATCGGCGGTGGGTTGCCG > NC_002947/4770419‑4770699 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |