Blast command line
0
0
Entering edit mode
2.2 years ago

Hello all,

  • I would like to run a blastx on a fasta file in order to perform a quick functionnal annotation on swissprot.

  • So I ran succesfully :

    blastx -query "${QUERY}" -db "${BANK}" -out ${OUT_FILE} -max_target_seqs 1 ${BLAST_PARAM} -num_threads $NCPUS -outfmt "6 qseqid sseqid sacc pident length evalue bitscore"
    
  • However the output for my first line looks like this :

       TRINITY_DN80826_c0_g1_i1        PIAS1_HUMAN     PIAS1_HUMAN     58.879  107         2.52e-36        132
    

First don't understand why there is no hit accession in column 3 and : it is supposed to be given by the option sacc, and also there is one in blast web version !

Second, when I run web blast on web version I have as a result in description column :

RecName: Full=E3 SUMO-protein ligase PIAS1; AltName: Full=DEAD/H box-binding protein 1; AltName: Full=E3 SUMO-protein transferase PIAS1; AltName: Full=Gu-binding protein; Short=GBP; AltName: Full=Protein inhibitor of activated STAT protein 1; AltName: Full=RNA helicase II-binding protein [Homo sapiens]

How to have this information command line and not only a gene name like I do ?

Thank you for your help !

RNA-Seq software error sequence • 506 views
ADD COMMENT
1
Entering edit mode

Take a look at the options for -outfmt 6 and make changes to options you are adding. Full list of options is here.

ADD REPLY
1
Entering edit mode

things like PIAS1_HUMAN are the swissprot accession IDs , so you actually got what you asked for. For more informative info do what @genomax said and change the output you request from the blast run (eg. use sallacc or the default output )

ADD REPLY
0
Entering edit mode

I have tried sallacc and it gives me also the swissprot accession. Is there any way to also get the ncbi accession (for instance O75925 for PIAS1_HUMAN) ?

ADD REPLY
1
Entering edit mode

is that a DB you created yourself locally?

if so (and seeing that adding sallacc only returns the swissprot ID) it might be you did not format it correctly ...

ADD REPLY
0
Entering edit mode

Does your db actually include such info? Maybe there's a map file you can later join with the blast output?

ADD REPLY

Login before adding your answer.

Traffic: 1369 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