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?
python script available here for converting 3 to 1 and 1 and 3. User input is a single column file.
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
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy