-
Notifications
You must be signed in to change notification settings - Fork 7
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
MN tag length and seq length don't match for 99.9% of reads #128
Comments
Hello @ddubocan, The problem comes from Dorado 0.5.1 would not set this value correctly. Normally I would say that you could manually fix the MN tag with something like pysam, but the barcoding itself was off (see the changelog). My collogues have told me that the best thing to do is re-basecall with Dorado 0.5.2 where the issue is fixed. Please let me know if this doesn't fix the error. A |
Thank you so much for the quick response! |
@ddubocan any luck? |
Running modkit(v0.2.5) directly on the output of dorado v0.5.1:
modkit adjust-mods --log-filepath log_modkit_adjust -t 10 --edge-filter 10,10 barcode12.bam barcode12.adjust.bam
only 30 out of 100,000 reads pass, the remaining reads all fail.
When I check the log file:
I do not have this occur when I use modkit (v0.13) -- nearly all reads get processed.
The text was updated successfully, but these errors were encountered: