Exome SAMtools Workflow

public public 1yr ago Version: Version 1 0 bookmarks

Exome SAMtools Workflow

Code Snippets

34
35
36
37
38
39
40
41
42
43
44
45
46
"""
echo "${meta.id}" > sample_name.list
bcftools mpileup \\
    --fasta-ref $fasta \\
    $options.args \\
    $bam \\
    | bcftools call --output-type v $options.args2 \\
    | bcftools reheader --samples sample_name.list \\
    | bcftools view --output-file ${prefix}.vcf.gz --output-type z $options.args3
tabix -p vcf -f ${prefix}.vcf.gz
bcftools stats ${prefix}.vcf.gz > ${prefix}.bcftools_stats.txt
echo \$(bcftools --version 2>&1) | sed 's/^.*bcftools //; s/ .*\$//' > ${software}.version.txt
"""

Login to post a comment if you would like to share your experience with this workflow.

Do you know this workflow well? If so, you can request seller status , and start supporting this workflow.

Free

Created: 1yr ago
Updated: 1yr ago
Maitainers: public
URL: https://github.com/inab/ipc_workflows/tree/main/exome/samtools
Name: exome-samtools
Version: Version 1
Badge:
workflow icon

Insert copied code into your website to add a link to this workflow.

Downloaded: 0
Copyright: Public Domain
License: Boost Software License 1.0
  • Future updates

Related Workflows

cellranger-snakemake-gke
snakemake workflow to run cellranger on a given bucket using gke.
A Snakemake workflow for running cellranger on a given bucket using Google Kubernetes Engine. The usage of this workflow ...