Open bernt-matthias opened 2 years ago
I run lofreq call-parallel with different numbers for --pp-threads and get different vcf output:
lofreq call-parallel
--pp-threads
lofreq call-parallel --pp-threads 1 --verbose --ref 'reference.fa' --out variants.vcf --min-cov 1 --max-depth 1000000 --min-bq 6 --min-alt-bq 6 --min-mq 0 --max-mq 255 --min-jq 0 --min-alt-jq 0 --def-alt-jq 0 --sig 0.01 --bonf dynamic reads.bam
##FILTER=<ID=min_snvqual_33,Description="Minimum SNV Quality (Phred) 33"> ##FILTER=<ID=min_snvqual_38,Description="Minimum SNV Quality (Phred) 38">
lofreq call-parallel --pp-threads 2 --verbose --ref 'reference.fa' --out variants.vcf --min-cov 1 --max-depth 1000000 --min-bq 6 --min-alt-bq 6 --min-mq 0 --max-mq 255 --min-jq 0 --min-alt-jq 0 --def-alt-jq 0 --sig 0.01 --bonf dynamic reads.bam
##FILTER=<ID=min_snvqual_38,Description="Minimum SNV Quality (Phred) 38"> ##FILTER=<ID=min_snvqual_38,Description="Minimum SNV Quality (Phred) 38">
Wondering if this is expected/desired.
I run
lofreq call-parallelwith different numbers for--pp-threadsand get different vcf output:lofreq call-parallel --pp-threads 1 --verbose --ref 'reference.fa' --out variants.vcf --min-cov 1 --max-depth 1000000 --min-bq 6 --min-alt-bq 6 --min-mq 0 --max-mq 255 --min-jq 0 --min-alt-jq 0 --def-alt-jq 0 --sig 0.01 --bonf dynamic reads.bamlofreq call-parallel --pp-threads 2 --verbose --ref 'reference.fa' --out variants.vcf --min-cov 1 --max-depth 1000000 --min-bq 6 --min-alt-bq 6 --min-mq 0 --max-mq 255 --min-jq 0 --min-alt-jq 0 --def-alt-jq 0 --sig 0.01 --bonf dynamic reads.bamWondering if this is expected/desired.