Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,845,583 | G→A | 7.5% | S17S (AGC→AGT) | PP_3394 ← | 3‑hydroxy‑3‑methylglutaryl‑CoA lyase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,845,583 | 0 | G | A | 7.5% | 151.2 / 2.3 | 53 | S17S (AGC→AGT) | PP_3394 | 3‑hydroxy‑3‑methylglutaryl‑CoA lyase |
Reads supporting (aligned to +/- strand): ref base G (20/29); new base A (2/2); total (22/31) | |||||||||||
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.35e-01 |
GCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGA > NC_002947/3845443‑3845730 | gCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCaa < 7:56433/149‑1 (MQ=255) caccacTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGc > 2:107943/1‑149 (MQ=255) caccacTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGc > 2:103505/1‑149 (MQ=255) ttCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCaa < 6:63023/141‑1 (MQ=255) ttCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCaa > 5:63023/1‑141 (MQ=255) cgcgGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATcgcg > 3:44824/1‑148 (MQ=255) gggCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAgg > 3:20359/1‑149 (MQ=255) gggCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAgg > 2:11509/1‑149 (MQ=255) aTCGGCCATCCGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGACATCGCGCAGGCCGa > 6:266525/1‑149 (MQ=255) tCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATGACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGAc < 7:192269/149‑1 (MQ=255) tCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGAc < 7:209479/149‑1 (MQ=255) gCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACtt < 2:97697/148‑1 (MQ=255) gCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACtt > 1:97697/1‑148 (MQ=255) gCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTc < 5:183591/149‑1 (MQ=255) gCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAgg < 7:28621/113‑1 (MQ=255) gCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAgg > 8:28621/1‑113 (MQ=255) cGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCattgattg > 7:214984/1‑149 (MQ=255) cGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCattgattg < 1:107943/149‑1 (MQ=255) cGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCattgattg < 1:103505/149‑1 (MQ=255) gCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGAGCCTGCAAgt > 7:142749/1‑112 (MQ=255) gCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGAGCCTGCAAgt < 8:142749/113‑2 (MQ=255) cGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATtgttgt < 8:172195/149‑1 (MQ=255) gTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTc < 4:20359/149‑1 (MQ=255) gCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATcgcg > 5:319216/1‑149 (MQ=255) gaacgaCGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGtct < 1:193665/149‑1 (MQ=255) gaacgaCGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGAGCCTGCAAgt < 8:150089/101‑2 (MQ=255) gaacgaCGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGAGCCTGCAAgt > 7:150089/1‑100 (MQ=255) gacgCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGtct < 6:3933/145‑1 (MQ=255) gacgCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGtct > 5:3933/1‑145 (MQ=255) aCTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGa < 3:219115/149‑1 (MQ=255) cTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGaa < 5:164670/149‑1 (MQ=255) ttCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAg > 1:182412/1‑149 (MQ=255) cATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCg < 6:187438/148‑1 (MQ=255) atatGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGtt > 2:264978/1‑149 (MQ=255) atGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCCGaa < 5:354684/142‑1 (MQ=255) atGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCCGaa > 6:354684/1‑142 (MQ=255) gcgcACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCtctg > 4:94751/1‑133 (MQ=255) gcgcACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCtctg < 3:94751/133‑1 (MQ=255) cgcACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCa < 5:266525/149‑1 (MQ=255) ccggcgtgcgcCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGa < 6:152335/142‑1 (MQ=255) cACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGa < 1:11509/149‑1 (MQ=255) cGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTgg < 6:319216/149‑1 (MQ=255) tAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAt > 2:126405/1‑149 (MQ=255) tCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGaaaa < 6:290082/149‑1 (MQ=255) cAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGa < 4:44824/149‑1 (MQ=255) gcggcgTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTaa > 1:348167/1‑123 (MQ=255) gcggcgTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTaa < 2:348167/123‑1 (MQ=255) gTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGttt > 6:123493/1‑149 (MQ=255) tgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑gggc > 6:142356/1‑147 (MQ=255) cATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATc < 1:126405/149‑1 (MQ=255) aTCATTTGCAGGTTCTGCATGCCTTCGCGAAGGCCTACTTCATTGATTTTAATTGTTGTCATCGGGTGCTCTCTGTCGTACGCGTTGCGCATGACTAATAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTc < 1:264978/148‑1 (MQ=255) aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTc > 1:299621/1‑145 (MQ=255) aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTc < 2:299621/145‑1 (MQ=255) cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt < 5:123493/149‑1 (MQ=255) cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt < 2:182412/149‑1 (MQ=255) gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCACCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGa < 1:8091/149‑1 (MQ=255) | GCCACCACTTCGCGGGCATCGGCCATCTGCGGCAGCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGA > NC_002947/3845443‑3845730 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |