Commandline BLAST – errors?

Commandline BLAST – errors?

0

Hi,

I’m running command line blastx and blastp against a number of databases. However, running the exact same script on the exact same input files against the exact same databases occasionally seems to output different filesizes.

I can only assume that this is because the task is failing – but how can I QC the outputs to tell which output file was successful? I know I can look at the log files, but some of these I am inheriting and I want to be able to double check if they were successful.

The outputs tend to be outfmt6 or asn then converted to asn (I’m following a specific methodology, this isn’t my choice)

blastx 
-query /input-file.fa 
-db nr 
-num_threads 24 
-max_target_seqs 10 
-outfmt 11 
-out output-file.asn

blast_formatter 
-archive output-file.asn 
-outfmt "6 qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore" > 
output-file.outfmt6

Thanks!


line


BLAST


RNA-seq


command

• 12 views

updated 19 minutes ago by

★

1.0k

written 2 hours ago by

0

Read more here: Source link