This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Use Blast To Find Exact Matches Of Short Sequences?

Hi, I am using tblastn (under blast 2.2.25+) for exact peptide mapping (no gaps).
I want to map few peptides (about 6 to 50 AAs in length) to genome.
However, as I test a known peptide of 6 AAs,tblastn failed to mapped this peptide.
I have read the doc of blast, but failed to find a solution. What did I miss?
Thank you!
PS. I have also tried PGM (ProteogenomicMapping). This tool can map the known peptide tested above correctly, but it's slow in my computer which is impossible for large scale mapping.

blast
• 16,713 views • link
updated by SRKR SRKR • written by Free Man Free Man

1 answer

In the command line BLAST there is an option -perc_identity. You can use this, keep it as 100 and then run the blast. With that setting you will be able to get hits only if there is 100% identity. The command would be like this:

blastn -db dbname -query input_file -out output_file -perc_identity 100

you can try this and I believe it will work. you can also use word-size to get hits with even shorter peptides, like your case. It's value should be a minimum of 2 in case of tblastn

-word_size 3

you can always get to know all the options available by typing -h (brief) or -help (detailed) after the blast type

tblastn -help

hope this helps...

• 146 views • link
written by SRKR SRKR

Hi, I got error: "Error: (CArgException::eInvalidArg) Unknown argument: "perc_identity"".
I did not find something like 'perc_identity' in the help doc for tblastn. It seems it is only avaliable for blastn. So what version are you using?

• 156 views • link
written by Free Man Free Man

yeah I am sorry, just now noticed that -perc_identity is not available with tbalstn. The best option that seems to be the case is to use -ungapped, which will avoid gaps, but still it might result in mismatches.

• 130 views • link
written by SRKR SRKR

What is your genome size? If it isn't too big a script can be useful to you to get the positions. Just have to six frame translate the genome and search for your amino acid sequences in the translates. You will get the positions all through the genome.

• 0 views • link
written by SRKR SRKR

Thanks for you suggestion! After tedious attempts using various parameters, I got the solution for my project:
Key parameters: -comp_based_stats 0 -ungapped -matrix PAM 30 -seg no

• 0 views • link
written by Free Man Free Man

Log in to answer this question.

More posts like this