Attention:
The frontline support team will be unavailable to answer questions on April 15th and 17th 2019. We will be back soon after. Thank you for your patience and we apologize for any inconvenience!

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.