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

All reads for filtered out in createASVCF for single end bam files #61

Open
xfanhuang opened this issue May 17, 2024 · 0 comments
Open

Comments

@xfanhuang
Copy link

https://github.com/natsuhiko/rasqual/blob/5aa553cf1b6501cf7ecd61a5efb34f7f20d354c6/src/ASVCF/createASVCF.sh#L121C15-L121C35

I encountered an issue with the AS reads being all zero for a single end bam file, and realized all reads were filtered out by this filter. My understanding is that only properly paired reads will have an "=" sign in column 7, so this filter throws out all single end reads which are marked by "*" in column 7. I feel like this line of filter should be modified for single end mode. Additionally, if the number of reads in the bam file is zero or if AS reads are alarmingly low, a warning could be shown.

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

1 participant