This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Blast - Formatting Output

Hi,

I've been using the blastn (version 2.2.28+) standalone tool against a custom formatted genome via:

blastn -db BLASTDB -word_size 7 -query input.fa -out filename -perc_identity 100 -outfmt 6 -max_target_seqs 2

To discard non-perfect hits and show only the 2 top hits.

The output file has a great format however is there a way to add an extra column that contains the actual target-seq (sequence of the matched hit)? Such that the fields are:

query id, subject id, % identity, alignment length, mismatches, gap opens, q. start, q. end, s. start, s. end, evalue, bit score, sequence

Thanks!

  • TJC
blast format output
• 185,487 views • link
updated by katieostrouchov • katieostrouchov • written by timjoncooper timjoncooper

thre is a solution to see the sequence(query) in alignment result?

• 0 views • link
written by midox midox

all the valid fields are listed in the help

• 0 views • link
written by Istvan Albert Istvan Albert

I know but in:

       qseqid means Query Seq-id
          qgi means Query GI
         qacc means Query accesion
      qaccver means Query accesion.version
         qlen means Query sequence length
       sseqid means Subject Seq-id

there is one that show the query(sequence)?

• 0 views • link
updated by Ram Ram • written by midox midox

Obviously none of these - after all none of those descriptions indicates that it would. Keep looking.

• 0 views • link
written by Istvan Albert Istvan Albert

Hi!! Do you know how to see the sequence (query) in your blast result?

• 0 views • link
written by figuerm • figuerm

Hi, Is there any way to look for the sequence variation using the above command?

• 0 views • link
written by apekshasharma7777 • apekshasharma7777

2 answers

Run blastn -help then look for the field called outfmt

*** Formatting options
 -outfmt <String>
   alignment view options:
     0 = pairwise,
     1 = query-anchored showing identities,
     2 = query-anchored no identities,
     3 = flat query-anchored, show identities,
     4 = flat query-anchored, no identities,
     5 = XML Blast output,
     6 = tabular,
     7 = tabular with comment lines,
     8 = Text ASN.1,
     9 = Binary ASN.1,
    10 = Comma-separated values,
    11 = BLAST archive format (ASN.1) 

   Options 6, 7, and 10 can be additionally configured to produce
   a custom format specified by space delimited format specifiers.
   The supported format specifiers are:
           qseqid means Query Seq-id
              qgi means Query GI
             qacc means Query accesion
          qaccver means Query accesion.version
             qlen means Query sequence length
           sseqid means Subject Seq-id
        sallseqid means All subject Seq-id(s), separated by a ';'
              sgi means Subject GI
           sallgi means All subject GIs
             sacc means Subject accession
          saccver means Subject accession.version
          sallacc means All subject accessions
             slen means Subject sequence length
           qstart means Start of alignment in query
             qend means End of alignment in query
           sstart means Start of alignment in subject
             send means End of alignment in subject
             qseq means Aligned part of query sequence
             sseq means Aligned part of subject sequence
           evalue means Expect value
         bitscore means Bit score
            score means Raw score
           length means Alignment length
           pident means Percentage of identical matches
           nident means Number of identical matches
         mismatch means Number of mismatches
         positive means Number of positive-scoring matches
          gapopen means Number of gap openings
             gaps means Total number of gaps
             ppos means Percentage of positive-scoring matches
           frames means Query and subject frames separated by a '/'
           qframe means Query frame
           sframe means Subject frame
             btop means Blast traceback operations (BTOP)
          staxids means Subject Taxonomy ID(s), separated by a ';'
        sscinames means Subject Scientific Name(s), separated by a ';'
        scomnames means Subject Common Name(s), separated by a ';'
       sblastnames means Subject Blast Name(s), separated by a ';'
                (in alphabetical order)
       sskingdoms means Subject Super Kingdom(s), separated by a ';'
                (in alphabetical order) 
           stitle means Subject Title
       salltitles means All Subject Title(s), separated by a '&lt;&gt;'
          sstrand means Subject Strand
            qcovs means Query Coverage Per Subject
          qcovhsp means Query Coverage Per HSP
• 1,013 views • link
updated by Ram Ram • written by Istvan Albert Istvan Albert

outfmt 7 or 10 works perfect

• 0 views • link
updated by Ram Ram • written by H@rry • H@rry

Thank you! Sorted it out now.

• 0 views • link
updated by Ram Ram • written by timjoncooper timjoncooper

Hi!!
I have the same question and I don´t know how you sort it out? Was it that you used oufmt 7 or that you use -outfmt "6 qlen" ??

• 0 views • link
updated by Ram Ram • written by figuerm • figuerm

To clarify by "space delimited format specifiers", it means write it as -outfmt "6 qacc sacc qseq sseq..."

• 0 views • link
updated by Ram Ram • written by ostrokach ostrokach

To add, one of format specifiers is std, which add there default set. It means that -outfmt "6 std qlen" prints standard and query length.

• 0 views • link
updated by Ram Ram • written by kamiljaron kamiljaron

How can I get the description ( first column in the figure) when I run the command line blastp?

blastp

• 1 views • link
written by BDK_compbio BDK_compbio

The description can be added by "stitle".

• 0 views • link
written by mitchez • mitchez

How to give a mismatch parameter in blastn. I was to perform alignment allowing 1 mismatch. I'm going through a lot of parameters but can't find this one.

• 0 views • link
written by glady glady

Hi, Is there a way to find query strand information as well? Thanks

• 0 views • link
written by Neu • Neu

I believe that strand is the relative position of subject to query, hence if sstrand is reverse, it that the query reverse complementary to the reference sequence.

• 0 views • link
written by kamiljaron kamiljaron

Hi . Can you please suggest How to use output format 8 . -outfmt 8 doesn't work. I am trying to do blastp and have tried using -m 8 and -outfmt 8. It doesn't work though it works with -outfmt 6.

• 0 views • link
written by shilpi.singh2811 • shilpi.singh2811

Just so you know, I was looking for this as well. The following did the job for me:

blastn -db <db_source> -query <query_source> -out <outfile> -outfmt "6 qseqid sseqid slen qstart qend length mismatch gapopen gaps sseq"  -word_size 5 -perc_identity 80

The option "sseq" will give the sequence that the query was aligned with. The option "qseq" will be the part of the query sequence.

• 0 views • link
written by sridhar.rg • sridhar.rg

Log in to answer this question.

More posts like this