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
Release history Release notifications | RSS feed
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)
Built Distribution
seqalign-0.1.15-py3-none-any.whl
(14.5 kB
view hashes)
Close
Hashes for seqalign-0.1.15-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 8fcbaa5d26d92a23728f11dc9fe05a987913740bc17d0908929826a317d34dce |
|
MD5 | d30c23ecaca7ead440044fba5d107373 |
|
BLAKE2b-256 | 2f5061490c3905f4fd2a11cc58171565ae1bfaec5dd662ead296090ce76e9982 |