Currently the "Log in with Google" button redirects you to a "Page not found." This is an issue that our forum vendors are working on fixing. In the meantime, while on the "Page not found" you can edit the URL to delete the second gatk, firecloud, or wdl (depending on what subforum you are acessing).
SelectVariants quick question about concordance/discordance
This is my code:
java -Xmx4g -jar GenomeAnalysisTK-2.7-2-g6bda569/GenomeAnalysisTK.jar \ -R genome.fa \ --filter_reads_with_N_cigar \ -T CombineVariants \ -V:NORMALSpost RNA-edit_SNPs_NORMALS-post_SNPeff_VA.vfc \ -V:OMNI 1000G_omni2.5.hg19.vcf \ -V:db137 dbsnp_137.hg19.vcf \ -V:Hapmap hapmap_3.3.hg19.vcf \ -V:ESP1,VCF ESP6500SI-V2-SSA137.updatedRsIds.Allchr.snps_indels_FIX.vcf \ --out NORMALS-post_TruePositives_raw.vcf && java -Xmx4g -jar GenomeAnalysisTK.jar \ -R genome.fa \ --filter_reads_with_N_cigar \ -T SelectVariants \ -V NORMALS-post_TruePositives_raw.vcf \ -select "set == 'NORMALSpost'" \ --out NORMALS-post_TruePositives.vcf
I hope I got this right: I now have SNPs (I have only called SNPs, not indels) that is found in my sample, 1000G, dbSNP and ESP. Right?
Can I easily make a VCF with all SNPs in MY sample only, not found in anything else by tweaking this code? Or must I do it some way else?