Hi all, anyone can help with the following from the ATACseqQC package:
objs <- splitBam(bamfile, tags=tags, outPath=outPath,
+ txs=txs, genome=genome,
+ conservation=phastCons100way.UCSC.hg19)
ERROR as follows:
[E::sam_parse1] unrecognized type N
[W::sam_read1] Parse error at line 26
[E::sam_parse1] unrecognized type N
[W::sam_read1] Parse error at line 26
Error in shiftGAlignmentsList(chunk0, positive = positive,negative = negative) : all(elementNROWS(gal) < 3) is not TRUE
this happens while I am trying to split the bam file into nucleosomefree, mono/di/trinucleosome areas
thank you beforehand
sorry for replying after so long I am sending you an email with the bam file
by the way the script I used was the one provided I changed nothing in order to see if it works in my case. Session info:
Hi, Have a try with newest version of ATACseqQC. The error should be handled by the program.
Jianhong.
Hi, I am getting the same error with mm10 with v1.8.5 downloaded 8/2/19. Is there any resolution or recommendation for this? Thanks!
Edit - it is giving the same error with 'splitBam' function in v1.9.7 downloaded from github. However, this error has disappeared when using shiftGAlignmentsList in v1.9.7. With v1.8.5 from Bioconductor, both functions gave the error.
you can update to the development version. Or set flag=scanBamFlag(isSecondaryAlignment = FALSE, isUnmappedQuery=FALSE, isNotPassingQualityControls = FALSE, isSupplementaryAlignment = FALSE) for readBamFile in released version to overcome the issue of all(elementNROWS(gal) < 3) is not TRUE
see https://support.bioconductor.org/p/122819/#123135
From: "ian.beddows [bioc]" noreply@bioconductor.org Reply-To: "reply+9ee280d2+code@bioconductor.org" <reply+9ee280d2+code@bioconductor.org> Date: Friday, August 2, 2019 at 9:40 AM To: "Jianhong Ou, Ph.D." jianhong.ou@duke.edu Subject: [bioc] C: atacseqQC error while splittin the bamfile
Activity on a post you are following on support.bioconductor.orghttps://urldefense.proofpoint.com/v2/url?u=https-3A__support.bioconductor.org&d=DwMCaQ&c=imBPVzF25OnBgGmVOlcsiEgHoG1i6YHLR0Sj_gZ4adc&r=PXg851DHXyo-Gs3eMIfeo49gUXVh-JSZu_MZDDxGun8&m=VtppQj4W1krIwaMiKKEZyV4S7oN-rHihMqROuu_fRio&s=QV61Y4n6aFg1XovygloDRAUEr9nemJtoQd7Pbin2nyw&e=
User ian.beddowshttps://urldefense.proofpoint.com/v2/url?u=https-3A__support.bioconductor.org_u_21527_&d=DwMCaQ&c=imBPVzF25OnBgGmVOlcsiEgHoG1i6YHLR0Sj_gZ4adc&r=PXg851DHXyo-Gs3eMIfeo49gUXVh-JSZu_MZDDxGun8&m=VtppQj4W1krIwaMiKKEZyV4S7oN-rHihMqROuu_fRio&s=h2KjbO0EffKD6uFncnz6W2MFhn9u-6oR34jtnO97CWQ&e= wrote Comment: atacseqQC error while splittin the bamfilehttps://urldefense.proofpoint.com/v2/url?u=https-3A__support.bioconductor.org_p_120707_-23123438&d=DwMCaQ&c=imBPVzF25OnBgGmVOlcsiEgHoG1i6YHLR0Sj_gZ4adc&r=PXg851DHXyo-Gs3eMIfeo49gUXVh-JSZu_MZDDxGun8&m=VtppQj4W1krIwaMiKKEZyV4S7oN-rHihMqROuu_fRio&s=C1XzjEp81u7RUI9Y8ODcLb3409zIpBPkXXD3gQZLGJg&e=:
Hi, I am getting the same error with mm10 with v1.8.5 downloaded 8/2/19. Is there any resolution or recommendation for this? Thanks!