Convert protein name to letter e.g. p.Arg135Trp to p.R135W
1
0
Entering edit mode
6 months ago
whb ▴ 20

How to convert snpeff output of protein change from p.Arg135Trp to p.R135W? so that I can import into maftools,

How can it be done in R or bash?

Thanks

aminoacid • 345 views
ADD COMMENT
0
Entering edit mode

python script available here for converting 3 to 1 and 1 and 3. User input is a single column file.

https://github.com/svsuresh/python_scripts/blob/master/hgvs_p_syntax_conversion.ipynb

ADD REPLY
2
Entering edit mode
6 months ago
tothepoint ▴ 600

Create a text file with all the input (Your_file.txt) and this will help you to get the desired output.

sed 's/Ala/A/g;s/Cys/C/g;s/Asp/D/g;s/Glu/E/g;s/Phe/F/g;s/Gly/G/g;s/His/H/g;s/Hid/H/g;s/Hie/H/g;s/Ile/I/g;s/Lys/K/g;s/Leu/L/g;s/Met/M/g;s/Asn/N/g;s/Pro/P/g;s/Gln/Q/g;s/Arg/R/g;s/Ser/S/g;s/Thr/T/g;s/Val/V/g;s/Trp/W/g;s/Tyr/Y/g;s/Mse/X/g' < Your_file.txt > desired_output.txt
ADD COMMENT

Login before adding your answer.

Traffic: 2240 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6