diamond icon indicating copy to clipboard operation
diamond copied to clipboard

--unal 1 incompatible with cigar column

Open bernt-matthias opened this issue 2 years ago • 2 comments

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?

bernt-matthias avatar Nov 25 '21 15:11 bernt-matthias

The cigar field is bugged for --unal 1. I have just commited a patch that should fix this.

bbuchfink avatar Nov 27 '21 12:11 bbuchfink

Thanks

bernt-matthias avatar Nov 29 '21 09:11 bernt-matthias