David Benjamin
- Total activity 560
- Last activity
- Member since
- Following 0 users
- Followed by 1 user
- Votes 2
- Subscriptions 197
Comments
Recent activity by David Benjamin Sort by recent activity-
Regarding 2), restricting your intervals to just the force-calling alleles (i.e. -L ${GGA_path}/vcf/${sample}_GGA_union.vcf.gz) might cause problems. Can you try applying interval padding (I would...
-
@Fia It is totally normal for FilterMutectCalls to filter 90% or more of Mutect2 variant calls. Mutect2 is designed to be very permissive and naive, leaving almost all the responsibility for filte...
-
When a sample has no read coverage near a force-calling allele Mutect2 can't emit it because the GATK engine skips over the region. Is that the case for your SNPs that were not called? Some discre...
-
M Carta I'm not sure this is the root of the issue but the lack of sequence dictionary lines eg ##contig=<ID=chr2,length=242193529> in the header of your "manual" VCFs worries me a bit. How about...
-
Melcar Collodetti A few comments on the workflow, which I hope will be helpful: The -L argument in GetPileupSummaries should usually be the same vcf file as the -V argument. We recommend using the...
-
Probably 1-4 days for code review and merging into the master branch, followed by roughly a month between minor releases.
-
Luka Culibrk This new PR: https://github.com/broadinstitute/gatk/pull/8717 removes germline events from the accounting of the ECNT annotation and should fix the unfortunate interaction between the...
-
Jyoti Mridha This looks like a custom-made filtering pipeline. It does not resemble our best practices. The hard filtering most likely does very little since it is filtering based on several anno...
-
Yes, that will work fine.
-
Also, here are the standard public service announcements regarding homemade panels of normals and tumor-only calling 1) Using the panels of normals in our public google bucket gs://gatk-best-practi...