-
Notifications
You must be signed in to change notification settings - Fork 2
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
Error with BAM file on Command 3 #13
Comments
Hi, it is reasonable to expect a few reads from FAST5 to not be present in BAM if basecalling fails for those reads. However, this should be tiny fraction of reads, and almost all reads should be successfully processed. Are you getting the "read not found" error for all the reads or only a fraction of them, and is there a BAM output produced by DeepMod2? And is the number of reads with this error the same as the number of reads in FAST5 files? Another possibility could be a mix up in samples provided to DeepMod2. For instance, the BAM file is from a different sample than FAST5 or the read names have been edited in the BAM in some way? You can check it as follows: Run this command to get read names from the BAM file: Run this command to get read names from the FAST5 files using sequencing_summary.txt file in the basecalled folder: Then check if the read_names_from_bam are the same as read_names_from_basecalls. |
I have indexed and aligned my BAM file still it gives the above error. When I let the command run it produced an empty file.
The text was updated successfully, but these errors were encountered: