Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,824,317 | G→C | 12.1% | S2306R (AGC→AGG) | pvdL ← | non‑ribosomal peptide synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,824,317 | 0 | G | C | 12.1% | 117.1 / 6.3 | 42 | S2306R (AGC→AGG) | pvdL | non‑ribosomal peptide synthetase |
Reads supporting (aligned to +/- strand): ref base G (18/18); new base C (3/2); total (22/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 = 9.41e-01 |
ACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGG > NC_002947/4824174‑4824461 | aCCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTggcg > 2:123233/1‑150 (MQ=255) cGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTggcggc > 1:642625/1‑150 (MQ=255) ggCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGt < 2:67713/150‑1 (MQ=255) aGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGCTCATTTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCtt < 2:286210/150‑1 (MQ=255) cATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAg < 2:642625/150‑1 (MQ=255) gttCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGc > 1:404004/3‑150 (MQ=255) cGCCATGGCTCACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGCCCAGCGTTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTgg < 2:332738/150‑1 (MQ=255) caccacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAgg > 1:177324/1‑150 (MQ=255) caccacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAgg < 1:123233/150‑1 (MQ=255) cacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAg > 1:288661/1‑150 (MQ=255) ccTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTTCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggcc < 1:105618/107‑2 (MQ=255) ccTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGc > 1:842410/1‑150 (MQ=255) ccTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggcc > 2:105618/1‑106 (MQ=255) tGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCa < 2:351389/150‑1 (MQ=255) ggTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAgg < 1:125969/150‑1 (MQ=255) gTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggcctaca > 1:57581/1‑100 (MQ=255) gTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggcctaca < 2:57581/105‑6 (MQ=255) tttGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGAgt > 1:570936/1‑150 (MQ=255) cGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCa > 1:598166/1‑128 (MQ=255) cGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCa < 2:598166/128‑1 (MQ=255) tGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACg > 1:808802/1‑150 (MQ=255) tGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACg > 1:301875/1‑150 (MQ=255) tgtAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctggc < 2:177324/150‑1 (MQ=255) gATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGag < 1:742272/150‑1 (MQ=255) gATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGag > 1:129262/1‑150 (MQ=255) aTCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGagc > 1:54265/1‑150 (MQ=255) tCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGagca < 2:301875/150‑1 (MQ=255) gTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCaa > 2:3335/1‑150 (MQ=255) tAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAg > 1:109752/1‑150 (MQ=255) ccTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGa > 2:511751/1‑150 (MQ=255) cGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTc < 1:3335/150‑1 (MQ=255) cAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTc < 1:511751/150‑1 (MQ=255) ttcTGAGGTTGTCCAGCGGTACCTCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCCTAAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGa < 2:180996/148‑1 (MQ=255) ggTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTACAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGATCAGCAAGCCGATGCGGCTGTCTTCGATCATGt > 1:81879/1‑150 (MQ=255) gTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTa < 2:54265/150‑1 (MQ=255) cAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAg > 2:367793/1‑150 (MQ=255) aGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAgg < 2:842410/150‑1 (MQ=255) tGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTc < 2:808802/150‑1 (MQ=255) atgccgccAGTGTGGCGGCATCGTCTTCCAGGCTGCAGCTGGTCAGGCCTTCCGGCAGCTAGCCCAGACTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATACGGCTGGCTTCGATCATATAGCCCAGGCGGTCGAGCGGGTACt > 2:162478/3‑150 (MQ=255) tgccgccAGCCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTc > 1:805644/2‑150 (MQ=255) ccgccAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCCACCCATTCCTTGTGAGCAGCAAGCCGAGGCGGCTGTCTTCGATTATGTAGCGCAAGCGGTCGAGCGGGTACTCgg > 1:444482/1‑150 (MQ=255) gccAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGAGCATGTAGCGCAGGCGGTCGAGCGGTTACTCgggg < 2:288661/150‑1 (MQ=255) | ACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGG > NC_002947/4824174‑4824461 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |