Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,188,018 | Δ1 bp | 16.0% | intergenic (+79/‑131) | PP_3683 → / → PP_3684 | hypothetical protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,188,018 | 0 | T | . | 16.0% | 136.3 / 23.3 | 50 | intergenic (+79/‑131) | PP_3683/PP_3684 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (18/23); new base . (3/5); total (22/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01 |
TAGCGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTGTTTTTCGCGTG > NC_002947/4187872‑4188167 | tAGCGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAg > 1:991795/1‑150 (MQ=255) cGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCgg < 2:1247071/150‑1 (MQ=255) ggttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGC‑AAGCGGCTTg < 2:458763/150‑1 (MQ=255) ggttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGC‑AAGCGGCTTg < 2:946544/150‑1 (MQ=255) ggttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGC‑AAGCGGCTTg < 2:880653/150‑1 (MQ=255) gttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTgg < 1:1259785/150‑1 (MQ=255) ttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGC‑AAGCGGCtt < 2:428627/147‑1 (MQ=255) ttgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTTAGC‑AAGCGGCtt > 1:428627/1‑147 (MQ=255) ttccgatcTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGGCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGc < 1:1372874/143‑1 (MQ=255) tgATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCAGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGcc < 1:893439/150‑1 (MQ=255) aTCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTACCGGGc < 2:991795/150‑1 (MQ=255) tGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCt < 2:1237504/150‑1 (MQ=255) tGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCt > 2:1064935/1‑150 (MQ=255) gCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTc > 2:1420669/1‑150 (MQ=255) ccGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTcg < 1:436942/150‑1 (MQ=255) cGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGC‑GGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCg > 1:592833/1‑150 (MQ=255) gATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCg < 1:792689/150‑1 (MQ=255) cAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCaa > 2:770156/1‑150 (MQ=255) gACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTgcgcg > 2:573046/1‑150 (MQ=255) gACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTgcgag > 2:573150/1‑148 (MQ=255) aCAAGCCGGTCTAGCATTCGCTTTTTCACCCCGCGGCGGTCCCCCATGGCGCGCCCCTGGCGGGGCCGGCTTGCGTGAACCAAGCCGCTT‑GCTAGGCGGCTTGGCCCTATCGGGCTTTCATCGTCGCGCCTCCGGCCAAGGCCTgcgcgc < 2:570446/150‑1 (MQ=255) ccTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGc < 2:931292/150‑1 (MQ=255) gCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCaaa < 1:1285165/150‑1 (MQ=255) gCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCaaa < 2:921318/150‑1 (MQ=255) agccGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATc < 1:113479/147‑1 (MQ=255) tctcAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTc < 1:432148/150‑1 (MQ=255) tcAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTc > 1:505945/1‑127 (MQ=255) tcAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTc < 2:505945/127‑1 (MQ=255) tcAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGc < 1:360405/134‑1 (MQ=255) tcAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGc > 2:360405/1‑134 (MQ=255) ggCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCAtc > 1:892384/1‑142 (MQ=255) ggCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCAtc < 2:892384/142‑1 (MQ=255) ggCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATa > 1:925442/1‑150 (MQ=255) ggCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATa > 2:1387489/1‑150 (MQ=255) gTTACCTCATGGCGCGCCCCTGGC‑GGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTc < 2:592833/150‑1 (MQ=255) ccTCATGGCGCGCACCTGGCGGTGGCGTGTGTCTTGAACCAAGCCGATT‑GCTAAGCGGTTTCGCCCTCGTGTTGTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCa < 1:1096310/150‑1 (MQ=255) cTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAg > 1:826313/1‑150 (MQ=255) gcCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACgtgt < 2:826313/150‑1 (MQ=255) cccTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCg > 2:231948/1‑150 (MQ=255) cTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCt > 1:1218275/1‑150 (MQ=255) ggCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGAt > 2:1316679/1‑150 (MQ=255) gCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATc < 2:1218275/150‑1 (MQ=255) gCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATc > 1:593431/1‑150 (MQ=255) cGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGcc > 1:868336/1‑150 (MQ=255) gAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTg < 1:231948/150‑1 (MQ=255) ggccAAGCCGCTTAGC‑AAGCGGCTTGGTCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTc > 2:426436/3‑136 (MQ=255) ggccAAGCCGCTTAGC‑AAGCGGCTTGGTCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTc < 1:426436/134‑1 (MQ=255) aCCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTgcg < 1:1420669/150‑1 (MQ=255) cAAGCCGCTTAGC‑AAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCt > 2:855506/1‑150 (MQ=255) aGCCGCTT‑AGCAAGCGGCTTGGTCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTGt > 2:205225/1‑150 (MQ=255) tAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGTTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTGTTTTTCGCGTg < 1:1064935/150‑1 (MQ=255) | TAGCGTTGGTTGATGCCTCATCGGTAGCTGCCGATCACGATGCCAAGGATGAAGACAAGCCGTCCTAGCAGTCGCTTTCTCAGCTCTCGGCGTTACCTCATGGCGCGCCCCTGGCGGGGCCGGCTGGCGTGAACCAAGCCGCTT‑GCTAAGCGGCTTGGCCCTAACGGGCTTTCATCCTCGCGCCTCCGGCCAAGGCCTGCGCGCTCCGCTTGCCAAATCAAAAGCTCATCTCTAGATAGCTCCCAGTACAACGTGTCGCTTCGCCTCCGATCAAGCCTTGCGCTGTTTTTCGCGTG > NC_002947/4187872‑4188167 |
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: ‑ |