This is a test version of Biostars. For the public version, visit https://www.biostars.org.
using max_target 1,max_hsps 1 for blastn and subsequent compilation of utr files

1) I performed a blastn for a list of mouse genes against the gene database and obtained an output consisting of several thousand hit, I needed only the best hit , which was first on the tabular output

so I tried restricting the result using max_target_seq 1 and max_hsps 1

I got a single hit result. but when I compared it to the unrestricted files I find that the results or the hit shown have different table values such as align.lemgth and %identity although the transcript itself is the same

Can anyone suggest a reason for this and any solutions

2) I used the original to get the first hits for the above

So I have a list of transcript that I obtained from performing a blastn execution. and I have to pull out the corresponding sequences from a source file of 3utr/5utr. but there are a few hundred transcripts and performing it manually is time consuming and error prone

Can anyone suggest a way or command line to perform it using terminal in linux

In summary a method to selectively pull out certain sequence from the gen database file using the ensembl transcript as identifying aspect

max_hsps max_target_seq blast utr
• 4,581 views • link
updated by Ram Ram • written by vigneshprbh37 • vigneshprbh37

1 answer

Hi,

Well, the reason you get different alignment and % identity values is because you have restricted the number of HSPs per hit to 1. Whether you really want to do this or not it depends on what you are doing afterwords with your results. As far as the second question goes, a quick solution would be something like this :

cat blastout | perl -lne '/^(.*?)\t/; $hash{$1} = 1}{ open(IN,"<","seqfile.fa"); while(<IN>){chomp;if(/>(.*)/){$l=($hash{$1} == 1)?(1):(0)}; print "$_" if $l == 1}' > output.fa

but it would be better if you provided a small example. blastout is your tabular output of the blastn search and seqfile.fa is your fasta formatted sequence file with >ID being the header. But as I said it would be better if you could provide some sample files.

m

• 1 views • link
updated by Ram Ram • written by mxs mxs

Log in to answer this question.

More posts like this