Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,720,208 | C→A | 7.4% | R154R (CGC→CGA) | paaY → | phenylacetate degradation pathway detoxifying thioesterase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,720,208 | 0 | C | A | 7.4% | 169.5 / 3.8 | 54 | R154R (CGC→CGA) | paaY | phenylacetate degradation pathway detoxifying thioesterase |
Reads supporting (aligned to +/- strand): ref base C (24/26); new base A (2/2); total (26/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01 |
GCGCCCATGTAGCACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGTCAGGCCCAAGGGCCAGGCAAGCGCATGAGCCTGGC > NC_002947/3720063‑3720353 | gcgcCCATGTAGCACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCgg < 4:307340/148‑1 (MQ=255) ccATGTAGCACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGAcgg < 1:317450/149‑1 (MQ=255) aTGTAGCACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGcgcg < 1:171337/141‑1 (MQ=255) aTGTAGCACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGcgcg > 2:171337/1‑141 (MQ=255) aTGTAGCACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGAGGGAcggcg < 3:129777/149‑1 (MQ=255) aTGTAGCACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGAGGGAcggcg < 3:134995/149‑1 (MQ=255) aTGTAGCACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGAGGGAcggcg > 4:134995/1‑149 (MQ=255) aTGTAGCACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGAGGGAcggcg > 4:129777/1‑149 (MQ=255) aGCACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAg > 3:94874/1‑148 (MQ=255) cACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTAc > 5:133626/1‑149 (MQ=255) aCCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTAc > 8:19665/1‑148 (MQ=255) cTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCAcc > 4:240261/1‑149 (MQ=255) tcgtcgCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGAcgg > 6:280189/1‑129 (MQ=255) tcgtcgCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGAcgg < 5:280189/129‑1 (MQ=255) gtcgCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCa < 8:78223/149‑1 (MQ=255) ccACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGc < 4:94874/149‑1 (MQ=255) cACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCAc > 1:34427/1‑133 (MQ=255) cACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCAc < 2:34427/133‑1 (MQ=255) cACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCa > 3:273657/1‑149 (MQ=255) cACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCa > 6:75154/1‑149 (MQ=255) cACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCa < 4:273657/149‑1 (MQ=255) tCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAg > 3:141094/1‑148 (MQ=255) tCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAg < 3:225880/131‑1 (MQ=255) tCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAg > 4:225880/1‑131 (MQ=255) aaGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTc > 6:330635/1‑149 (MQ=255) tGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAAt > 7:19565/1‑149 (MQ=255) cTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGcc < 7:164614/149‑1 (MQ=255) cAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGa < 5:227301/86‑1 (MQ=255) cAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGa > 6:227301/1‑86 (MQ=255) gagCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAgg > 6:348572/1‑146 (MQ=255) gCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCagag < 5:348572/149‑1 (MQ=255) ccTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGc > 6:67306/1‑149 (MQ=255) gTGATGGGGGCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCggg < 5:67306/149‑1 (MQ=255) gTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCcgc > 5:246457/1‑149 (MQ=255) gCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCa > 8:283571/1‑149 (MQ=255) ccAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCAt < 6:133626/149‑1 (MQ=255) cAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATg < 4:141094/149‑1 (MQ=255) aGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATg > 5:343274/1‑148 (MQ=255) gCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATgg > 4:86411/1‑148 (MQ=255) gCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGa < 7:19665/149‑1 (MQ=255) ccAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGa < 4:130350/147‑1 (MQ=255) ccAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGa > 3:130350/1‑147 (MQ=255) aGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGa > 7:2262/1‑148 (MQ=255) aGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGa > 4:203558/1‑148 (MQ=255) aaaCGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCg < 8:19565/148‑1 (MQ=255) aaCGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCAAGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCggg < 3:240261/149‑1 (MQ=255) cGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGt < 6:343274/149‑1 (MQ=255) cGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGt < 2:5464/149‑1 (MQ=255) ccACTGAGCTAGCAGGAGCTGGCCTTGAAACAGCGCGGGACGGCGGAGTACGAGCAACTGGAGCAGCTCTGCATGATCAGCTTTGTCGAATGCCGACCGCTGGCCGAGGCTGAGCCGGGGCGGCCGCGCATGGAAGACACCGGTGTctg < 2:116156/149‑3 (MQ=255) agcgagcAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGTCAGGccc < 5:75154/147‑1 (MQ=255) cgagcAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGTCAGGCCCAAgg < 6:151427/149‑1 (MQ=255) cgagcAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGTCAGGCCCAAgg > 8:315220/1‑149 (MQ=255) gAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGTCAGGCCCAAGGGCCAgg > 6:43352/1‑148 (MQ=255) gAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGTCAGGCCCAAGGGCCAGGc < 8:2262/149‑1 (MQ=255) gCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGTCAGGCCCAAGGGCCAGGCaa < 7:283571/149‑1 (MQ=255) cTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGTCAGGCCCAAGGGCCAGGCAAGCGCAt < 3:86411/149‑1 (MQ=255) gcgcGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGTCAGGCCCAAGGGCCAGGCAAGCGCATGAGCCTGGc < 5:43352/149‑1 (MQ=255) | GCGCCCATGTAGCACCTCGCTGCATCGTCGCGGCCACCGCCTTCGTCAAAGCTGGCTTCGAATGTGCGGCGCAGAGCCTGGTGATGGGGTCGCCAGCCCAGGTCAAACGCCCACTGAGCGAGCAGGAGCTGGCCTGGAAACAGCGCGGGACGGCGGAGTACCAGCACCTGGCGCAGCGCTGCATGAACAGCATGGTCGAATGCCCACCGCTGGCCGAGGCAGAGCCGGGGCGCCCGCGCATGGAAGACACCGGGGTCAGGCCCAAGGGCCAGGCAAGCGCATGAGCCTGGC > NC_002947/3720063‑3720353 |
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: ‑ |