Skip to content

Commit

Permalink
remove bcftools sample flag
Browse files Browse the repository at this point in the history
  • Loading branch information
jaleezyy committed Jun 6, 2024
1 parent 4282516 commit a426b3c
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions Snakefile
Original file line number Diff line number Diff line change
Expand Up @@ -647,9 +647,9 @@ rule run_freebayes:
done
# apply ambiguous variants first using IUPAC codes. this variant set cannot contain indels or the subsequent step will break
bcftools consensus -s '-' -f {input.reference} -I {params.out}.ambiguous.norm.vcf.gz > {params.out}.ambiguous.fasta
bcftools consensus -f {input.reference} -I {params.out}.ambiguous.norm.vcf.gz > {params.out}.ambiguous.fasta
# apply remaninng variants, including indels
bcftools consensus -s '-' -f {params.out}.ambiguous.fasta -m {params.out}.mask.txt {params.out}.fixed.norm.vcf.gz | sed s/MN908947\\.3.*/{wildcards.sn}/ > {output.consensus}
bcftools consensus -f {params.out}.ambiguous.fasta -m {params.out}.mask.txt {params.out}.fixed.norm.vcf.gz | sed s/MN908947\\.3.*/{wildcards.sn}/ > {output.consensus}
"""

rule consensus_compare:
Expand Down

0 comments on commit a426b3c

Please sign in to comment.