Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,777,456:1 | +T | 8.9% | intergenic (+19/+24) | betA‑II → / ← PP_5065 | choline dehydrogenase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,777,456 | 1 | . | T | 8.9% | 203.0 / 15.7 | 56 | intergenic (+19/+24) | betA‑II/PP_5065 | choline dehydrogenase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base . (28/23); new base T (3/2); total (31/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 = 7.99e-01 |
GAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCC > NC_002947/5777318‑5777596 | gagaAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGcccc > 1:679898/1‑150 (MQ=255) gaTTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCa > 2:1080097/1‑150 (MQ=255) cGGACAAGTTCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGttt < 2:335606/150‑1 (MQ=255) gACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGc > 2:1048349/1‑150 (MQ=255) gACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGc > 1:382480/1‑150 (MQ=255) aCAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCt > 1:4686/1‑150 (MQ=255) gATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAgg < 1:1080097/150‑1 (MQ=255) gATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAgg < 1:995280/150‑1 (MQ=255) cGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCggg < 2:648304/150‑1 (MQ=255) gTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCTGC‑GCGGCCCCtga < 2:380847/131‑3 (MQ=255) gTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑TGCGCGGCCCCtga > 1:380847/1‑129 (MQ=255) gTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGc > 1:407746/1‑150 (MQ=255) tGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCa < 1:949234/150‑1 (MQ=255) tGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCa > 1:868546/1‑150 (MQ=255) ggTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAg < 2:382480/150‑1 (MQ=255) gTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAgc < 1:1017084/150‑1 (MQ=255) gCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCa > 2:115639/1‑150 (MQ=255) ccACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCaaaaaa < 2:868546/150‑1 (MQ=255) cTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtgt < 2:419016/150‑1 (MQ=255) cTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtgt < 2:662577/150‑1 (MQ=255) tGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtgtg < 2:679898/150‑1 (MQ=255) cgcgCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTgcgc > 2:348038/1‑150 (MQ=255) cGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTc > 1:106597/1‑114 (MQ=255) cGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTc < 2:106597/114‑1 (MQ=255) tactacGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTc < 1:981039/150‑1 (MQ=255) cGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGacac < 2:452586/150‑1 (MQ=255) ggggATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTgc > 1:562941/1‑119 (MQ=255) ggggATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTgc < 2:562941/119‑1 (MQ=255) ggggATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCg > 1:698852/1‑150 (MQ=255) gCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGc > 1:956086/1‑150 (MQ=255) cACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGcc > 1:1034796/1‑150 (MQ=255) aaaGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTg > 1:331326/1‑150 (MQ=255) ccGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGaa < 1:357264/150‑1 (MQ=255) gATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGccc < 1:912925/54‑1 (MQ=255) gATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGccc > 2:912925/1‑54 (MQ=255) aTGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGaacaa > 1:929364/1‑150 (MQ=255) aTGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGaacaa < 2:331326/150‑1 (MQ=255) gtgaagtgaaACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGt > 2:999588/1‑150 (MQ=255) tgaagtgaaACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCCTTCTGGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTc < 2:1137548/150‑1 (MQ=255) gtgaaACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCtgt > 1:170512/1‑150 (MQ=255) gtgaaACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATAGAACAACTGTCCtgt < 1:115639/150‑1 (MQ=255) tgaaACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCtgtg > 2:41568/1‑150 (MQ=255) cTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTcac < 2:4686/150‑1 (MQ=255) gCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTcacgcacg > 1:1034319/1‑150 (MQ=255) gCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATACTCGACACCCAGCCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTcacgcacg > 2:301721/1‑150 (MQ=255) ccTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACgg > 2:1065261/1‑150 (MQ=255) cTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACTCACGGc < 2:956086/150‑1 (MQ=255) cgtgaaaGGGGCCGCTGC‑GCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGaac < 2:1273259/114‑1 (MQ=255) cgtgaaaGGGGCCGCTGC‑GCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGaac > 1:1273259/7‑120 (MQ=255) gtgaaaGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCt > 1:694721/6‑150 (MQ=255) tgaaaGGGGCCGCTGC‑GCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCtt > 1:335125/5‑150 (MQ=255) aTCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTc > 1:234416/1‑150 (MQ=255) aTCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAACGAAAACCGCATCCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTc > 2:1277187/1‑150 (MQ=255) aaGGGGCCGCTGC‑GCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTccc > 1:381184/2‑150 (MQ=255) ggggCCGC‑TGCGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAgc < 2:607254/44‑1 (MQ=25) ggggCCGC‑TGCGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAgc > 1:607254/1‑44 (MQ=25) | GAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCC > NC_002947/5777318‑5777596 |
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: ‑ |