Skip to main content

Manage sequence alignments

Project description

seqalign

Manage sequence alignments

Installation

pip install seqalign

or

pip install --user seqalign

If you have not already installed picardtools, you will also need to run:

picardtools-download

Examples

with SequenceAlignment(<path to input BAM or FASTQ file>) as sa:
    sa.cleans_up_bam = False
    sa.remove_supplementary_alignments()
    sa.samtools_sort(memory_limit=10)
    sa.samtools_index()
    sa.write(<path to output BAM file>)

Notes

The input_file argument to SequenceAlignment() should be a string for single-end reads or for data that is already aligned. For raw paired-end reads, it should be a tuple containing two strings giving the paths to the two FASTA / FASTQ files.

Command line usage

seqalign includes a couple of pre-built pipelines for aligning ChIP-seq and ATAC-seq data. They rely on certain environment variables being set, as described here To see the basic command line documentation, simply enter seqalign at the terminal:

seqalign
There are two commands for aligning sequencing data:
`seqalign-chip-seq` and `seqalign-atac-seq`. See their
respective documentation by running:

seqalign-chip-seq --help
seqalign-atac-seq --help

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

seqalign-0.1.15.tar.gz (12.7 kB view hashes)

Uploaded Source

Built Distribution

seqalign-0.1.15-py3-none-any.whl (14.5 kB view hashes)

Uploaded Python 3

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page