Error running the addCoverage command in QSEA with a custom BSGenome reference
1
0
Entering edit mode
@aptuniversity-15979
Last seen 6.4 years ago

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
qsea • 968 views
ADD COMMENT
0
Entering edit mode
@matthias-lienhard-6292
Last seen 10 months ago
Max Planck Institute for molecular Gene…

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 COMMENT

Login before adding your answer.

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