ATACseqQC shiftGAlignmentsList leads to (all(elementNROWS)>3) is not TRUE error
0
0
Entering edit mode
mbasam • 0
@mbasam-21086
Last seen 19 months ago

All other parts of the ATACseqQC package works, but when the shiftGAlignmentsList command is ran to shift the reads and split the BAM files, the command fails with the following error: (all(elementNROWS)>3) is not TRUE. However, in R, when I run that line, it returns as TRUE. What is the exact error the script is running into? Is it a possible bug?

ATACseqQC ATAC ATACseq software error Tutorial • 464 views
ADD COMMENT
0
Entering edit mode

That step is checking whether or not all the paired reads are two reads (left and right). elementNROWS>3 means there are some of the pairs contain three or more reads in one pair.

Just in case there are anything unexpected, could you please share me your reads filtering codes and a small piece of bam which can repeat your error?

Jianhong.

ADD REPLY
0
Entering edit mode

Hi Jianhong,

How can I send you these files and code? Can please provide me with an email address so that I could send you a link to the file and code?

ADD REPLY
0
Entering edit mode
ADD REPLY
0
Entering edit mode

Hi Jianhong,

Here is the filter code: http://s3.amazonaws.com/mudaser/bamfiles/ATACseq_heatmap.R And here is subset of the bam file: http://s3.amazonaws.com/mudaser/bamfiles/chr21-filt.bam

ADD REPLY
0
Entering edit mode

Hi Jianhong,

Here is the filter code: http://s3.amazonaws.com/mudaser/bamfiles/ATACseq_heatmap.R And here is subset of the bam file: http://s3.amazonaws.com/mudaser/bamfiles/chr21-filt.bam

ADD REPLY
0
Entering edit mode

Hi,

Thank you for sharing your bam file with me. Please have a try to change the line

gal <- readBamFile(bamfile, tag=tags, asMates=TRUE, bigFile=TRUE)

to

gal <- readBamFile(bamfile, tag=tags, asMates=TRUE, bigFile=TRUE, flag = scanBamFlag(isSecondaryAlignment = FALSE, isUnmappedQuery = FALSE, isNotPassingQualityControls = FALSE, isSupplementaryAlignment = FALSE))

Because current ATACseqQC can not handle supplementary alignments.

Jianhong.

ADD REPLY
0
Entering edit mode

Thank you for your valuable insight. The test worked. Best regards!

ADD REPLY

Login before adding your answer.

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