This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Map Peptides To Genomic Locations (Hg19)

Is there a reliable tool to map the following sample peptides to genomic locations:

NMGGPYGGGNYGPGGSGGSGGYGGR
SQSSGSSATHPISVPGAR
HTGPNSPDTANDGFVR
NMGGPYGGGNYGPGGSGGSGGYGGR
SKSQSSGSSATHPISVPGAR
SHHAPMSPGSSGGGGQPLAR
GSGTASDDEFENLR

..... .............(many other peptides)

I am interested in the genomic coordinates such as chromosome and the start stop location of the peptide. Thanks

map genomic coordinates

2 answers

There are two strategies: 1. Search your peptides againts human proteins with blastp or blat, then convert the peptide coordinates to the genomic coordinates using exon annotation. 2. Search your peptides translated to nucleotides with tblastn against the reference genome.

Simply do a blastx of your genome with the peptide sequences.

  1. Convert the peptide sequences into multifasta format (you can use excel and microsoft word for this).
  2. Install command line blast. 3.Convert the peptide sequences into a blast database using the makeblastdb tool.
  3. Do a blastx of the genome against the created database.

You will get the coordinates of your peptide sequences mapped on to the genome.

Log in to answer this question.