Use bismark files to run methylKit
1
0
Entering edit mode
16 months ago
anasjamshed ▴ 120

Can anyone provide me with some bam or sam bismark files to run methylKit on it?

I want to execute processBismarkAln()

bismark R • 1.5k views
ADD COMMENT
1
Entering edit mode
16 months ago
Basti ★ 2.0k

There are test datasets in the methylKit package : https://www.rdocumentation.org/packages/methylKit/versions/0.99.2/topics/processBismarkAln

my.file=system.file("extdata", "test.fastq_bismark.sorted.min.sam", 
                              package = "methylKit")
obj=processBismarkAln(my.file,"test",assembly="hg18",save.folder=NULL,
                 save.context="CpG",read.context="CpG")
ADD COMMENT
0
Entering edit mode

Thanks. How can I use it?

ADD REPLY
0
Entering edit mode

I do not understand your question, I provided the code to test processBismarkAln with a sam test file

ADD REPLY
0
Entering edit mode

but where are these files?

ADD REPLY
1
Entering edit mode

In methylKit package, in the extdata folder as I provided an example of sam file with my code :

my.file=system.file("extdata", "test.fastq_bismark.sorted.min.sam", 
                              package = "methylKit")

Complete list of all test datasets available for this package can be shown with :

list.files(system.file("extdata",package="methylKit"))
 [1] "bismark_coverage_CpG.bed"                 "bismark_cytosineReport.txt"              
 [3] "control1.myCpG.txt"                       "control2.myCpG.txt"                      
 [5] "cpgi.hg18.bed.txt"                        "ctrl.bismark_paired_end.sorted.bam"      
 [7] "ctrl.bismark_paired_end.sorted.wh.sam"    "ctrl1.txt.bgz"                           
 [9] "ctrl2.txt.bgz"                            "generic1.CpG.txt"                        
[11] "refseq.hg18.bed.txt"                      "test.bismark_single_end.sorted.bam"      
[13] "test.bismark_single_end.sorted.wh.sam"    "test.fastq_bismark.sorted.min.sam"       
[15] "test.fastq_bismark.unsorted_chr.min.sam"  "test.fastq_bismark.unsorted.min.sam"     
[17] "test1_test2_ctrl1_ctrl2_diffMeth.txt.bgz" "test1_test2_ctrl1_ctrl2.txt.bgz"         
[19] "test1.myCpG.txt"                          "test1.txt.bgz"                           
[21] "test2.myCpG.txt"                          "test2.txt.bgz"            
ADD REPLY
0
Entering edit mode

ok. I tried and get this error: enter image description here

ADD REPLY
0
Entering edit mode

This is not an error message

ADD REPLY
0
Entering edit mode

ok, but my question is when I try to run my own file, it gives me an error:

Trying to process:
test.sorted.bam
Using htslib.
Error: Unexpected cigar: M

Why I get this error?

ADD REPLY

Login before adding your answer.

Traffic: 2677 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6