Hi all
I am getting some odd values in my blast outputs. The evalues in the outputs don't seem to correspond properly with the percent_ids (and number of mismatches). I am using Blast 2.2.25 and I have been submitting the jobs using the following command:
#Resources
#$ -pe orte 4
#Run this command:
blastx -query xaa -db /home/my_dir/refseq_dbs/refseq_proteins_eukaryotes.fasta -evalue
0.01 -outfmt 10 -max_target_seqs 1 -out out_a.txt -num_threads 4
I have also tried using the -num_descriptions and the -num_alignments flags instead of the -max_target_seqs flag but this has not fixed the problem.
Here is a few lines from the blastoutput:
qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore
contig4822 281371382 1 80 79 0 326 87 729 808 5.00E-016 86
contig5762 326676144 3 96 91 1 348 67 691 786 4.00E-033 142
contig209 190194343 1 71 70 0 357 145 5 75 1.00E-021 92
contig320 432885829 2 117 115 0 352 2 275 391 4.00E-047 189
contig3304 348520766 1 156 145 3 466 20 1084 1237 3.00E-024 113
Am I missing something?
Thanks in advance!
I am getting some odd values in my blast outputs. The evalues in the outputs don't seem to correspond properly with the percent_ids (and number of mismatches). I am using Blast 2.2.25 and I have been submitting the jobs using the following command:
#Resources
#$ -pe orte 4
#Run this command:
blastx -query xaa -db /home/my_dir/refseq_dbs/refseq_proteins_eukaryotes.fasta -evalue
0.01 -outfmt 10 -max_target_seqs 1 -out out_a.txt -num_threads 4
I have also tried using the -num_descriptions and the -num_alignments flags instead of the -max_target_seqs flag but this has not fixed the problem.
Here is a few lines from the blastoutput:
qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore
contig4822 281371382 1 80 79 0 326 87 729 808 5.00E-016 86
contig5762 326676144 3 96 91 1 348 67 691 786 4.00E-033 142
contig209 190194343 1 71 70 0 357 145 5 75 1.00E-021 92
contig320 432885829 2 117 115 0 352 2 275 391 4.00E-047 189
contig3304 348520766 1 156 145 3 466 20 1084 1237 3.00E-024 113
Am I missing something?
Thanks in advance!
Comment