diamond icon indicating copy to clipboard operation
diamond copied to clipboard

Feature Request: Recovery of untranslated reads after translated alignment

Open BenjaminJPerry opened this issue 3 years ago • 1 comments

Hello,

Is it possible to recovery to an output file the original DNA/RNA reads that have aligned to a protein database?

For example, I have a known protein and want to aligned a paired-end sequencing read dataset to the protein and recover the original read that aligned after translation.

Thanks!

BenjaminJPerry avatar Jun 10 '21 01:06 BenjaminJPerry

You can retrieve the read from your query file using the sequence id that diamond reports. You can also have diamond report the query sequence in the tabular output using the full_qseq output field.

bbuchfink avatar Jun 12 '21 07:06 bbuchfink