Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,818,324 | A→G | 8.9% | V4304A (GTC→GCC) | pvdL ← | non‑ribosomal peptide synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,818,324 | 0 | A | G | 8.9% | 116.1 / 2.5 | 45 | V4304A (GTC→GCC) | pvdL | non‑ribosomal peptide synthetase |
Reads supporting (aligned to +/- strand): ref base A (20/21); new base G (2/2); total (22/23) | |||||||||||
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.04e-01 |
CAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAGCTGCTTCTGCACCCGCGAAGCGATCTGCGTGG > NC_002947/4818178‑4818462 | cAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCt > 2:559238/1‑150 (MQ=255) aaCGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGtc < 2:284790/150‑1 (MQ=255) aCGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGtca < 2:989516/150‑1 (MQ=255) cGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCatcatc > 1:144738/1‑150 (MQ=255) gACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATcc < 2:52365/150‑1 (MQ=255) tGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAgc < 1:559238/150‑1 (MQ=255) tGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAgc > 1:445172/1‑150 (MQ=255) ccTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCCCTTCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAgcgc > 1:674073/1‑150 (MQ=255) gattgatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACt > 1:636515/1‑150 (MQ=255) gatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGccc > 1:1100010/1‑150 (MQ=255) atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGtg > 2:548311/1‑129 (MQ=255) atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGtc < 1:548311/130‑1 (MQ=255) atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACGGCCCt > 2:1056283/1‑150 (MQ=255) tCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTg > 1:430728/1‑150 (MQ=255) tGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCa < 2:667676/150‑1 (MQ=255) aaTCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGccc < 2:1100010/150‑1 (MQ=255) aTCCGGATGACCGCCTCAGCCGCGGTGTGCGCAGCCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCGCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGcccc < 1:1245418/150‑1 (MQ=255) tCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCt < 2:445172/150‑1 (MQ=255) cGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCATTGCCCTGCAGCCCCTGt < 2:438399/150‑1 (MQ=255) ggATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGGGAAACTCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTa > 1:91482/1‑150 (MQ=255) cGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTc < 2:723496/150‑1 (MQ=255) tgCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATTCACGGCCAACGCCTCCAcc < 2:636515/150‑1 (MQ=255) gACCAACGAAGGGGTCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACt < 2:91482/150‑1 (MQ=255) gAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAg > 2:9045/1‑87 (MQ=255) gAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAg < 1:9045/87‑1 (MQ=255) gAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCCCCCGTTACGCACT‑CAGGTCACTGCGCCCGTCGACTTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACAt < 2:144738/150‑1 (MQ=255) aaaGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATc > 1:431840/1‑150 (MQ=255) aaGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCg > 1:847415/1‑150 (MQ=255) gcgcCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGCCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGc > 2:204421/1‑150 (MQ=255) ccAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGc > 2:55693/1‑150 (MQ=255) cTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCa < 1:55693/150‑1 (MQ=255) gCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGc < 2:491233/150‑1 (MQ=255) cATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGGCCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc < 2:1188676/99‑1 (MQ=255) cATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGGCCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc < 2:1188497/99‑1 (MQ=255) cATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGGCCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc > 1:1188676/1‑99 (MQ=255) cATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGGCCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc > 1:1188497/1‑99 (MQ=255) cacaGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAgc < 2:430728/150‑1 (MQ=255) cacaGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCGCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAgc < 1:1056283/150‑1 (MQ=255) acaGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAgct > 1:1051171/1‑150 (MQ=255) ccccTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAGCTGCtt > 1:853835/1‑150 (MQ=255) cccTCCCGTTTC‑CCATCCAGGGCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAGCTGCTTc < 2:674073/150‑1 (MQ=255) aGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCt < 2:331098/46‑1 (MQ=255) aGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCt > 1:331098/1‑46 (MQ=255) cACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCCCCGTGCTGCACTCGAACATCGCCAGCAGCGGCACGTTCAGCTGTAGCTGCTTCTGCACCCGCGAGGCGACCTTCg > 1:177910/1‑150 (MQ=255) tGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAGCTGCTTCTGCACCCGCGAAGCGATCTGCGTgg > 1:1045947/1‑150 (MQ=255) | CAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAT‑CAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAGCTGCTTCTGCACCCGCGAAGCGATCTGCGTGG > NC_002947/4818178‑4818462 |
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: ‑ |