We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

Filter in DepthOfCoverage 3.3-0

HerderHerder NetherlandsMember

Hello GATK,

We've been using the DepthOfCoverage tool in our pipeline and noticed that the =>Q30 score per gene provided by GATK is significantly different than the =>Q30 score we calculate ourselves without any filters. On your manual page for DOC, you mention using several filters (although it isn't specified if this applies to all versions of GATK or just the last one):
NotPrimaryAlignmentFilter
FailsVendorQualityCheckFilter
DuplicateReadFilter
UnmappedReadFilter
MalformedReadFilter
BadCigarFilter

In your general manual page, you mention there is a parameter for disabling filter (--disable_read_filter, -drf). This doesn't seem to be implemented yet in 3.3-0, seeing as I get an error that the argument isn't defined and the argument isn't showing up in the help page either. Is there another way to disable filters in 3.3-0?

Kind regards,
Mark

Tagged:

Answers

  • SheilaSheila Broad InstituteMember, Broadie ✭✭✭✭✭

    @Herder
    Hi Mark,

    Can you post the exact command you ran and the exact log output with the error?

    Thanks,
    Sheila

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    The --disable_read_filter functionality is only enabled for a subset of read filters (at the moment, just DuplicateReadFilter). This is done by design as the other filters are all necessary for the smooth operation of the tools. If you want to do an apples-to-apples comparison you can run PrintReads with the same filters to generate a file that you can use for your calculations.

Sign In or Register to comment.