We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
https://github.com/OpenGene/ctdna-pipeline/blob/master/pipeline.sh#L24 says:
java -jar VarScan.v2.3.8.jar mpileup2snp outdir/test.dedup.mpileup --min-coverage 4 --min-reads2 2 --min-avg-qual 20 --min-var-freq 0.001 --min-freq-for-hom 90 --output-vcf 1 >outdir/test.snp.vcf
Now I'm not sure whether the value of --min-freq-for-hom 90 is OK. Because now the new version of VarScan (v2.4.3) says:
--min-freq-for-hom 90
USAGE: java -jar VarScan.jar mpileup2cns [pileup file] OPTIONS mpileup file - The SAMtools mpileup file OPTIONS: ... --min-freq-for-hom Minimum frequency to call homozygote [0.75] ...
Note, the default value is 0.75 rather than 75 (%). Therefore perhaps we should use value 0.9 rather than 90 (5)
The text was updated successfully, but these errors were encountered:
No branches or pull requests
https://github.com/OpenGene/ctdna-pipeline/blob/master/pipeline.sh#L24 says:
java -jar VarScan.v2.3.8.jar mpileup2snp outdir/test.dedup.mpileup --min-coverage 4 --min-reads2 2 --min-avg-qual 20 --min-var-freq 0.001 --min-freq-for-hom 90 --output-vcf 1 >outdir/test.snp.vcf
Now I'm not sure whether the value of
--min-freq-for-hom 90
is OK. Because now the new version of VarScan (v2.4.3) says:Note, the default value is 0.75 rather than 75 (%). Therefore perhaps we should use value 0.9 rather than 90 (5)
The text was updated successfully, but these errors were encountered: