Is it possible to print DNA sequences of subject organism in tblastn?
0
0
Entering edit mode
5.7 years ago
Kumar ▴ 120

I did tblastn analysis by following command

tblastn -query t3ss.fasta -subject org1.fasta -outfmt '6  qseqid qstart qend sseqid sstart send pident length evalue sseq' -out org1.xls

It prints the subject coordinates and translated amino acid sequences etc., Is there any way to print the DNA sequences of subject organism in tblastn?

BLAST alignment tblastn sequence analysis script • 1.2k views
ADD COMMENT
1
Entering edit mode

You could use the accessions you get back with blastdbcmd to retrieve the DNA sequences, but as mentioned, there’s no way to do this with blast itself.

ADD REPLY
0
Entering edit mode

No. Not if you want to use functionality defined in blast+.

ADD REPLY
0
Entering edit mode

Sorry genomax, I did not get. What is called functionality in blast+

ADD REPLY
1
Entering edit mode

what genomax means to say is that you can not do this with blast+ itself (as in there is no option to request this directly form the blast output).

You could however parse the coordinates you get from the blast output and extract those regions from the original sequences if you want to obtain the DNA sequence from the matched region

ADD REPLY

Login before adding your answer.

Traffic: 2308 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6