The Frontline Support team will be slow to respond December 17-18 due to an institute-wide retreat and offline December 22- January 1, while the institute is closed. Thank you for your patience during these next few weeks. Happy Holidays!
Bug in HaplotypeCaller: Lines meant to STDERR go to STDOUT
I use HaplotypeCaller v. 3.6-0-g89b7209. I want to compress the output before writing it to the disk and use command like this:
gatk -T HaplotypeCaller ... -o /dev/stdout | grep -v -e --- -e WARN | bgzip -s - > output.gvcf.gz
The grep command is needed to get around a bug that sends lines meant to STDERR to STDOUT. This is similar to a bug I reported in v. 3.3: