Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,582,005 | G→T | 10.7% | E26* (GAG→TAG) | benA → | benzoate 1,2‑dioxygenase subunit alpha |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,582,005 | 0 | G | T | 10.7% | 155.2 / 7.6 | 55 | E26* (GAG→TAG) | benA | benzoate 1,2‑dioxygenase subunit alpha |
Reads supporting (aligned to +/- strand): ref base G (25/24); new base T (2/4); total (27/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.69e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01 |
AATCCGGATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATGGGGCGCCAGCCGATTTT > NC_002947/3581857‑3582144 | aaTCCGGATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGAGTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCggg > 2:354954/1‑149 (MQ=255) cGGATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTCCCGCTTGCAGCGGTAGAt < 7:114965/149‑1 (MQ=255) ggATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTCCCGCTTGCAGCGGTAGATg < 3:340534/149‑1 (MQ=255) gATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGt > 8:10957/1‑149 (MQ=255) aTAGCCCCCCGCCAGCCGGCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTCCCGCTTGCAGCGGTAGATgcc < 4:314326/149‑3 (MQ=255) tAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTc < 2:133706/149‑1 (MQ=255) ccccccGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCAcc > 1:261737/1‑149 (MQ=255) cacccGTCTCCCTAACCTGACCACGTCTAAACAATAACCAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCCAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACg < 6:111687/146‑1 (MQ=255) aGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACgg < 5:8816/149‑1 (MQ=255) tctcCCTAACCAGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATACGACTACCTCAATGCCATGCTCGAGGACGACCGTGATAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGtt > 2:248305/1‑149 (MQ=255) ctcCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTc < 5:350055/149‑1 (MQ=255) ccTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGa > 1:188994/1‑148 (MQ=255) ccTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGa < 2:188994/148‑1 (MQ=255) ccTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGAc > 1:328127/1‑149 (MQ=255) aCCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCACTGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGTAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCt < 5:29747/147‑1 (MQ=255) cTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTg > 3:220882/1‑146 (MQ=255) aCCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGAt < 2:328127/148‑1 (MQ=255) aCCAGGTCTAAACAATAACAAGGGAGACACTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGATCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGAt > 7:265203/1‑148 (MQ=255) ccAGGTCTAAACAATAACAAGGGAGACCCTGGCCATTGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGgag > 1:54426/1‑113 (MQ=255) ccAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGgag < 2:54426/112‑1 (MQ=255) cAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGgaga > 1:156295/1‑145 (MQ=255) cAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGgaga < 2:156295/145‑1 (MQ=255) aGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGaaa < 7:303956/149‑1 (MQ=255) tAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAAcaca > 4:129097/1‑148 (MQ=255) tAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACAt > 7:92101/1‑149 (MQ=255) aaCAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTAGACTACATCAATGCCATGCTCGAGGACGAACGTGAAAAGGGCATCAACCCCTGCAAGCGGTAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCt > 1:27355/1‑149 (MQ=255) aaGGGAGACCCTGGCCA‑TGTCCCTGGGATTTGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGc > 4:13389/1‑148 (MQ=255) agaCCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTgg < 2:261737/149‑1 (MQ=255) tCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCGCCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGcc < 8:251962/147‑1 (MQ=255) tCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCGCCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGcc > 7:251962/1‑147 (MQ=255) ccTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCAc > 5:347786/1‑148 (MQ=255) ggATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGaaa < 4:247688/148‑1 (MQ=255) gATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGc < 6:347786/149‑1 (MQ=255) ttCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCa < 4:220882/149‑1 (MQ=255) tagcTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACCGGATCTACCTGGCCCACGAAAGCCAGATccc > 4:299326/4‑149 (MQ=255) ccTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTGCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATccc < 8:79012/147‑1 (MQ=255) cTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCGACTGGATCTACCTGGCCCACGAAAGCCAGATcccc > 8:236744/1‑147 (MQ=255) cTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATcccc > 5:54427/1‑147 (MQ=255) tCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATcccc > 3:183561/1‑146 (MQ=255) cAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACg < 4:55977/131‑1 (MQ=255) cAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACg > 3:55977/1‑131 (MQ=255) aTGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGagaag < 2:27355/149‑1 (MQ=255) ccATGGTCGAGGACGACCGTTAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAAt < 8:265203/149‑1 (MQ=255) ccATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAAt < 6:54427/149‑1 (MQ=255) tCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTc > 5:322205/1‑149 (MQ=255) cGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTcc > 4:15769/1‑149 (MQ=255) gAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTc < 1:354954/147‑1 (MQ=255) acgacCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCAcc < 4:183561/149‑1 (MQ=255) cGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGAc < 7:10957/149‑1 (MQ=255) ggCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATGGGgcg > 1:277717/1‑149 (MQ=255) cATCTCCCGCTTGCAGCGGTAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCg > 2:85873/1‑113 (MQ=255) cATCTCCCGCTTGCAGCGGTAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCg < 1:85873/113‑1 (MQ=255) tCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGCTCCCCGAGAAGAATGAGTTCCTCACCCTGACGATGGGGCGcccg > 4:113818/1‑147 (MQ=255) tCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATg > 3:149998/1‑140 (MQ=255) tCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATg < 4:149998/140‑1 (MQ=255) aCCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGaaa > 8:3342/1‑96 (MQ=255) ccGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGCAACACCTCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGCCTTCCTCACCCTGCCCATGGGGCGCCAGCCGa > 7:294092/1‑149 (MQ=255) cAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCTACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATGGGGCGCCAGCCGAtttt < 7:236744/147‑1 (MQ=255) | AATCCGGATAGCCCCCCGCCAGCCGTCTCCCTAACCTGACCAGGTCTAAACAATAACAAGGGAGACCCTGGCCA‑TGTCCCTGGGATTCGACTACCTCAATGCCATGCTCGAGGACGACCGTGAAAAGGGCATCTACCGCTGCAAGCGGGAGATGTTCACCGACCCACGGCTGTTCGACCTGGAGATGAAACACATCTTCGAGGGCAACTGGATCTACCTGGCCCACGAAAGCCAGATCCCCGAGAAGAATGACTTCCTCACCCTGACCATGGGGCGCCAGCCGATTTT > NC_002947/3581857‑3582144 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |