diffbind analyze for wheat
Entering edit mode
zgyu12 • 0
Last seen 16 days ago
[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 • 181 views
Entering edit mode
Last seen 1 day 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.

Entering edit mode

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


Login before adding your answer.

Traffic: 239 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6