Hi!
We have a series of mutations that are referenced to HG19. We would like to know how these mutations change the HG19 Amino Acids. We would like to have a resulting file with both the original and the new amino acid created by the mutation.
We are starting with files that have 1000s of mutations and would like to process these files all at once.
Is there a way to use the "aaChanges" function to accomplish this? We tried to do this but the aaChanges will apparently only return known SNP aa changes (as opposed to aa changes caused by novel mutations).
Thank you!
I am using this as my input - but it contains mutations and SNPs. But my output only has the SNP positions in it.
chr2 179395067 179395068 C/G
chr2 179440029 179440030 G/A
chr2 179485919 179485920 G/A
chr2 179506790 179506791 A/T
chr2 179528807 179528808 T/A
chr2 179535087 179535088 T/C
chr2 179536884 179536885 C/T
chr2 179543705 179543706 G/A
chr2 179554305 179554306 C/T
chr2 179566038 179566039 C/T
chr2 179566051 179566052 C/T
chr2 179611205 179611206 G/T
chr2 179635605 179635606 A/G
chr2 179641951 179641952 G/A
chr2 179643566 179643567 A/G
chr2 179650701 179650702 C/T
chr2 179666830 179666831 G/A
chr2 179667090 179667091 C/T
Do you have the coordinates of the SNP?
Hi! Thank you for your quick response!
Yes I am using this as my input - but it contains mutations and SNPs. But my output only has the SNP positions in it.
chr2 179395067 179395068 C/G
chr2 179440029 179440030 G/A
chr2 179485919 179485920 G/A
chr2 179506790 179506791 A/T
chr2 179528807 179528808 T/A
chr2 179535087 179535088 T/C
chr2 179536884 179536885 C/T
chr2 179543705 179543706 G/A
chr2 179554305 179554306 C/T
chr2 179566038 179566039 C/T
chr2 179566051 179566052 C/T
chr2 179611205 179611206 G/T
chr2 179635605 179635606 A/G
chr2 179641951 179641952 G/A
chr2 179643566 179643567 A/G
chr2 179650701 179650702 C/T
chr2 179666830 179666831 G/A
chr2 179667090 179667091 C/T
Can you extend your original question with this information. Thanks.