This is a test version of Biostars. For the public version, visit https://www.biostars.org.
varscan not taking opstions

I am running varscan 2.4.3, mpileup2snp is not taking the settings as you can below, it works but with default settings any suggestion how to change them?

varscan mpileup2snp  sam1 –-min-coverage 10 –-min-var-freq 0.05 –-p-value 0.05 --output-vcf 1 > test6var005.vcf
Only SNPs will be reported
Warning: No p-value threshold provided, so p-values will not be calculated
Min coverage:   8
Min reads2: 2
Min var freq:   0.2
Min avg qual:   15
P-value thresh: 0.01
Reading input from sam1

thanks

Marco

varscan2 mpileup2snp setting
• 1,649 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by fabbri.marco • fabbri.marco

1 answer

varscan is a java program. It must be invoked with something like: 

USAGE: java -jar VarScan.jar  [COMMAND] [OPTIONS]

I imagine, you're using a bash wrapper around the command above. Don't use this wrapper and invoke varscan directly with java.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this