Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.We appreciate your help!

Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

SelectVariants

Hi, I'm using SelectVariants to process VCF file downloaded from ftp-mouse.sanger.ac.uk/current_snps/. This vcf contains SNPs for many mouse strains, and I need to extract information for only two strains for later use of ASEReadCounter. This is my command:
java -Xmx4g -jar GenomeAnalysisTK.jar -T SelectVariants -R GRCm38_68.fa -V mgp.v5.merged.snps_all.dbSNP142.vcf -o mgp.v5.merged.snps_all.dbSNP142_subset.vcf -sn CAST_EiJ -sn 129S1_SvImJ -selectType SNP -restrictAllelesTo BIALLELIC

First lines in logs look ok, but they are followed by a number of "@^@^@^@^@^@^@^":

INFO 11:43:44,646 HelpFormatter - ---------------------------------------------------------------------------------
INFO 11:43:44,650 HelpFormatter - The Genome Analysis Toolkit (GATK) v3.7-0-gcfedb67, Compiled 2016/12/12 11:21:18
INFO 11:43:44,651 HelpFormatter - Copyright (c) 2010-2016 The Broad Institute
INFO 11:43:44,651 HelpFormatter - For support and documentation go to https://software.broadinstitute.org/gatk
INFO 11:43:44,651 HelpFormatter - [Tue Feb 28 11:43:44 EST 2017] Executing on Linux 2.6.32-642.3.1.el6.x86_64 amd64
INFO 11:43:44,652 HelpFormatter - Java HotSpot(TM) 64-Bit Server VM 1.8.0_45-b14
INFO 11:43:44,657 HelpFormatter - Program Args: -T SelectVariants -R GRCm38_68.fa -V mgp.v5.merged.snps_all.dbSNP142.vcf -o mgp.v5.merged.snps_all.dbSNP142_subset.vcf -sn CAST_EiJ -sn 129S1_SvImJ -selectType SNP -restrictAllelesTo BIALLELIC
INFO 11:43:44,668 HelpFormatter - Executing as [email protected] on Linux 2.6.32-642.3.1.el6.x86_64 amd64; Java HotSpot(TM) 64-Bit Server VM 1.8.0_45-b14.
INFO 11:43:44,669 HelpFormatter - Date/Time: 2017/02/28 11:43:44
INFO 11:43:44,669 HelpFormatter - ---------------------------------------------------------------------------------
INFO 11:43:44,669 HelpFormatter - ---------------------------------------------------------------------------------
INFO 11:43:44,696 GenomeAnalysisEngine - Strictness is SILENT
INFO 11:43:44,833 GenomeAnalysisEngine - Downsampling Settings: Method: BY_SAMPLE, Target Coverage: 1000

And I don't get any output. What am I doing wrong? Thank you.

Tagged:

Answers

Sign In or Register to comment.