Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,824,317 | G→C | 12.6% | 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.6% | 147.0 / 12.2 | 57 | S2306R (AGC→AGG) | pvdL | non‑ribosomal peptide synthetase |
Reads supporting (aligned to +/- strand): ref base G (23/26); new base C (4/3); total (28/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.00e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01 |
GATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTC > NC_002947/4824170‑4824463 | gATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCt < 1:254172/150‑1 (MQ=255) tCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGCACGCCCCGCCCTCCCGCACGCTGCCGCTCTCCCGCGCCGCACCGCCGCAGGCCGCCCCGCTgg > 1:104702/1‑150 (MQ=255) cccGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTggcg < 2:784661/149‑1 (MQ=255) aCCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGTCGCCAGGCTggcg > 2:601351/1‑150 (MQ=255) cTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATc < 2:465180/150‑1 (MQ=255) tGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCg < 1:1073601/150‑1 (MQ=255) tGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCTGGGTAGCCGCCGGCTCAGTAAAACCAAAATTGTCCAGCGGTGCATCGCTGTAGGCCGACAGGCTGGCGGCATCGTCTTCCAGGc > 2:213451/1‑150 (MQ=21) tGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGc > 2:730115/1‑150 (MQ=255) tGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggcc > 1:469592/1‑136 (MQ=255) tGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggcc < 2:469592/137‑2 (MQ=255) gCCATGGCTGACCACCCCGCCCCGGGGGTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAgacacg < 1:587873/150‑6 (MQ=255) gCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGc < 2:418212/150‑1 (MQ=255) accaccACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCa < 1:455592/150‑1 (MQ=255) accaccACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCa < 1:984648/150‑1 (MQ=255) caccacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAgg < 2:51313/150‑1 (MQ=255) caccacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggcctaca < 1:544957/119‑6 (MQ=255) caccacGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggcctaca > 2:544957/1‑114 (MQ=255) ggTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAgg < 2:716164/150‑1 (MQ=255) gTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAggg < 1:934217/150‑1 (MQ=255) gAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCagcagc < 2:1084299/150‑1 (MQ=255) gAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCagcagc < 1:828945/150‑1 (MQ=255) gttgtAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctg < 1:1035138/148‑1 (MQ=255) gtAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCAGCCAGCGTGCGCTGCTGCACACcctcggc < 2:854174/150‑5 (MQ=255) gATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGag > 1:945318/1‑150 (MQ=255) gATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGag > 1:767054/1‑150 (MQ=255) gATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGCGTGCGCTGCTGCACACcctggg < 2:279492/146‑2 (MQ=255) aTCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGagc < 1:236357/150‑1 (MQ=255) gTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctgg > 2:948864/1‑139 (MQ=255) gTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACgctgg < 1:948864/139‑1 (MQ=255) aGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGcc > 1:984410/1‑150 (MQ=255) ccTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGACTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGa > 2:597813/1‑150 (MQ=255) tGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATg > 2:609970/1‑150 (MQ=255) tGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCACGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATg > 2:606211/1‑150 (MQ=255) ggTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGCAGGCTGCCAGGCTGGCGGCATCGTCTTCAAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGc < 2:146155/150‑1 (MQ=255) gctgcGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGc < 2:104702/150‑1 (MQ=255) ctgcGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTgcag > 1:887560/1‑114 (MQ=255) ctgcGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTgcag < 2:887560/114‑1 (MQ=255) gCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCtt < 1:730115/150‑1 (MQ=255) gAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCAt > 1:536049/1‑150 (MQ=255) ggcctacAGCGATGCACCGCTGTAGGCCGCCAAGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTa > 1:113231/7‑150 (MQ=255) gTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTa < 2:537266/150‑1 (MQ=255) cAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAg > 2:224696/1‑150 (MQ=255) cGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCGCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCg > 2:591886/1‑150 (MQ=255) gTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCCGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGt > 1:342477/1‑150 (MQ=255) gTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGt > 1:123964/1‑150 (MQ=255) tCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGc > 2:864371/1‑150 (MQ=255) gCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCgg < 1:597813/150‑1 (MQ=255) gatgccgccAGCCTGGCGGCCTCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGagcagc < 2:532151/96‑1 (MQ=255) gatgccgccAGCCTGGCGGCCTCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGagcagc > 1:532151/4‑99 (MQ=255) atgccgccAGCCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACt > 1:973411/3‑150 (MQ=255) ggCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTc > 2:109980/1‑150 (MQ=255) ggCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGCCTATCGGGTACTc > 1:227785/1‑150 (MQ=255) cgccAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCggg < 2:945318/150‑1 (MQ=255) gccAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCaa > 1:226959/1‑95 (MQ=255) gccAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCaa < 2:226959/95‑1 (MQ=255) ccAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGGCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGt > 2:264593/1‑150 (MQ=255) cAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTc > 2:374751/1‑150 (MQ=255) | GATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCTGGCCACGCCTTCAGGCAGCTCGCCCAGGGTGTGCAGCAGCGCACGCTGGCTGAGCAGCAAGCCGATGCGGCTGTCTTCGATCATGTAGCGCAGGCGGTCGAGCGGGTACTCGGGGTC > NC_002947/4824170‑4824463 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |