Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,883,651 | C→G | 13.5% | A60G (GCC→GGC) | PP_5157 → | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,883,651 | 0 | C | G | 13.5% | 95.6 / 7.6 | 37 | A60G (GCC→GGC) | PP_5157 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (15/17); new base G (3/2); total (18/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.60e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01 |
CTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTGC > NC_002947/5883506‑5883800 | cTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGGCCTg > 2:642673/1‑150 (MQ=255) aTCTGTCCATTGGTTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGcc > 2:416280/1‑150 (MQ=255) cAGTCCATTGGCTCAGGGGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCGAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCaa < 2:156970/150‑1 (MQ=255) aGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGGCCTGGCCaaa < 1:642673/150‑2 (MQ=255) gTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAgg > 2:1178263/1‑150 (MQ=255) ggCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGa < 2:648510/150‑1 (MQ=255) tCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGAttt > 1:182284/1‑150 (MQ=255) gggcggcGCTAACAGACGAAAACGCCAAGCGCTGCAGACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTgc < 2:12803/149‑1 (MQ=255) gATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATATAAGGTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGtt > 1:865887/1‑150 (MQ=255) gacgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTgc > 2:1076679/1‑136 (MQ=255) gacgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTgc < 1:1076679/136‑1 (MQ=255) gaAATCGCCGACCGGGGCGAACTTCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGaa < 2:58599/150‑1 (MQ=255) ccGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCgg < 2:1392677/150‑1 (MQ=255) ccGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTAc < 1:188213/150‑1 (MQ=255) aTTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGgt > 1:451012/1‑150 (MQ=255) aTTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGgt > 1:1254313/1‑150 (MQ=255) tGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGgtgt < 2:655767/150‑1 (MQ=255) cGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCga > 1:285595/1‑150 (MQ=255) cGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCaa < 1:416280/150‑1 (MQ=255) cGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCaa > 1:164172/1‑150 (MQ=255) gagaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGGGTCGAGAGTGGCAAg > 1:233581/1‑150 (MQ=255) cTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGGCCTGGCCaa > 1:1108098/1‑72 (MQ=255) cTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGGCCTGGCCaa < 2:1108098/72‑1 (MQ=255) gCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACAt > 1:1091313/1‑150 (MQ=255) aCGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCAt < 2:285595/150‑1 (MQ=255) gcacgACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTc > 2:878043/1‑150 (MQ=255) gcaagACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCgcc < 2:1091313/150‑3 (MQ=255) gacgacCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGt < 2:1206069/150‑1 (MQ=255) gacCACCTGCAAGGTTTCGACATCGAATTTGGCCCGGCGCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGaa < 2:233581/150‑1 (MQ=255) ccTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGa < 1:863885/63‑1 (MQ=255) ccTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGa > 2:863885/1‑63 (MQ=255) tCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACg < 2:1254313/150‑1 (MQ=255) aCATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGcc > 2:1158287/1‑150 (MQ=255) tCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGa > 2:1154389/1‑123 (MQ=255) tCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGa < 1:1154389/123‑1 (MQ=255) aTTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGa < 2:182284/150‑1 (MQ=255) ttGGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATg > 1:231997/1‑150 (MQ=255) cccTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTGc > 2:1189154/1‑150 (MQ=255) | CTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTGC > NC_002947/5883506‑5883800 |
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: ‑ |