-
Notifications
You must be signed in to change notification settings - Fork 6
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
IGV visualization of bedMethyl file #104
Comments
Hello @MaestSi, Sorry for the delay and thanks for the question. Probably the easiest thing is to convert the bedMethyl into a bedGraph or bigWig file. You would have to decide if you want to combine the top and bottom strands into one track or not, as well as whether or not to combine the modification calls. |
Dear all,
I would like to know if there is a way to properly visualize methylation information obtained from modkit in IGV genome browser (testing IGV 2.17 snapshot).
In particular, I ran Dorado with
--modified-bases 5mC_5hmC 6mA
parameter, I converted the unmapped bam file to fastq withsamtools fastq -T '*’ $UNMAPPED_BAM > $FASTQ
and I mapped the reads to the reference with minimap2 retaining all information in the reads' headers (-y
parameter). Finally, I ranmodkit pileup $BAM modkit.bedMethyl
. I have been able to upload to IGV the bam file, and to color bases by modification. However, I would like to make the bedMethyl file more informative in a visual way. For example, I would like to color the track by %Showing methylation, is it possible? Any advice for better visualization would be greatly apprecciated.Thanks,
SM
The text was updated successfully, but these errors were encountered: