Bug Bulletin: we have identified a bug that affects indexing when producing gzipped VCFs. This will be fixed in the upcoming 3.2 release; in the meantime you need to reindex gzipped VCFs using Tabix.

Compare Cohort 1 and Cohort 2: Make file with SNPs that are only in Cohort 1, not in Cohort 2

myoglumyoglu Posts: 39Member

Hi! The title says it all. Im trying to narrow down my vcf. Which tools is good to answer this question?

I have tried this:

java -Xmx4g -jar GenomeAnalysisTK.jar \ -R genome.fa \ --filter_reads_with_N_cigar \ -T CombineVariants \ -V:COHORT1 path \ -V:COHORT2 path \ --filteredAreUncalled \ -o path

But the output seams incomprehensibly massive!

Thanks!

Tagged:

Best Answer

Sign In or Register to comment.