To celebrate the release of GATK 4.0, we are giving away free credits for running the GATK4 Best Practices pipelines in FireCloud, our secure online analysis portal. It’s first come first serve, so sign up now to claim your free credits worth $250. Sponsored by Google Cloud. Learn more at https://software.broadinstitute.org/firecloud/documentation/freecredits

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 sv111@clarinet002-069 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.