Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,668,027 | G→T | 8.3% | intergenic (+12/‑29) | PP_mr59 → / → ahcY | SAH_riboswitch regulatory region/adenosylhomocysteinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,668,027 | 0 | G | T | 8.3% | 164.5 / 2.9 | 56 | intergenic (+12/‑29) | PP_mr59/ahcY | SAH_riboswitch regulatory region/adenosylhomocysteinase |
Reads supporting (aligned to +/- strand): ref base G (33/18); new base T (2/3); total (35/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.52e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCTGATGGGCCTGCGTCGCAAGTAC > NC_002947/5667895‑5668175 | gTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGa < 1:1073517/150‑1 (MQ=255) ccACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCt > 2:809133/1‑150 (MQ=255) cacaGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGGTGCAGCAGGCTAGGCCTGTCAGGCTAGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTg < 1:1158426/150‑1 (MQ=255) gCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCaaa < 1:1184604/150‑1 (MQ=255) aTGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCctg > 2:909081/1‑150 (MQ=255) aTGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCctg > 1:1183125/1‑150 (MQ=255) tGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCctgc > 1:1035910/1‑150 (MQ=255) gCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCctgct < 2:1183125/150‑1 (MQ=255) cccGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCctgctg > 1:843772/1‑150 (MQ=255) ccGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGc < 2:803962/120‑1 (MQ=255) ccGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGc > 1:803962/1‑120 (MQ=255) cGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGt > 1:184363/1‑150 (MQ=255) ccACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGttt < 1:538412/150‑1 (MQ=255) ccACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGttt < 2:325913/150‑1 (MQ=255) ccTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACTAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTAc > 2:1086900/1‑150 (MQ=255) tCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGAt > 1:869197/1‑150 (MQ=255) ccGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGAtt > 2:288081/1‑150 (MQ=255) cGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAGCGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTa > 1:819241/1‑150 (MQ=255) cGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTa > 1:614374/1‑150 (MQ=255) gAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTAc > 1:499017/1‑150 (MQ=255) gcgcTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGa < 1:1160351/99‑1 (MQ=255) gcgcTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGa > 2:1160351/1‑99 (MQ=255) cTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGc > 1:282497/1‑150 (MQ=255) cagcagGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGa < 2:184363/150‑1 (MQ=255) aGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCt < 1:964630/150‑1 (MQ=255) gCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACTAATGGAGGCTCTCATGa < 2:274250/93‑1 (MQ=255) gCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACTAATGGAGGCTCTCATGa > 1:274250/1‑93 (MQ=255) cAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCt < 2:869197/150‑1 (MQ=255) aGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACTAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTg < 1:1086900/150‑1 (MQ=255) gCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACCACGGCGCCCATTCGCACATGACTATCCGGGTCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTggg < 1:326992/150‑1 (MQ=255) gCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTggg < 1:909081/150‑1 (MQ=255) cGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGcc < 1:799714/150‑1 (MQ=255) ggCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTgggg > 2:249723/1‑141 (MQ=255) ggCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTgggg < 1:249723/141‑1 (MQ=255) gCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAGCGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGGTTACAAAGTCGCCGACATCTCCCTGACCGCCTGAGGCCGTCGCGaa > 2:374936/1‑150 (MQ=255) gCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCTATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGaa > 2:462997/1‑150 (MQ=255) gTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAAc > 1:89821/1‑150 (MQ=255) tcgctcCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACcatca > 1:194617/1‑150 (MQ=255) ctcCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTgg > 2:672167/1‑120 (MQ=255) ctcCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTgg < 1:672167/120‑1 (MQ=255) tcCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGc > 2:1089053/1‑150 (MQ=255) tcCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGc > 1:1175366/1‑150 (MQ=255) ccGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCg > 1:733323/1‑150 (MQ=255) ccGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCg > 2:357968/1‑150 (MQ=255) cgcgAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGa > 1:1170154/1‑150 (MQ=255) acaACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTgcg > 2:586569/1‑150 (MQ=255) ggCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTACCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCTGAt < 2:1035910/150‑1 (MQ=255) gcCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCTGATggg > 1:632758/1‑150 (MQ=255) gcCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCTGATggg > 1:632672/1‑150 (MQ=255) ccATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCTGATGGGcc > 2:201096/1‑150 (MQ=255) cATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCTGATGGGCCt < 2:1175366/150‑1 (MQ=255) gCACATTACGAATGGAGGCTCTCATCAGCCCAGCAAGCAGGCCTGCTGGTTTTACCGCTTCCAAAGTCGCCGAAATGAGATCGACGAGCACACGTCTGCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCGGAGGGGCCGGCGTc > 1:1178333/1‑150 (MQ=255) aTTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCTGATGGGCCTGCGTCGCaa > 2:888478/1‑150 (MQ=255) ttACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCTGATGGGCCTGCGTCGCAAg > 2:1118465/1‑150 (MQ=255) ttACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCAGACCTCAAGCGCGCCGCCTTGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCTCATGGGCCTGCGTCGCAAg > 2:644574/1‑150 (MQ=255) cGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCTGATGGGCCTGCGTCGCAAGTAc < 2:1178333/150‑1 (MQ=255) | GTTGCCCACAGCGGTGCCATGCCCGCCACCTGTCCGAGGGGCGCTGCAGCAGGCTAGGCCTGTCAGGCTCGGATGGGGCGTTGCTCGCTCCCGCGAGCGCTTAACGCACAACGGCGCCCATTCGCACATTACGAATGGAGGCTCTCATGAGCGCTGCAAACATGCCTGCTGGTTTTACCGATTACAAAGTCGCCGACATCTCCCTGGCCGCCTGGGGCCGTCGCGAAACCATCATTGCCGAATCGGAAATGCCTGCGCTGATGGGCCTGCGTCGCAAGTAC > NC_002947/5667895‑5668175 |
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: ‑ |