`view qseq` uses letters `ARND` instead of `ACTG`
Opened this issue · 5 comments
imLew commented
User diamond version 2.1.7 diamond view --daa ... --outfmt qseq
produces sequences that contain the letters ARND
instead of ACTG
.
When I run diamond --outfmt qseq
on the fastq files used to produce the daa file this does not happen. Furthermore the sequences have the same lengths and is a 1 to 1 correspondence between the letters (save a few N
s in the sequences using ACTG
), so I don't think they are ambiguity codes.
imLew commented
example.zip
Here is an example running diamond view --outfmt 6 qseq -a
this should produce sequences with ARND
.
terrycojones commented
@bbuchfink Hi Benjamin. Do you have any idea about this issue? I hope you're well! Terry
bbuchfink commented
It will be fixed in the next release.
terrycojones commented
Thanks!
imLew commented
Thank you!