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

BRASS filter generation fails with an uninformative error when *.bas not present #71

Closed
edawson opened this issue Jul 25, 2018 · 1 comment

Comments

@edawson
Copy link
Contributor

edawson commented Jul 25, 2018

I was working through generating a PON and got a strange error when running the brassI_np_in.pl script:

> brassI_np_in.pl `pwd` 1 /dat/redacted.20k.bam
Starting: /opt/wtsi-cgp/bin/bamcollate2 outputformat=sam exclude=PROPER_PAIR,UNMAP,MUNMAP,SECONDARY,QCFAIL,DUP,SUPPLEMENTARY mapqthres=6 classes=F,F2 T=/home/ubuntu/redacted/tmpMap/bamcollate2 filename=/dat/redacted.20k.bam | /usr/bin/perl /opt/wtsi-cgp/bin/brassI_prep_bam.pl -b /dat/redacted.20k.bam.bas -np | /opt/wtsi-cgp/bin/bamsort inputformat=sam verbose=0 index=1 md5=1 tmpfile=/home/ubuntu/redacted/tmpMap/bamsort md5filename=/home/ubuntu/redacted/redacted.brm.bam.md5 indexfilename=/home/ubuntu/redacted/redacted.brm.bam.bai O=/home/ubuntu/redacted/redacted.brm.bam
ScramDecoder::readAlignment(): failed to read alignment without reaching EOF

/opt/wtsi-cgp/bin/../lib/libmaus2.so.2(libmaus2::util::StackTrace::StackTrace()+0x54)[0x7f42b94b2a54]
/opt/wtsi-cgp/bin/bamsort(libmaus2::exception::LibMausException::LibMausException()+0x20)[0x454550]
/opt/wtsi-cgp/bin/bamsort(libmaus2::bambam::ScramDecoder::readAlignmentInternal(bool)+0x2b3)[0x47b7b3]
/opt/wtsi-cgp/bin/bamsort(bamsort(libmaus2::util::ArgInfo const&)+0x2519)[0x44e309]
/opt/wtsi-cgp/bin/bamsort(main+0x177e)[0x44724e]
/lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5)[0x7f42b7da2f45]
/opt/wtsi-cgp/bin/bamsort()[0x44819f]

WARNING: SAM header designates more than one PG tree root by PP tags.
WARNING: PG line with id bwamem.2 has multiple children referencing it by PP tags.
WARNING: PG line with id bwamem.3 has multiple children referencing it by PP tags.
WARNING: PG line with id bwamem.4 has multiple children referencing it by PP tags.
WARNING: PG line with id bwamem.8 has multiple children referencing it by PP tags.

I grabbed the command line invocations from the first line, started running them individually and found that the second command fails because there's no .bas file:

/opt/wtsi-cgp/bin/bamcollate2 outputformat=sam exclude=PROPER_PAIR,UNMAP,MUNMAP,SECONDARY,QCFAIL,DUP,SUPPLEMENTARY mapqthres=6 classes=F,F2 T=/home/ubuntu/redacted/tmpMap/bamcollate2 filename=/dat/redacted.20k.bam | /usr/bin/perl /opt/wtsi-cgp/bin/brassI_prep_bam.pl -b /dat/redacted.20k.bam.bas -np
File not found: /dat/redacted.20k.bam.bas

Creating a .bas file fixes the issue.

Is it possible to propagate the underlying error better?

@sb43
Copy link
Member

sb43 commented Oct 19, 2018

This is now fixed.

@sb43 sb43 closed this as completed Nov 2, 2018
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

2 participants