(bio)-informatics, data processing and visualization

Monday, May 17, 2021

BLAST-N Plus run parameters and tab-delimited output



==================================================

makeblastdb -in Genome_DNA.Fasta -out Genome_DNA.bdb -dbtype nucl -input_type fasta -max_file_sz 2GB -hash_index -parse_seqids 

makeblastdb \
    -in Genome_DNA.Fasta \
    -out Genome_DNA.bdb \
    -dbtype nucl \
    -input_type fasta \
    -max_file_sz 2GB \
    -hash_index \
    -parse_seqids

==================================================

blastn -task blastn -query Query_Seqs.Fasta -db Genome_DNA.bdb -out Query_Seqs.vs.Genome_DNA.blastplus.out.m0 -outfmt 0 -evalue 1e-120 -dbsize 1000000 -dust no -word_size 24 -xdrop_ungap 50 -xdrop_gap 500 -xdrop_gap_final 1000 -max_target_seqs 240 -line_length 100 -num_threads 6

blastn \
    -task blastn \
    -query Query_Seqs.Fasta \
    -db Genome_DNA.bdb \
    -out Query_Seqs.vs.Genome_DNA.blastplus.out.m0 \
    -outfmt 0 \
    -evalue 1e-120 \
    -dbsize 1000000 \
    -dust no \
    -word_size 24 \
    -xdrop_ungap 50 \
    -xdrop_gap 500 \
    -xdrop_gap_final 1000 \
    -max_target_seqs 240 \
    -line_length 100 \
    -num_threads 6

==================================================

blastn -task blastn -query Query_Seqs.Fasta -db Genome_DNA.bdb -out Query_Seqs.vs.Genome_DNA.blastplus.out.m7 -outfmt '7 qseqid sseqid evalue pident score length nident mismatch gaps frames qstart qend sstart send qcovhsp qlen slen' -evalue 1e-120 -dbsize 1000000 -dust no -word_size 24 -xdrop_ungap 50 -xdrop_gap 500 -xdrop_gap_final 1000 -max_target_seqs 240 -num_threads 6 

blastn \
    -task blastn \
    -query Query_Seqs.Fasta \
    -db Genome_DNA.bdb \
    -out Query_Seqs.vs.Genome_DNA.blastplus.out.m7 \
    -outfmt '7 qseqid sseqid evalue pident score length nident mismatch gaps frames qstart qend sstart send qcovhsp qlen slen' \
    -evalue 1e-120 \
    -dbsize 1000000 \
    -dust no \
    -word_size 24 \
    -xdrop_ungap 50 \
    -xdrop_gap 500 \
    -xdrop_gap_final 1000 \
    -max_target_seqs 240 \
    -num_threads 6   

# Fields: 

query id - 1 [0] (qseqid)
subject id - 2 [1] (sseqid)
evalue - 3 [2] (evalue)
% identity - 4 [3] (pident)
score - 5 [4] (score)
alignment length - 6 [5] (length)
identical - 7 [6] (nident)
mismatches - 8 [7] (mismatch)
gaps - 9 [8] (gaps)
query/sbjct frames - 10 [9] (frames)
q. start - 11 [10] (qstart)
q. end - 12 [11] (qend)
s. start - 13 [12] (sstart)
s. end - 14 [13] (send)
% query coverage per hsp - 15 [14] (qcovhsp)
query length - 16 [15] (qlen)
subject length - 17 [16] (slen)

==================================================


No comments:

Post a Comment