.. _`bio/bwa/samse`:

BWA SAMSE
=========


.. image:: https://img.shields.io/github/issues-pr/snakemake/snakemake-wrappers/bio/bwa/samse?label=version%20update%20pull%20requests
   :target: https://github.com/snakemake/snakemake-wrappers/pulls?q=is%3Apr+is%3Aopen+label%3Abio/bwa/samse

Map single-end reads with bwa samse. For more information about BWA see `BWA documentation <http://bio-bwa.sourceforge.net/bwa.shtml>`_.


Example
-------

This wrapper can be used in the following way:

.. code-block:: python

    rule bwa_samse:
        input:
            fastq="reads/{sample}.1.fastq",
            sai="sai/{sample}.1.sai",
            idx=multiext("genome", ".amb", ".ann", ".bwt", ".pac", ".sa"),
        output:
            "mapped/{sample}.bam",
        params:
            extra=r"-r '@RG\tID:{sample}\tSM:{sample}'",  # optional: Extra parameters for bwa.
            sort="none",  # optional: Enable sorting. Possible values: 'none', 'samtools' or 'picard'`
            sort_order="queryname",  # optional: Sort by 'queryname' or 'coordinate'
            sort_extra="",  # optional: extra arguments for samtools/picard
        log:
            "logs/bwa_samse/{sample}.log",
        wrapper:
            "v3.0.1/bio/bwa/samse"

Note that input, output and log file paths can be chosen freely.

When running with

.. code-block:: bash

    snakemake --use-conda

the software dependencies will be automatically deployed into an isolated environment before execution.


Software dependencies
---------------------

* ``bwa=0.7.17``
* ``samtools=1.18``
* ``picard-slim=3.1.1``




Authors
-------

* Julian de Ruiter


Code
----

.. code-block:: python

    """Snakemake wrapper for bwa sampe."""

    __author__ = "Julian de Ruiter"
    __copyright__ = "Copyright 2017, Julian de Ruiter"
    __email__ = "julianderuiter@gmail.com"
    __license__ = "MIT"


    from os import path

    from snakemake.shell import shell

    index = snakemake.input.get("idx", "")
    if isinstance(index, str):
        index = path.splitext(snakemake.input.idx)[0]
    else:
        index = path.splitext(snakemake.input.idx[0])[0]

    # Extract arguments.
    extra = snakemake.params.get("extra", "")

    sort = snakemake.params.get("sort", "none")
    sort_order = snakemake.params.get("sort_order", "coordinate")
    sort_extra = snakemake.params.get("sort_extra", "")

    log = snakemake.log_fmt_shell(stdout=False, stderr=True)

    # Determine which pipe command to use for converting to bam or sorting.
    if sort == "none":
        # Simply convert to bam using samtools view.
        pipe_cmd = "samtools view -Sbh -o {snakemake.output[0]} -"

    elif sort == "samtools":
        # Sort alignments using samtools sort.
        pipe_cmd = "samtools sort {sort_extra} -o {snakemake.output[0]} -"

        # Add name flag if needed.
        if sort_order == "queryname":
            sort_extra += " -n"

        # Use prefix for temp.
        prefix = path.splitext(snakemake.output[0])[0]
        sort_extra += " -T " + prefix + ".tmp"

    elif sort == "picard":
        # Sort alignments using picard SortSam.
        pipe_cmd = (
            "picard SortSam {sort_extra} INPUT=/dev/stdin"
            " OUTPUT={snakemake.output[0]} SORT_ORDER={sort_order}"
        )

    else:
        raise ValueError("Unexpected value for params.sort ({})".format(sort))

    # Run command.
    shell(
        "(bwa samse"
        " {extra}"
        " {index}"
        " {snakemake.input.sai}"
        " {snakemake.input.fastq}"
        " | " + pipe_cmd + ") {log}"
    )


.. |nl| raw:: html

   <br>