diffbind analyze for wheat
1
0
Entering edit mode
zgyu12 • 0
@584cc641
Last seen 7 months ago
China
[E::hts_idx_push] Region 536872557..536872641 cannot be stored in a bai index. Try using a csi index with min_shift = 14, n_lvls >= 6
Error: Error processing one or more read files. Check warnings().
In addition: Warning messages:
1: In mclapply(arglist, fn, ..., mc.preschedule = TRUE, mc.allow.recursive = TRUE) :
  all scheduled cores encountered errors in user code
2:   file: bam/JW_1.uniq.sort.rmdup.1A.bam 



# include your problematic code here with any corresponding output 
# please also include the results of running the following in an R session 

dbObj <- dba(sampleSheet=samples)
dbObj <- dba.count(dbObj, bUseSummarizeOverlaps=TRUE)

I use the hg38 genome for diffbind is OK, while wheat genome is so big, how can use diffbind to analyze the wheat project?

DiffBind • 617 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 19 hours ago
United States

That error doesn't come from DiffBind, but instead it comes from samtools/htslib.

You may be able to index your files using a csi index, but I am not sure that it will then work with DiffBind, as IIRC Rsamtools requires a bai index.

0
Entering edit mode

Thanks, I use cis index for bam, but it raise the same error. May be Rsamtools not support csi index.

ADD REPLY
0
Entering edit mode

Hi, did you ever solve this issue? I have the same issue using axolotl attic-seq data. Curious to know if there's a solution for diff peak analyses using csi indexing.

ADD REPLY

Login before adding your answer.

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