Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,133,672 | G→T | 7.8% | A114D (GCC→GAC) | PP_3637 ← | sulfonate ABC transporter ATP‑binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,133,672 | 0 | G | T | 7.8% | 157.1 / 3.5 | 52 | A114D (GCC→GAC) | PP_3637 | sulfonate ABC transporter ATP‑binding protein |
Reads supporting (aligned to +/- strand): ref base G (24/23); new base T (2/2); total (27/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 = 6.68e-01 |
CCAGCAACAAGACTTCGGGCTGGCGTGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCAGACGCGGGTGCGGCCCCTGCAGCGGCTGGCCGTGCAAGCTCACGCT > NC_002947/4133523‑4133820 | ccAGCAACAAGACTTCGGGCTGGCGTGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCgg > 1:987631/1‑150 (MQ=255) caaGACTTCGGGCTGGCGTGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCGATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATCGTCCATCAGTACACGACCATCGAGATCCGc > 2:68238/1‑150 (MQ=255) caaGACTTCGGGCTGGCGTGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGc > 2:1206924/1‑150 (MQ=255) aCTTCGGGCTGGCGTGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGc < 2:711735/150‑1 (MQ=255) ttCGGGCTGGCGTGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCt > 1:298002/1‑150 (MQ=255) ggCTGGCGTGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGc < 1:1016152/150‑1 (MQ=255) gCTGGCGTGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCt > 2:714693/1‑150 (MQ=255) ggCGTGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGAc < 1:595885/150‑1 (MQ=255) cGTGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACgg < 2:716261/150‑1 (MQ=255) tGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCt > 2:418994/1‑138 (MQ=255) tGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCt < 1:418994/138‑1 (MQ=255) cAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGaa > 2:729502/1‑150 (MQ=255) aGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGTCATCCGCCAGCCTGGGCTGACGGGCGAAg > 1:98960/1‑150 (MQ=255) aGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGTCATCCGCCAGCCTGGGCTGACGGGCGAAg > 1:99002/1‑150 (MQ=255) gCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGGCGGGCGAAGTCCAg < 1:584953/150‑1 (MQ=255) gggCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTAATGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGc < 2:1092089/150‑1 (MQ=255) cAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGaa > 2:1132618/1‑150 (MQ=255) tGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGNCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGt > 1:598119/1‑150 (MQ=255) ccATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCa > 2:1091596/1‑150 (MQ=255) cATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAg > 2:1153334/1‑150 (MQ=255) cATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAg < 1:729502/150‑1 (MQ=255) aTGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCAGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAgg < 2:203754/150‑1 (MQ=255) gccAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAg > 2:764211/1‑150 (MQ=255) ccAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGc > 2:544034/1‑150 (MQ=255) cAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGcc < 1:826140/150‑1 (MQ=255) aCAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGTCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATg < 2:99002/150‑1 (MQ=255) aCAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGTCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATg < 2:98960/150‑1 (MQ=255) aCAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATg < 2:713113/150‑1 (MQ=255) cAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATgg < 1:1132618/150‑1 (MQ=255) aaCTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGgc > 1:65616/1‑150 (MQ=255) ccGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGc > 1:169708/1‑150 (MQ=255) ccGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGc > 1:846510/1‑150 (MQ=255) gggTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTg < 1:1091596/150‑1 (MQ=255) cTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATc < 2:37007/150‑1 (MQ=255) cACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGaa > 1:262650/1‑150 (MQ=255) cACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGaa > 1:262918/1‑150 (MQ=255) ggCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGc < 1:1160954/150‑1 (MQ=255) gCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGcc < 1:88848/150‑1 (MQ=255) ccAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCa > 1:971003/1‑150 (MQ=255) cAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCAg > 2:1179833/1‑150 (MQ=255) cacaGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCAc > 2:252298/1‑138 (MQ=255) cacaGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCAc < 1:252298/138‑1 (MQ=255) cacaGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCAGACGCgg > 1:1016004/1‑150 (MQ=255) gCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCAGACGCGGGTGc < 2:65616/150‑1 (MQ=255) cTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCAGACGCGGGTGCgg > 1:174909/1‑150 (MQ=255) cGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCAGACGCGGGTGCGGCCCCTg > 2:605296/1‑150 (MQ=255) aTGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTtctc > 1:1011764/1‑81 (MQ=255) aTGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTtctc < 2:1011764/81‑1 (MQ=255) cGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATGCTTAAAAGCCACCGCCAGACGCGGGTGCGGACCATGCAGCGGCt < 2:840175/150‑1 (MQ=255) gATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCAGACGCGGGTGCGGCCCCTGCAGCGGCTg < 1:544034/150‑1 (MQ=255) gTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCAGACGCGGGTGCGGCCCCTGCAGTGGCTGGCCGt < 2:298002/150‑1 (MQ=255) aCGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTGTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCAGACGCGGGTGCGGCCCCTGCAGCGGCTGGCCGTGCaa > 1:92673/1‑150 (MQ=255) aCGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCACGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCAGACGCGGGTGCGGCCCCTGCAGCGGCTGGCCGTGCaa < 1:68238/150‑1 (MQ=255) ggCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGGGAAAGACTGACTCGGGTGCGGCCCCTGCAGCGGCTGGCCGTGCAAGCTCACGCt > 2:937737/1‑150 (MQ=255) | CCAGCAACAAGACTTCGGGCTGGCGTGCCAGGCAGCGGGCCAGCGCGGTGCGCTGGGCCATGCCGCCAGACAACTGCGCCGGGTACTGGCCACGGGCGTGGGCCAGGCCCACAGCCTCAATCGCATGGTCGATGCGTGCACGGCGTTCGGCATCCGCCAGCCTGGGCTGACGGGCGAAGTCCAGGCCGAAGGCGACGTTTTTCTCCAGGCTCAGCCATGGCAGCAGGCTGGGATCCTGAAAAGCCACCGCCAGACGCGGGTGCGGCCCCTGCAGCGGCTGGCCGTGCAAGCTCACGCT > NC_002947/4133523‑4133820 |
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: ‑ |