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.