diamond
diamond copied to clipboard
--unal 1 incompatible with cigar column
Running this command line (from Galaxy):
diamond blastp --quiet --threads "${GALAXY_SLOTS:-12}" --db ./database --query '/tmp/tmpeyzj5kt3/files/b/b/1/dataset_bb1793d6-dc15-4f75-bd67-09516886d297.dat' --no-self-hits --outfmt '6' qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore cigar scovhsp sskingdoms skingdoms sphylums --out '/tmp/tmpeyzj5kt3/files/b/8/3/dataset_b839f81d-4b20-48c9-9e51-76eaa3179022.dat' --compress '0' --matrix 'BLOSUM62' --comp-based-stats '1' --masking '1' --max-target-seqs '25' --evalue '0.001' --id '0' --query-cover '0' --subject-cover '0' --block-size '2.0' --unal 1
gives
Error: Invalid output field: cigar
in the latest diamond version (2.0.13) .. in 2.0.8 there is also a segfault
Can you also clarify if --unal 1
is compatible only with --outfmt 6
?
The cigar field is bugged for --unal 1
. I have just commited a patch that should fix this.
Thanks