• Total activity 64
  • Last activity
  • Member since
  • Following 0 users
  • Followed by 0 users
  • Votes 0
  • Subscriptions 21


Recent activity by vctrymao Sort by recent activity Recent activity Votes
  • MergeBamAlignment produces unknown error


    I am using MergeBamAlignment:  ``` java -Dsamjdk.compression_level=5 -Xms3000m -Xmx3200m -jar /n/data1/hms/dbmi/park/alon/software/picard.jar \MergeBamAlignment \VALIDATION_STRINGENCY=SILENT \EXPEC...

  • Funcotator table format

    Is it possible to use Funcotator for germline calls but format the output in a table format instead of a VCF? I understand that MAF isn't yet available for germline Funcotator, but is there any oth...

  • Getting HaplotypeCaller SNP support in tumor

    I have a list of germline SNPs from HaplotypeCaller. I want to get the read level information, etc. of these SNPs in the tumor BAM file and format it into a VCF (ie. to see the read level support f...

  • HaplotypeCaller and low-VAF variants

    I am using HaplotypeCaller to generate germline variants. I noticed that HaplotypeCaller calls some variants with very low VAF values (see below):  What is the reasoning behind this? In what situa...

  • Missing mutations

    I am using the latest version of Mutect2 ( and trying to call mutations from Mutect. I noticed that Mutect is not calling muations with clear read support in IGV. I tried some suggestions f...

  • Combining libraries for sample


    I am trying to follow this guide in order to combine libraries for a single sample for alignment. However, the article seems to focus primarily on combining lanes. How should I combine multiple lib...

  • Using GenomicsDB with intervals


    Hi, I am using GATK and trying to run GenomicsDB with some gVCFs I have. I noticed that it requires an interval file; however, I am interested in running this across the whole genome, and t...

  • MarkDuplicates returns error with Multilane samples?

    I am using GATK I am following this guide (https://gatk.broadinstitute.org/hc/en-us/articles/360035889471-How-should-I-pre-process-data-from-multiplexed-sequencing-and-multi-library-design...

  • Mutect tumor-only vs matched normal calls?


    I am using GATK This version includes the new FilterMutectCalls module.  I was wondering, is there anything different that Mutect2 does in tumor-only mode vs. in matched-normal mode? In th...

  • Germline mutations in Mutect2 vs HaplotypeCaller

    I am using GATK I was wondering, is it usual to observe variants in the Mutect2 output that are labeled PASS but are still found by HaplotypeCaller? And is the reverse true: is it usual to...