Comments (2)
I've tried doing this and there's an issue with the code and incompatibility with the way CIGAR strings are encoded for nanopore and pacbio bam data as compared to Bismark processed whole genome bisulfite sequencing bam data.
From what I've looked at it seems long read data uses the bam tags MM and ML.
While bisulfite sequencing and bismark tool encode methylation in bam tags XM.
So some changes in the code seem to need to be made to allow compatibility with MM and ML CIGAR strings instead of only XM CIGAR strings in the bam files. This doesn't seem like a hard thing to do, just someone needs to do it. Please let me know when this is done as I'd like to use it too!
From bismark: https://www.bioinformatics.babraham.ac.uk/projects/bismark/Bismark_User_Guide.pdf
(12) NM-tag (edit distance to the reference) (13) MD-tag (base-by-base mismatches to the reference) (14) XM-tag (methylation call string) (15) XR-tag (read conversion state for the alignment) (16) XG-tag (genome conversion state for the alignment)https://medium.com/@shlokanegi30/mm-and-ml-tags-in-mod-basecalled-bams-using-remora-9fc68a3cb72
MM and ML tags are specific to BAM files, generated by any DNA mod-basecalling algorithm (Eg- Guppy, Megalodon, Bonito, etc.) using nanopore sequencing long-reads.
This blog specifically focuses on understanding MM and ML tags in BAM files generated by Remora for detecting both 5mC and 5hmC mods (DUAL mode)The SAM tags encoding 5mC positions and scores (MM, ML) are added to all HiFi reads.
from wgbs_tools.
wgbstools now supports nanopore files
from wgbs_tools.
Related Issues (20)
- While running bam2pat on ML and MM tag containing bam file ,reads information is printed on the terminal. HOT 1
- bam2pat FileNotFoundError HOT 1
- issue running wgbstools mix_pat
- issue when use `beta_to_450k` HOT 2
- wgbstools convert is not working as expected HOT 2
- Workflow question HOT 1
- Question regarding generating bigger blocks in brain .beta files wgbstools segment
- segment sites mismatch HOT 1
- Failed compiling cview HOT 6
- init_genome error? HOT 1
- beta2bed warning - incomatible (sic) with current genome reference
- Intersect (wgbstools view -L Atlas.U25.l4.hg19.bed Left_Ventricle_STL001.IGF2.pat.gz) HOT 1
- Nanopore wgbstools vis HOT 1
- wgbstools for single base DNA methylation HOT 1
- All ttest values generated from `wgbstools find_markers` are NAs HOT 1
- beta_to_table: np.float deprecated in NumPy>=1.20 HOT 2
- modbam reads and hardclipping
- bam2pat-combine OT and OB?
- Galaxy wrapper
- Bioconda integration
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from wgbs_tools.