Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Rsamtools version not compatible with VN:1.4 BAM files? #34

Open
ToonRosseel opened this issue Oct 18, 2021 · 0 comments
Open

Rsamtools version not compatible with VN:1.4 BAM files? #34

ToonRosseel opened this issue Oct 18, 2021 · 0 comments

Comments

@ToonRosseel
Copy link

ToonRosseel commented Oct 18, 2021

Dear, I would really like to test DECoN. The installation was successful, but I get following error message when I try to run ReadInBams.R with our BAM files:

Error in value[[3L]](cond) : 
  failed to open BamFile: SAM/BAM header missing or empty
  file: 'file.bam'
Calls: getBamCounts ... tryCatch -> tryCatchList -> tryCatchOne -> <Anonymous>
In addition: Warning message:
In doTryCatch(return(expr), name, parentenv, handler) :
  [bam_header_read] invalid BAM binary header (this is not a BAM file).
Execution halted

According the error message, the headers are missing, but they are not! The upper header lines of my BAM files look like this:

@HD	VN:1.4	SO:coordinate
@SQ	SN:chr1	LN:248956422
@SQ	SN:chr2	LN:242193529
@SQ	SN:chr3	LN:198295559
@SQ	SN:chr4	LN:190214555
@SQ	SN:chr5	LN:181538259
@SQ	SN:chr6	LN:170805979
@SQ	SN:chr7	LN:159345973
@SQ	SN:chr8	LN:145138636
...

I seems to me that the error message is made by the Rsamtools package version 1.18.3. If I read my BAM file with a more recent Rsamtools version (eg. v2.8.0), I do not get this message. But as DECoN depends on an old R version (3.1.2), it's not just possibly to use a more recent Rsamtools version (and all it's dependencies). Does anyone know an easy (or less easy) fix for this issue?

@ToonRosseel ToonRosseel changed the title Rsamtools version not compatible with VN:1.4 Bam files? Rsamtools version not compatible with VN:1.4 BAM files? Oct 18, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant