We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 34161cb commit 47c0706Copy full SHA for 47c0706
SequenceAnalysis/resources/external/viral_consensus.sh
@@ -23,7 +23,7 @@ OUT_IUPAC=`basename $BAM .bam`".consensus.iupac.fasta"
23
24
# call variants
25
echo 'Calling variants'
26
-$BCFTOOLS mpileup -Ou -d 20000 -f $FASTA $BAM | $BCFTOOLS call -c --ploidy 1 --threads $THREADS -mv -Oz -o $VCF_CALLS
+$BCFTOOLS mpileup -Ou -d 20000 -f $FASTA $BAM | $BCFTOOLS call -c --ploidy 1 --threads $THREADS -v -Oz -o $VCF_CALLS
27
$BCFTOOLS index -t $VCF_CALLS
28
COUNT1=`$BCFTOOLS view -H $VCF_CALLS | wc -l`
29
echo 'Variants called: '$COUNT1
0 commit comments