Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,830,392 | C→T | 13.0% | E281E (GAG→GAA) | pvdL ← | non‑ribosomal peptide synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,830,392 | 0 | C | T | 13.0% | 112.6 / 7.9 | 46 | E281E (GAG→GAA) | pvdL | non‑ribosomal peptide synthetase |
Reads supporting (aligned to +/- strand): ref base C (18/22); new base T (3/3); total (21/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01 |
CAGGCCATAACTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGCGACATCAGCACGCAAGGCACGCCACTGAAG > NC_002947/4830245‑4830539 | caggccaTAACTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTc < 1:1015876/150‑1 (MQ=255) aaCTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGAccc < 2:1000207/150‑1 (MQ=255) aCTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGTCAGCGAGGCTTCGCTggccgg < 1:89825/150‑6 (MQ=255) gaagaaGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCcgctcgct > 1:229064/1‑150 (MQ=255) gaagaaGCGCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCcgctcgct < 1:994606/150‑1 (MQ=255) tctGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACCGGc > 2:497638/1‑150 (MQ=255) tctGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCCGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGc < 2:329416/150‑1 (MQ=255) tGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGa > 2:7446/1‑150 (MQ=255) gAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGaa < 2:229064/150‑1 (MQ=255) aGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGt < 1:7446/150‑1 (MQ=255) cacaGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTc > 1:244187/1‑148 (MQ=255) cacaGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTc < 2:244187/148‑1 (MQ=255) acaGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCggg < 2:890080/150‑1 (MQ=255) tGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATa > 1:175981/1‑131 (MQ=255) tGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATa < 2:175981/131‑1 (MQ=255) aaCTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTg > 1:375087/1‑150 (MQ=255) gCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGccc > 1:379557/1‑149 (MQ=255) ccTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCt < 2:817869/150‑1 (MQ=255) tGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGc > 2:240796/1‑150 (MQ=255) gATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCa < 1:240796/150‑1 (MQ=255) gATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCa < 2:375087/150‑1 (MQ=255) tAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAgc < 2:861462/150‑1 (MQ=255) aGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAgcg > 1:159320/1‑150 (MQ=255) gCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCg > 1:973172/1‑150 (MQ=255) gCCGGAAAACTCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCg < 1:47644/150‑1 (MQ=255) ccGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCgg < 2:629153/150‑1 (MQ=255) gAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGcc > 2:627468/1‑150 (MQ=255) aGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGa < 1:469199/150‑1 (MQ=255) ccACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGt > 2:658927/1‑150 (MQ=255) ccACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGt > 2:558767/1‑150 (MQ=255) aCGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTAc < 2:1063440/150‑1 (MQ=255) gcgcCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCt > 1:872333/1‑150 (MQ=255) ccAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTg > 1:613254/1‑114 (MQ=255) ccAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTg < 2:613254/114‑1 (MQ=255) cAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGc < 2:872333/150‑1 (MQ=255) cAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGc > 2:602062/1‑150 (MQ=255) cAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGCGACATCAg < 1:508381/150‑1 (MQ=255) ccGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAAc < 2:738142/99‑1 (MQ=255) ccGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAAc > 1:738142/1‑99 (MQ=255) gggtcAGCGAGGCTTCGCTGGCCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGa > 2:740043/5‑120 (MQ=255) gggtcAGCGAGGCTTCGCTGGCCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGa < 1:740043/116‑1 (MQ=255) gggtcAGCGAGGCTTCGCTGGCCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGCGACATCAGCACGCAAGGCa > 1:807183/5‑150 (MQ=255) ggCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGCGACATCAGCACGCAAGGCAc > 1:631865/1‑150 (MQ=255) gtcAGCGAGGCTTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATg < 1:1032839/64‑1 (MQ=255) gtcAGCGAGGCTTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATg > 2:1032839/3‑66 (MQ=255) gCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGCGACATCAGCACGCAAGGCACGCCACTGAAg < 2:973172/150‑1 (MQ=255) | CAGGCCATAACTGGCGAAGAAGCTCTGCGGGTCGAAGCCACAGGCCTGGAACTTGTCGGCAAAGGTGGCCAGGCTGTCCTGGCGGATAGGCTCGGAGCCGGAATAGGCCACGCGCCAGCGGCTCAGGTCGAGCCCGGCCAGCGAAGCCTCGCTGACCCGCTCGCTGCACAGGCGATAAGCGAAGTCGGGGCCACCGCTGATGGTGCCACCGTATTCGCTGATGGCCTGCAACCAGCGCAACGGCCGGGCCAGGAAGTACCTCGGCGACATCAGCACGCAAGGCACGCCACTGAAG > NC_002947/4830245‑4830539 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |