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

Break up long reads into 150mers to enable long read data for variant calling #40

Closed
zachary-foster opened this issue Nov 29, 2023 · 2 comments
Assignees

Comments

@zachary-foster
Copy link
Contributor

A collaborator has tested a lot of methods for variant calling with long read data and determined that breaking up reads into 150mers and using graphtyper works the best of the methods tested. This would be a easy way to accept long reads into the variant calling portion of the pipeline.

@masudermann
Copy link
Contributor

Here is an example command for doing this:
cat $infile | seqkit sliding -W 150 -s 150 -o ${strain}_nano150.fastq.gz

@zachary-foster zachary-foster self-assigned this Apr 22, 2024
@zachary-foster
Copy link
Contributor Author

This has been implemented in the dev branch

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

No branches or pull requests

2 participants