Calculating short read coverage at genome wide windows for create.set- MEDIPS
0
0
Entering edit mode
@oussamabadad-17646
Last seen 6.0 years ago

Hi all

i am trying to create set for my medipseq data , here is my command

                       

AKM-ANF_MeDIP = MEDIPS.createSet(file = "/scratch/obadad/MEDIP-Seq_BAMs/AKM-ANF-mapped.sorted.bam", BSgenome = "BSgenome.Oeuropaea.IOGC.v1", extend = 300, shift = 0, uniq = 1e-3, window_size = 100)

everything is good till it gets to   Calculating short read coverage at genome wide windows...    


All chromosomes in the reference BSgenome will be processed:
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] "16" "17" "18" "19" "20" "21" "22" "23"
Reading bam alignment AKM-ANF-mapped.sorted.bam
 considering  1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23  using bam index
Total number of imported short reads: 94574662
Extending reads...
Creating GRange Object...
Keep at most 2 read(s) mapping to the same genomic location
Number of remaining reads: 70758499
Calculating genomic coordinates...
Creating Granges object for genome wide windows...
Calculating short read coverage at genome wide windows...
Error in AKM - ANF_MeDIP = MEDIPS.createSet(file = "/scratch/obadad/MEDIP-Seq_BAMs/AKM-ANF-mapped.sorted.bam",  :
  object 'AKM' not found

i have been reading for a while but no clue, do i have to create a MEDIPset object to define AKM-ANF before create.set

do i have to have the reference genome in the same directory as the BAM files ?

i am a bit lost here , i ll appreciate your help.

cheers

oussama

medip-seq createset • 879 views
ADD COMMENT

Login before adding your answer.

Traffic: 729 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