[med-svn] [Git][med-team/mapdamage][upstream] New upstream version 2.2.1+dfsg
Steffen Möller
gitlab at salsa.debian.org
Mon Jun 29 23:27:32 BST 2020
Steffen Möller pushed to branch upstream at Debian Med / mapdamage
Commits:
a9e55a3f by Steffen Moeller at 2020-06-30T00:19:26+02:00
New upstream version 2.2.1+dfsg
- - - - -
4 changed files:
- README.md
- bin/mapDamage
- mapdamage/version.py
- setup.py
Changes:
=====================================
README.md
=====================================
@@ -1,29 +1,33 @@
## mapDamage
-[](http://bioconda.github.io/recipes/mapdamage2/README.html) [](https://anaconda.org/bioconda/mapdamage2/files)   [](https://www.repostatus.org/#inactive)
+[](http://bioconda.github.io/recipes/mapdamage2/README.html) [](https://anaconda.org/bioconda/mapdamage2/files)
+  [](https://www.repostatus.org/#inactive)
#### `bioconda` installation
-- python3 version **2.2.0**
-```
-conda install -c bioconda mapdamage2=2.2.0
-```
+* python3 version **2.2.1**
-- python3 version **2.2.0** **with** R and 4 mandatory packages for the Bayesian inference:
+```
+conda install -c bioconda mapdamage2=2.2.1
```
-conda install -c bioconda mapdamage2=2.2.0=pyr36_1
+
+* python3 version **2.2.1** **with** R and 4 mandatory packages for the Bayesian inference:
+
+```
+conda install -c bioconda mapdamage2=2.2.1=pyr36_1
```
---
### Important
-- From version `2.2.0` the `master` branch is requiring **python3** as `python2` is not supported from 2020-01-01.
+* From version `2.2.1` the `master` branch is requiring **python3** as `python2` is not supported from 2020-01-01.
-- Users with versions dating prior to June 12 2013 please update. A nasty bug that caused the statistical part of `mapDamage` to use half of the data for estimation of the damage parameters, sorry for the inconvenience..
+* Users with versions dating prior to June 12 2013 please update. A nasty bug that caused the statistical part of `mapDamage` to use half of the data for estimation of the damage parameters, sorry for the inconvenience..
### Introduction
+
Complete documentation, instructions, examples, screenshots and FAQ are available at [this address](http://ginolhac.github.io/mapDamage/).
[mapDamage2](https://geogenetics.ku.dk/publications/mapdamage2.0/) is a computational framework written in **Python3** and **R**, which tracks and quantifies DNA damage patterns
@@ -31,31 +35,30 @@ among ancient DNA sequencing reads generated by Next-Generation Sequencing platf
`mapDamage` was developed at the [Centre for GeoGenetics](https://geogenetics.ku.dk/) by the [Orlando Group ](https://geogenetics.ku.dk/research_groups/palaeomix_group/).
-
### Citation
+
If you use this program, please cite the following publication:
Jónsson H, Ginolhac A, Schubert M, Johnson P, Orlando L.
[mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters.
_Bioinformatics_ 23rd April 2013. doi: 10.1093/bioinformatics/btt193](http://bioinformatics.oxfordjournals.org/content/early/2013/05/17/bioinformatics.btt193)
-
The original `mapDamage1` was published in the following article:
Ginolhac A, Rasmussen M, Gilbert MT, Willerslev E, Orlando L.
[mapDamage: testing for damage patterns in ancient DNA sequences. _Bioinformatics_ 2011 **27**(15):2153-5
http://bioinformatics.oxfordjournals.org/content/27/15/2153](http://bioinformatics.oxfordjournals.org/content/27/15/2153)
-
### Test the no-stats part and rescaling
you can test `mapDamage` by running:
-```
+```
cd mapDamage/mapdamage/
python3 mp_test.py
```
should return
-```
+
+```
Started with the command: /usr/local/bin/mapDamage -i tests/test.bam -r tests/fake1.fasta -d tests/results --no-stats
Reading from 'tests/test.bam'
Writing results to 'tests/results/'
@@ -69,7 +72,7 @@ Ran 2 tests in 3.357s
OK
```
-
### Contact
+
Please report bugs and suggest possible improvements to Aurélien Ginolhac, Mikkel Schubert or Hákon Jónsson by email:
ginolhac at gmail.com, mikkelsch at gmail.com or jonsson.hakon at gmail.com.
=====================================
bin/mapDamage
=====================================
@@ -139,20 +139,6 @@ def _check_mm_option():
raise SystemExit("The --mapdamage-modules path (path=%s) must contain the mapdamage module" % path_to_mm)
return path_to_mm
-def is_paired(alignfile):
- """
- check if an alignment file contains paired end data by
- looking at the first 100 reads
- """
-
- n = 0
- for read in alignfile:
- if read.is_paired:
- return True
- n += 1
- if n == 100:
- return False
-
def main():
start_time = time.time()
@@ -241,14 +227,6 @@ def main():
else:
in_bam = pysam.Samfile(options.filename)
- # exit SAM/BAM is paired (look at top 100 reads)
- if is_paired(in_bam):
- logger.error("alignment file must be single-ended")
- sys.stderr.write("alignment file must be single-ended\n")
- sys.exit(1)
- # reset SAM/BAM to the begining
- in_bam = pysam.Samfile(options.filename)
-
reflengths = dict(list(zip(in_bam.references, in_bam.lengths)))
# check if references in SAM/BAM are the same in the fasta reference file
fai_lengths = mapdamage.seq.read_fasta_index(options.ref + ".fai")
=====================================
mapdamage/version.py
=====================================
@@ -2,4 +2,4 @@
try:
from mapdamage._version import __version__
except ImportError:
- __version__ = "2.2.0"
+ __version__ = "2.2.1"
=====================================
setup.py
=====================================
@@ -35,7 +35,7 @@ def setup_version():
version = subprocess.check_output(("git", "describe", "--always", "--tags", "--dirty"))
with open(os.path.join("mapdamage", "_version.py"), "w") as handle:
handle.write("#!/usr/bin/env python\n")
- handle.write("__version__ = %r\n" % (version.strip(),))
+ handle.write("__version__ = %r\n" % (version.decode("utf-8").strip(),))
except (subprocess.CalledProcessError, OSError) as error:
raise SystemExit("Could not determine mapDamage version: %s" % (error,))
View it on GitLab: https://salsa.debian.org/med-team/mapdamage/-/commit/a9e55a3f0ab6028446b146672cdcd8fa1f045427
--
View it on GitLab: https://salsa.debian.org/med-team/mapdamage/-/commit/a9e55a3f0ab6028446b146672cdcd8fa1f045427
You're receiving this email because of your account on salsa.debian.org.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://alioth-lists.debian.net/pipermail/debian-med-commit/attachments/20200629/dbd3b18b/attachment-0001.html>
More information about the debian-med-commit
mailing list