Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

feat: write a custom collapsing script #128

Open
deliaBlue opened this issue Dec 2, 2023 · 0 comments
Open

feat: write a custom collapsing script #128

deliaBlue opened this issue Dec 2, 2023 · 0 comments

Comments

@deliaBlue
Copy link
Collaborator

Right now, the workflow uses FASTX Toolkit to collapse the reads. In this processes, the names are changed by numbers but there's no track of which reads belong to which number.

Create a custom script that not only collapses reads, but also a lookup table that lists the new name for each read. In addition, this table can then be used during the uncollapsing step sam_uncollapse.pl to change the read names back to the original ones. Finally, in the step that converts the FASTQ into a FASTA one, names are changed by numbers due to the -r option; this option will have to be removed in order to avoid the issue described here.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant