diamond
diamond copied to clipboard
Feature Request: Recovery of untranslated reads after translated alignment
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!
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.