Off topic:How createSet from reduced "paire-ed" sequencing on MeDIPS?
0
0
Entering edit mode
pertille • 0
@pertille-9104
Last seen 5.2 years ago
Sweden

I'm trying to create a set with my .BAM file that was generated from reduced "paired end" sequencing (using restriction enzyme), I can easily run it as "single end", but when I add "paired = T" in the command, this message appears:

> MeDIP_RJF_AV163=MEDIPS.createSet(file=RJF_AV163, BSgenome=genome, extend=extend, shift=shift, uniq=uniq, window_size=ws, chr.select=chr_all, paired=T)
Reading bam alignment RJF_AV163_cleaned_sorted.bam
considering  chr1 chr2 chr3 chr4 chr5 chr6 chr7 chr8 chr9 chr10 chr11 chr12 chr13 chr14 chr15 chr16 chr17 chr18 chr19 chr20 chr21 chr22 chr23 chr24 chr25 chr26 chr27 chr28 chr32 chrW chrZ  using bam index
Total number of imported first mate reads in properly mapped pairs: 3855090
scanBamFlag: isPaired = T, isProperPair=TRUE , hasUnmappedMate=FALSE, isUnmappedQuery = F, isFirstMateRead = T, isSecondMateRead = F
Mean insertion size: 148.6603 nt
SD of the insertion size: 71.82463 nt
Max insertion size: 582 nt
Min insertion size: 21 nt
Creating GRange Object...
Error in .Call2("solve_user_SEW0", start, end, width, PACKAGE = "IRanges") :
  solving row 2010: negative widths are not allowed

Some one knows how to figure it out this problem?

 

 

medips R input files • 550 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 655 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