align BS-Seq reads with dragmap
Project description
# dragmap_meth (dragmap_meth.py)
Alignment of BS-Seq reads using [dragmap](https://github.com/Illumina/DRAGMAP).
## Intro
This works for single-end reads and for paired-end reads from the directional protocol (most common).
## Installation
### Install dragmap-meth via pip
dragmap-meth is available on https://pypi.org/project/dragmap-meth/
` pip install dragmap-meth `
## Install dragmap-meth via conda
` conda env create -n env4dragmap-meth --file environment.yaml python=3 git clone https://github.com/xie186/dragmap-meth.git cd dragmap-meth/ `
## QuickStart
The commands:
` python dragmap-meth.py buildhashtable -r ref.fa -o ref/ python dragmap-meth.py dragmap -ht ref/ -r1 t_R1.fastq.gz -r2 t_R2.fastq.gz |samtools view -bS - -o dragmap-meth.bam `
will create dragmap-meth.bam. To align single end-reads, specify only 1 file: -r1 some_read.fastq.gz
### Dependencies
` Python 3.10.1 dragmap `
## Reference:
Introducing DRAGMAP, the new genome mapper in DRAGEN-GATK: https://gatk.broadinstitute.org/hc/en-us/articles/4410953761563-Introducing-DRAGMAP-the-new-genome-mapper-in-DRAGEN-GATK
Demystifying the versions of GRCh38/hg38 Reference Genomes, how they are used in DRAGEN™ and their impact on accuracy: https://www.illumina.com/science/genomics-research/articles/dragen-demystifying-reference-genomes.html
https://github.com/DavidStreid/bio_docker/tree/main/dragmap
https://github.com/HudsonAlpha/CSL_public_benchmark
## Acknowledgement
Special thanks to [bwa-meth](https://github.com/brentp/bwa-meth). Part of the codes were adapted from bwa-meth.
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
Built Distribution
Hashes for dragmap_meth-0.1.2-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | c13a17c52aa1b9d449cb98db20822cbac1f2dfdb9bad57055303e5aa6da6810f |
|
MD5 | 4f3bd56b16bfee90c7817da1906a47be |
|
BLAKE2b-256 | f24adf8fb6b583e39a85ea93d52132b41534f27c749a48dca6c69a6647ba942a |