Hello everyone
I am working with VCF files. I want to extract peptides from mutations that have occurred. for example I have one mutation p.I2447_I2448del and three NM accessions for that in refseq column (NM_015046.7;NM_001351528.2;NM_001351527.1) . I get protein sequence by those accession numbers. now I have three protein sequence and I apply this mutation p.I2447_I2448del for all of them, but in 2447_2448 positions there are SP amino acids instead of II in second protein sequence from NM_001351528.2. So what should I do?
Thanks a lot
1 answer
The three transcripts seem to be from the SETX gene, of which the canonical transcript ID is the first one: NM_015046.7, which translates to the protein NP_055861.3, which has II at positions 2447-2448. You are looking at a different protein product, which is why the co-ordinates don't match. Always look at the protein product from the MANE/canonical transcript first.
A small educational note: if an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one answer if they work. This will help future users that might find this post find the right answer.
Log in to answer this question.
-
How to retrieve STRING preferredName from amino acid sequence using the STRING API?
written by anadomingues •
-
Uniprot Database for download
written by SHN •
-
Getting Intron Positions from Amino Acid Sequences
written by fafad046 •
-
Any script that can produce different mutations of protein sequences at specified positions?
written by maa4618 •
-
Some questions about the CPTAC data processing
written by tujuchuanli
-
Uniprot Accession and Position to peptide sequence
written by desc •
-
mutations' effect on secondary structure
written by lait
-
Associating mutation effect according to mutation
written by Bioaln
-
How does RefSeq Calculate their Positions?
written by pwg46
-
How Is Epitope Mapping Done
written by WinnyH •
You have RNA accession ids for a protein change. What is the NP accession? Could you check your naming here: https://mutalyzer.nl/ ? Also, HGVS uses three letter amino acid instead of single ones.
Thanks for your response. Actually i am using MAF file (other type of vcf file). there is not NP accession number in it. I meant HGVSp_Short.
what did you use to map the VCF to the hgvs?