Search
Question: Error running the addCoverage command in QSEA with a custom BSGenome reference
0
gravatar for apt.university
4 months ago by
apt.university0 wrote:

Hi,

I am trying to run the following command: 

qseaSet=addCoverage(qseaSet, uniquePos=TRUE, paired=TRUE)

but I keep on getting the error below.

I am running version sea_1.6.0 and running on a custom BSGenome reference. Any suggestions are greatly appreciated.

Madi

Reading bam alignment combined_bams/S4_Apo1_ALL_sorted.bam

Number of imported sequencing fragments: 13908900

Calculating short read coverage for genome wide windows.

Error in data.frame(seqnames = ReadsL[[idx]]$rname[keep], start = pos[keep] +  :

  arguments imply differing number of rows: 0, 1

In addition: Warning messages:

1: In min(isize) : no non-missing arguments to min; returning Inf

2: In max(isize) : no non-missing arguments to max; returning -Inf
ADD COMMENTlink modified 4 months ago by Matthias Lienhard140 • written 4 months ago by apt.university0
0
gravatar for Matthias Lienhard
4 months ago by
Max Planck Institute for molecular Genetics, Berlin, Germany
Matthias Lienhard140 wrote:

Hi,

this is due to a bug, resulting in the error message if a chromosome/conting is included that does not contain any valid alignment. It is fixed in qsea version 1.7.1.

As a workaround you could also exclude chromosomes without alignments from the chr.select parameter.

Best, Matthias

 

ADD COMMENTlink written 4 months ago by Matthias Lienhard140
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 361 users visited in the last hour