[med-svn] [Git][med-team/mapdamage][master] 13 commits: New upstream version 2.2.0+dfsg
Andreas Tille
gitlab at salsa.debian.org
Fri Dec 13 13:06:25 GMT 2019
Andreas Tille pushed to branch master at Debian Med / mapdamage
Commits:
2d75da74 by Andreas Tille at 2019-12-13T10:54:26Z
New upstream version 2.2.0+dfsg
- - - - -
81bd3d52 by Andreas Tille at 2019-12-13T10:54:26Z
routine-update: New upstream version
- - - - -
3e895fc8 by Andreas Tille at 2019-12-13T10:54:27Z
Update upstream source from tag 'upstream/2.2.0+dfsg'
Update to upstream version '2.2.0+dfsg'
with Debian dir 7101accc542f8ec0430dafc97b4c2f37229122db
- - - - -
ddcb5eca by Andreas Tille at 2019-12-13T10:57:08Z
Refresh patches
- - - - -
99e21e7e by Andreas Tille at 2019-12-13T11:53:34Z
R-U: Trailing whitespace in debian/changelog
- - - - -
b2ef1bb5 by Andreas Tille at 2019-12-13T11:53:34Z
R-U: autopkgtest: s/ADTTMP/AUTOPKGTEST_TMP/g
- - - - -
d9c80f43 by Andreas Tille at 2019-12-13T11:53:39Z
Set upstream metadata fields: Bug-Database.
- - - - -
9032df80 by Andreas Tille at 2019-12-13T12:09:46Z
Fix and enhance autopkgtest
- - - - -
406d885e by Andreas Tille at 2019-12-13T12:32:52Z
Fix and enhance build time and autopkgtest
- - - - -
28bfdf6b by Andreas Tille at 2019-12-13T12:38:18Z
More fixes in autopkgtest
- - - - -
8fd6fda8 by Andreas Tille at 2019-12-13T12:43:24Z
Respect DEB_BUILD_OPTIONS
- - - - -
6dd5c804 by Andreas Tille at 2019-12-13T13:01:25Z
Drop old tests since upstream is not providing according data any more
- - - - -
90380206 by Andreas Tille at 2019-12-13T13:04:27Z
Upload to unstable
- - - - -
15 changed files:
- README.md
- debian/changelog
- debian/control
- debian/docs
- debian/install
- debian/patches/use_debian_packaged_seqtk.patch
- debian/rules
- debian/tests/run-unit-test
- debian/upstream/metadata
- mapdamage/mp_test.py
- mapdamage/rescale.py
- + mapdamage/tests/probs/Stats_out_MCMC_correct_prob.csv
- + mapdamage/tests/test.rescaled.correct.sam
- mapdamage/version.py
- setup.py
Changes:
=====================================
README.md
=====================================
@@ -1,19 +1,35 @@
## mapDamage
-[](http://bioconda.github.io/recipes/mapdamage2/README.html) [](https://anaconda.org/bioconda/mapdamage2/files)
+[](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.0** **with** R and 4 mandatory packages for the Bayesian inference:
+```
+conda install -c bioconda mapdamage2=2.2.0=pyr36_1
+```
---
### Important
-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.
+
+- From version `2.2.0` 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.
### Introduction
Complete documentation, instructions, examples, screenshots and FAQ are available at [this address](http://ginolhac.github.io/mapDamage/).
-[mapDamage2](http://geogenetics.ku.dk/publications/mapdamage2.0/) is a computational framework written in **Python3** and **R**, which tracks and quantifies DNA damage patterns
+[mapDamage2](https://geogenetics.ku.dk/publications/mapdamage2.0/) is a computational framework written in **Python3** and **R**, which tracks and quantifies DNA damage patterns
among ancient DNA sequencing reads generated by Next-Generation Sequencing platforms.
-`mapDamage` was developed at the [Centre for GeoGenetics](http://geogenetics.ku.dk/) by the [Orlando Group ](http://geogenetics.ku.dk/research/research_groups/palaeomix_group/).
+`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
@@ -29,15 +45,30 @@ Ginolhac A, Rasmussen M, Gilbert MT, Willerslev E, Orlando L.
http://bioinformatics.oxfordjournals.org/content/27/15/2153](http://bioinformatics.oxfordjournals.org/content/27/15/2153)
-### Test
+### Test the no-stats part and rescaling
-in the package, you can test `mapDamage` by running:
+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/'
+pdf tests/results/Fragmisincorporation_plot.pdf generated
+additional tests/results/Length_plot.pdf generated
+Successful run
+.
+----------------------------------------------------------------------
+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:
=====================================
debian/changelog
=====================================
@@ -1,9 +1,19 @@
+mapdamage (2.2.0+dfsg-1) unstable; urgency=medium
+
+ * New upstream version
+ * Remove trailing whitespace in debian/changelog
+ * autopkgtest: s/ADTTMP/AUTOPKGTEST_TMP/g
+ * Set upstream metadata fields: Bug-Database.
+ * Fix and enhance build time and autopkgtest
+
+ -- Andreas Tille <tille at debian.org> Fri, 13 Dec 2019 14:01:31 +0100
+
mapdamage (2.1.1+dfsg-1) unstable; urgency=medium
* Team upload.
- * New upstream version
+ * New upstream version
- eliminates two patches
- + na.rm
+ + na.rm
+ import
- rescale_test.py no longer available - simpler build process
* Standards-Version: 4.4.1
=====================================
debian/control
=====================================
@@ -4,10 +4,15 @@ Uploaders: Andreas Tille <tille at debian.org>
Section: science
Priority: optional
Build-Depends: debhelper-compat (= 12),
- python3-all,
+ python3,
dh-python,
python3-setuptools,
- python3-pysam <!nocheck>
+ python3-pysam <!nocheck>,
+ seqtk <!nocheck>,
+ r-base-core <!nocheck>,
+ r-cran-gam <!nocheck>,
+ r-cran-ggplot2 <!nocheck>,
+ r-cran-inline <!nocheck>
Standards-Version: 4.4.1
Vcs-Browser: https://salsa.debian.org/med-team/mapdamage
Vcs-Git: https://salsa.debian.org/med-team/mapdamage.git
=====================================
debian/docs
=====================================
@@ -1,4 +1,3 @@
debian/README.test
debian/tests/run-unit-test
debian/tests/run-r-packages-check
-debian/tests/ref_T.fa
=====================================
debian/install
=====================================
@@ -1,2 +1 @@
-#No longer available?
-#mapdamage/rescale_test usr/share/doc/mapdamage/tests
+mapdamage/tests/probs usr/lib/python3/dist-packages/mapdamage/tests
=====================================
debian/patches/use_debian_packaged_seqtk.patch
=====================================
@@ -3,13 +3,11 @@ Last-Update: Thu, 28 Jul 2016 15:13:14 +0200
Bug-Debian: https://bugs.debian.org/859090
Description: Use Debian packaged seqtk
-Index: mapdamage/setup.py
-===================================================================
---- mapdamage.orig/setup.py
-+++ mapdamage/setup.py
-@@ -6,21 +6,6 @@ from distutils.command.install import in
- import os
- import subprocess
+--- a/setup.py
++++ b/setup.py
+@@ -11,21 +11,6 @@ if sys.version_info < (3, 5):
+ print("At least Python 3.5 is required.\n", file=sys.stderr)
+ exit(1)
-def compile_seqtk():
- """Compiling the seqtk toolkit"""
@@ -29,7 +27,7 @@ Index: mapdamage/setup.py
def setup_version():
if not os.path.exists(".git"):
# Release version, no .git folder
-@@ -40,17 +25,8 @@ class compileInstall(DistutilsInstall):
+@@ -45,17 +30,8 @@ class compileInstall(DistutilsInstall):
def run(self):
self.record=""
setup_version()
@@ -47,19 +45,17 @@ Index: mapdamage/setup.py
setup(
cmdclass={'install': compileInstall},
-@@ -59,7 +35,7 @@ setup(
- author='Aurélien Ginolhac, Mikkel Schubert, Hákon Jónsson',
+@@ -64,7 +40,7 @@ setup(
+ author='Aurélien Ginolhac, Mikkel Schubert, Ãkon Jónsson',
author_email='MSchubert at snm.ku.dk, jonsson.hakon at gmail.com',
packages=['mapdamage'],
-- package_data={'mapdamage': ['Rscripts/*.R','Rscripts/stats/*.R','seqtk/seqtk']},
-+ package_data={'mapdamage': ['Rscripts/*.R','Rscripts/stats/*.R']},
+- package_data={'mapdamage': ['Rscripts/*.R','Rscripts/stats/*.R','tests/*','seqtk/seqtk']},
++ package_data={'mapdamage': ['Rscripts/*.R','Rscripts/stats/*.R','tests/*']},
scripts=['bin/mapDamage'],
url='https://github.com/ginolhac/mapDamage',
license='LICENSE.txt',
-Index: mapdamage/bin/mapDamage
-===================================================================
---- mapdamage.orig/bin/mapDamage
-+++ mapdamage/bin/mapDamage
+--- a/bin/mapDamage
++++ b/bin/mapDamage
@@ -163,7 +163,7 @@ def main():
sys.path.insert(0,path_to_mm)
import mapdamage
@@ -69,10 +65,8 @@ Index: mapdamage/bin/mapDamage
if not (os.path.isfile(fpath_seqtk) and os.access(fpath_seqtk, os.X_OK)):
sys.stderr.write("Seqtk executable not accessible; mapDamage has not\n"
"been intalled properly or current user does not\n"
-Index: mapdamage/mapdamage/composition.py
-===================================================================
---- mapdamage.orig/mapdamage/composition.py
-+++ mapdamage/mapdamage/composition.py
+--- a/mapdamage/composition.py
++++ b/mapdamage/composition.py
@@ -33,7 +33,7 @@ def get_base_comp(filename,destination=F
Gets the basecomposition of all the sequences in filename
and returns the value to destination if given.
=====================================
debian/rules
=====================================
@@ -13,11 +13,13 @@ override_dh_fixperms:
find debian -name checkLibraries.R -exec chmod +x \{\} \;
find debian -name runGeneral.R -exec chmod +x \{\} \;
-# rescale_test.py seems to have left the building
-#override_dh_install:
-# dh_install
-# mv `find debian -name rescale_test.py` $(TESTDIR)
-# sed -i 's/^import rescale/from mapdamage &/' $(TESTDIR)/rescale_test.py
+override_dh_auto_test:
+ifeq (,$(filter nocheck,$(DEB_BUILD_OPTIONS)))
+ cd mapdamage && \
+ PATH=$(PATH):$(CURDIR)/bin \
+ PYTHONPATH=$(shell pybuild --print build_dir --interpreter python3) \
+ python3 mp_test.py
+endif
override_dh_python3:
dh_python3 --no-ext-rename
=====================================
debian/tests/run-unit-test
=====================================
@@ -2,20 +2,23 @@
pkg=mapdamage
-if [ "$ADTTMP" = "" ] ; then
- ADTTMP=$(mktemp -d /tmp/${pkg}-test.XXXXXX)
- trap "rm -rf $ADTTMP" 0 INT QUIT ABRT PIPE TERM
+if [ "$AUTOPKGTEST_TMP" = "" ] ; then
+ AUTOPKGTEST_TMP=$(mktemp -d /tmp/${pkg}-test.XXXXXX)
+ trap "rm -rf $AUTOPKGTEST_TMP" 0 INT QUIT ABRT PIPE TERM
fi
-cd $ADTTMP
+cd $AUTOPKGTEST_TMP
+
+cp -a /usr/lib/python3/dist-packages/mapdamage/tests .
+cp /usr/lib/python3/dist-packages/mapdamage/mp_test.py .
+python3 mp_test.py
+
+exit 0
-cp -a /usr/share/doc/${pkg}/tests/* .
cp -a /usr/share/doc/${pkg}/ref_T.fa .
gunzip -r *
-# FIXME: rescale_test.py throws errors
-# will discuss with upstream
-# python3 rescale_test.py
+# these old test so not work any more since rescale_test/ was removed upstream
mapDamage -i rescale_test/pe_test/pe.sam -r ref_T.fa
pdf_result="results_pe/Fragmisincorporation_plot.pdf"
=====================================
debian/upstream/metadata
=====================================
@@ -1,27 +1,28 @@
Reference:
- - Author: >
+- Author: >
Hákon Jónsson and Aurélien Ginolhac and Mikkel Schubert and
Philip Johnson and Ludovic Orlando
- Title: >
+ Title: >
mapDamage2.0: fast approximate Bayesian estimates of ancient DNA
damage parameters
- Journal: Bioinformatics
- Year: 2013
- Volume: 29
- Number: 13
- Pages: 1682-4
- DOI: 10.1093/bioinformatics/btt193
- PMID: 23613487
- URL: >
+ Journal: Bioinformatics
+ Year: 2013
+ Volume: 29
+ Number: 13
+ Pages: 1682-4
+ DOI: 10.1093/bioinformatics/btt193
+ PMID: 23613487
+ URL: >
http://bioinformatics.oxfordjournals.org/content/early/2013/04/23/bioinformatics.btt193.abstract
- eprint: >
+ eprint: >
http://bioinformatics.oxfordjournals.org/content/early/2013/04/23/bioinformatics.btt193.full.pdf+html
Registry:
- - Name: SciCrunch
- Entry: SCR_001240
- - Name: bio.tools
- Entry: NA
- - Name: OMICtools
- Entry: OMICS_02099
- - Name: conda:bioconda
- Entry: mapdamage2
+- Name: SciCrunch
+ Entry: SCR_001240
+- Name: bio.tools
+ Entry: NA
+- Name: OMICtools
+ Entry: OMICS_02099
+- Name: conda:bioconda
+ Entry: mapdamage2
+Bug-Database: https://github.com/ginolhac/mapDamage/issues
=====================================
mapdamage/mp_test.py
=====================================
@@ -1,12 +1,28 @@
import unittest
import subprocess
-
+import pysam
+import mapdamage
+import optparse
+import filecmp
def read_nocomment(dnacomp):
with open(dnacomp, 'r') as f:
a = f.readlines()
return [x for x in a if not x.startswith('#')]
+def mock_options(filename,rescale_out,folder):
+ """Make the options object with nice values for testing"""
+ return optparse.Values({
+ "filename":filename,
+ "rescale_out":rescale_out,
+ "verbose":True,
+ "folder":folder,
+ "rescale_length_5p":12, # default values as in --seq-length
+ "rescale_length_3p":12, # default values as in --seq-length
+ "quiet":True
+ })
+
+
class testCases(unittest.TestCase):
def test_no_stats(self):
@@ -14,5 +30,31 @@ class testCases(unittest.TestCase):
subprocess.run(["mapDamage", "-i", "tests/test.bam", "-r", "tests/fake1.fasta", "-d", "tests/results", "--no-stats"], check = True)
self.assertTrue(read_nocomment("tests/dnacomp.txt") == read_nocomment("tests/results/dnacomp.txt"))
+class testRescaling(unittest.TestCase):
+ def test_single_end_file(self):
+ """Test, rescaling BAM file"""
+ #
+ # The expected substition frequencies before and after scaling using the scaled qualities as probalities:
+ # CT 0.06226411977920493 0.04163524443356556
+ # TC 0.020395286584806528 0.020395286584806528
+ # GA 0.04400459948304954 0.03794905109091021
+ # AG 0.05355350777642983 0.05355350777642983
+ # Quality metrics before and after scaling
+ # CT-Q0 5 5
+ # CT-Q10 5 2
+ # CT-Q20 3 2
+ # CT-Q30 3 2
+ # CT-Q40 0 0
+ # GA-Q0 5 5
+ # GA-Q10 5 4
+ # GA-Q20 1 0
+ # GA-Q30 1 0
+ # GA-Q40 0 0
+ options = mock_options("tests/test.bam","tests/test.rescaled.sam","tests/probs/")
+ ref = pysam.Fastafile("tests/fake1.fasta")
+ mapdamage.rescale.rescale_qual(ref,options,debug=True)
+ self.assertTrue(filecmp.cmp("tests/test.rescaled.sam","tests/test.rescaled.correct.sam"))
+
+
if __name__=='__main__':
unittest.main()
=====================================
mapdamage/rescale.py
=====================================
@@ -19,8 +19,8 @@ def phred_char_to_pval(ch):
def get_corr_prob(folder, rescale_length_5p, rescale_length_3p):
"""
- Reads the damage probability correction file, returns a
- dictionary with this structure
+ Reads the damage probability correction file, returns a
+ dictionary with this structure
position (one based) - CT - probability
- GA - probability
"""
@@ -28,21 +28,21 @@ def get_corr_prob(folder, rescale_length_5p, rescale_length_3p):
if not os.path.isfile(full_path):
sys.exit("Missing file, the file \n\tStats_out_MCMC_correct_prob.csv\nshould be in the folder\n\t"+folder+"\nDid you run the MCMC estimates of the parameters?")
try:
- fi_handle = csv.DictReader(open(full_path))
- corr_prob = {}
- for line in fi_handle:
- if (line["Position"] in corr_prob):
- sys.exit('This file has multiple position definitions %s, line %d: %s' % \
- (folder, fi_handle.line_num, corr_prob[line["Position"]]))
- else:
- corr_prob[int(line["Position"])] = {'C.T':float(line["C.T"]), 'G.A':float(line["G.A"])}
-
- # Exclude probabilities for positions outside of user-specified region
- for key in list(corr_prob.keys()):
- if key < -rescale_length_3p or key > rescale_length_5p:
- corr_prob.pop(key)
+ with open(full_path) as fi:
+ fi_handle = csv.DictReader(fi)
+ corr_prob = {}
+ for line in fi_handle:
+ if (line["Position"] in corr_prob):
+ sys.exit('This file has multiple position definitions %s, line %d: %s' % \
+ (folder, fi_handle.line_num, corr_prob[line["Position"]]))
+ else:
+ corr_prob[int(line["Position"])] = {'C.T':float(line["C.T"]), 'G.A':float(line["G.A"])}
- return corr_prob
+ # Exclude probabilities for positions outside of user-specified region
+ for key in list(corr_prob.keys()):
+ if key < -rescale_length_3p or key > rescale_length_5p:
+ corr_prob.pop(key)
+ return corr_prob
except csv.Error as e:
sys.exit('File %s, line %d: %s' % (os.path.join(folder,"Stats_out_MCMC_correct_prob.csv"), \
fi_handle.line_num, e))
@@ -50,11 +50,11 @@ def get_corr_prob(folder, rescale_length_5p, rescale_length_3p):
def corr_this_base(corr_prob, nt_seq, nt_ref, pos, length,direction="both"):
"""
- The position specific damaging correction, using the input
- corr_prob dictionary holding the damage correcting values
- nt_seq nucleotide in the sequence
+ The position specific damaging correction, using the input
+ corr_prob dictionary holding the damage correcting values
+ nt_seq nucleotide in the sequence
nt_ref nucleotide in the reference
- pos relative position from the 5' end
+ pos relative position from the 5' end
length length of the sequence
direction which end to consider the rescaling
returns the correction probability for this particular set
@@ -71,8 +71,8 @@ def corr_this_base(corr_prob, nt_seq, nt_ref, pos, length,direction="both"):
else:
# other transitions/transversions are not affected by damage
return 0
-
- back_pos = pos-length-1
+
+ back_pos = pos-length-1
# position from 3' end
if pos in corr_prob:
@@ -90,7 +90,7 @@ def corr_this_base(corr_prob, nt_seq, nt_ref, pos, length,direction="both"):
if direction == "forward":
return p5_corr
elif direction == "backward":
- return p3_corr
+ return p3_corr
elif direction == "both":
if pos < abs(back_pos) :
# then we use the forward correction
@@ -151,7 +151,7 @@ def record_subs(subs,nt_seq,nt_ref,nt_qual,nt_newqual,prob_corr):
else:
sub_type = "NN"
if (sub_type != "NN"):
- # record only transitions
+ # record only transitions
subs[sub_type+"-before"][int(ord(nt_qual))-33] += 1
subs[sub_type+"-after"][int(ord(nt_newqual))-33] += 1
if (nt_ref in ["A","C","G","T"]):
@@ -175,7 +175,7 @@ def print_subs(subs):
if subs["C"]!=0:
# the special case of no substitutions
print(("\tCT\t"+str(subs["CT-pvals_before"]/subs["C"])+"\t\t"+str(subs["CT-pvals"]/subs["C"])))
- else:
+ else:
print("\tCT\tNA\t\tNA")
if subs["T"]!=0:
print(("\tTC\t"+str(subs["TC-pvals"]/subs["T"])+"\t\t"+str(subs["TC-pvals"]/subs["T"])))
@@ -200,19 +200,19 @@ def print_subs(subs):
print(("\tGA-Q20 \t"+str(subs["GA-before-Q20"])+"\t\t"+str(subs["GA-after-Q20"])))
print(("\tGA-Q30 \t"+str(subs["GA-before-Q30"])+"\t\t"+str(subs["GA-after-Q30"])))
print(("\tGA-Q40 \t"+str(subs["GA-before-Q40"])+"\t\t"+str(subs["GA-after-Q40"])))
-
+
def rescale_qual_read(bam, read, ref, corr_prob,subs, debug = False,direction="both"):
"""
bam a pysam bam object
read a pysam read object
ref a pysam fasta ref file
- reflengths a dictionary holding the length of the references
- subs a dictionary holding the corrected number of substition before and after scaling
+ reflengths a dictionary holding the length of the references
+ subs a dictionary holding the corrected number of substition before and after scaling
corr_prob dictionary from get_corr_prob
returns a read with rescaled quality score
-
- Iterates through the read and reference, rescales the quality
+
+ Iterates through the read and reference, rescales the quality
according to corr_prob
"""
if not debug:
@@ -238,7 +238,7 @@ def rescale_qual_read(bam, read, ref, corr_prob,subs, debug = False,direction="b
pos_on_read = 0
number_of_rescaled_bases = 0.0
for (i, nt_seq, nt_ref, nt_qual) in zip(range(length_align), seq, refseq, qual):
- # rescale the quality according to the triplet position,
+ # rescale the quality according to the triplet position,
# pair of the reference and the sequence
if ((nt_seq == "T" and nt_ref =="C") or (nt_seq == "A" and nt_ref =="G")):
# need to rescale this subs.
@@ -250,13 +250,13 @@ def rescale_qual_read(bam, read, ref, corr_prob,subs, debug = False,direction="b
else:
# don't rescale, other bases
newp = 1 - phred_char_to_pval(nt_qual)
- newq = nt_qual
+ newq = nt_qual
if pos_on_read < length_read:
- new_qual[pos_on_read] = newq
+ new_qual[pos_on_read] = newq
record_subs(subs,nt_seq,nt_ref,nt_qual,new_qual[pos_on_read],newp)
if nt_seq != "-":
pos_on_read += 1
- # done with the aligned portion of the read
+ # done with the aligned portion of the read
else:
if not debug:
logger.warning("Warning: The aligment of the read is longer than the actual read %s",(read.qname))
@@ -266,7 +266,7 @@ def rescale_qual_read(bam, read, ref, corr_prob,subs, debug = False,direction="b
if read.is_reverse:
new_qual = new_qual[::-1]
if (read.cigar[0][0] == 4):
- # check for soft clipping at forward end
+ # check for soft clipping at forward end
new_qual = read.qual[0:read.cigar[0][1]] + new_qual
if (read.cigar[-1][0] == 4):
# the same backwards
@@ -285,7 +285,7 @@ def rescale_qual_read(bam, read, ref, corr_prob,subs, debug = False,direction="b
def rescale_qual(ref, options,debug=False):
- """
+ """
ref a pysam fasta ref file
bam_filename name of a BAM/SAM file to read
fi file containing the csv with correction probabilities
@@ -319,20 +319,20 @@ def rescale_qual(ref, options,debug=False):
pass
elif not hit.qual and not debug:
logger.warning("Cannot rescale base PHRED scores for read '%s'; no scores assigned." % hit.qname)
- elif hit.is_paired :
+ elif hit.is_paired :
if first_pair and not debug:
- # assuming the ends are non-overlapping
+ # assuming the ends are non-overlapping
logger.warning("Warning! Assuming the pairs are non-overlapping, facing inwards and correctly paired.")
first_pair=False
#5p --------------> 3p
#3p <-------------- 5p
# pair 1 (inwards)
- #5p ---->
+ #5p ---->
# <---- 5p
# A B
- # pair 2 (outwards), this happens if the reference is RC this is not supported
+ # pair 2 (outwards), this happens if the reference is RC this is not supported
# ----> 3p
- #3p <----
+ #3p <----
# A B
# Correct outwards pairs from the 3p and inwards pairs with the 5p end
if ((not hit.is_reverse) and hit.mate_is_reverse and (hit.pnext>hit.pos) and hit.tid==hit.mrnm):
=====================================
mapdamage/tests/probs/Stats_out_MCMC_correct_prob.csv
=====================================
@@ -0,0 +1,25 @@
+"","Position","C.T","G.A"
+"1",1,0.653616777018436,0
+"2",2,0.574821404505089,0
+"3",3,0.524660717195785,0
+"4",4,0.491676371805506,0
+"5",5,0.46919312975022,0
+"6",6,0.453361082263795,0
+"7",7,0.441882878430246,0
+"8",8,0.433345603590333,0
+"9",9,0.426853067214069,0
+"10",10,0.4218187429273,0
+"11",11,0.417847465539491,0
+"12",12,0.414666179579865,0
+"13",-12,0,0.427459648772127
+"14",-11,0,0.430678304009036
+"15",-10,0,0.434694598982698
+"16",-9,0,0.439783285536968
+"17",-8,0,0.446341431282045
+"18",-7,0,0.454957458647015
+"19",-6,0,0.466528745659899
+"20",-5,0,0.482466321045583
+"21",-4,0,0.505055143009322
+"22",-3,0,0.538101706832429
+"23",-2,0,0.588153916474448
+"24",-1,0,0.666277684522864
=====================================
mapdamage/tests/test.rescaled.correct.sam
=====================================
@@ -0,0 +1,11 @@
+ at SQ SN:fake1 LN:201
+ at PG ID:bwa PN:bwa VN:0.7.17-r1188 CL:bwa samse -f input.sam fake1.fasta - input.fastq
+fake1:50-100 0 fake1 69 37 31M * 0 0 CCATGTCGGGCAGGCTGGTCTCGAACTCCTG ////////E6/EE/E/A////E//AAAE/EE XT:A:U NM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:31 MR:f:0
+fake1:50-100b 0 fake1 69 37 31M * 0 0 CTATGTCGGGCAGGCTGGTCTCGAACTCCTG /#//////E6/EE/E/A////E//AAAE/EE XT:A:U NM:i:1 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:1C29 MR:f:0.57482
+fake1:50-100c 4 * 0 0 * * 0 0 CGGTAGAGATGGAGTTTCACCATGTCGGGCAAG //////////////A////////////E6/EE/
+fake1:50-100d 4 * 0 0 * * 0 0 TAATAGAGATGGAGTTTCACCATGTCGTGCAAG //////////////A////////////E6/EE/
+fake1:70-120 0 fake1 70 37 38M * 0 0 TATGTCGGGCAGGCTGGTCTCGAACTCCTGACCTCAGA #////6EE//AA6E//E66E/EAEEEEEAEEEAAEEA# XT:A:U NM:i:2 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:0C36G0 MR:f:1.31989
+fake1:80-130 16 fake1 80 25 41M * 0 0 GAGCTGGTCTCGAACTCCTGACCTCAGGCGATCTGCCTGTC AAABBCCDDFEEEEEAEEEEA/EE/E/AAE//A//E/A/// XT:A:U NM:i:3 X0:i:1 X1:i:0 XM:i:3 XO:i:0 XG:i:0 MD:Z:0A0G37C1 MR:f:0
+fake1:80-130b 16 fake1 80 25 51M * 0 0 AGGCCGGTCTCGAACTCCTGACCTCAGGCGATCTGCCTGCCTTAACCTCCC AAABBCCDDFEEEEEAEEEEA/EE/E/AAE//A//E/A///E/A$/EE66/ XT:A:U NM:i:3 X0:i:1 X1:i:0 XM:i:3 XO:i:0 XG:i:0 MD:Z:4T37C1G6 MR:f:0.44188
+fake1:80-130c 16 fake1 80 37 51M * 0 0 AGGCTGGTCTCGAACTCCTGACCTCAGGCGATCTGCCTGCCTTAGCCCCCC AAABBCCDDFEEEEEAEEEEA/EE/E/AAE//A//E/A///E/A//EE66/ XT:A:U NM:i:2 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:42C4T3 MR:f:0
+fake1:80-130d 16 fake1 80 25 51M * 0 0 AGGCTGGTCTCGAACTCCTGACCTCAGACGATCTGCCTGCCTTAGCCCCCC AAABBCCDDFEEEEEAEEEEA/EE/E/AAE//A//E/A///E/A//EE66/ XT:A:U NM:i:3 X0:i:1 X1:i:0 XM:i:3 XO:i:0 XG:i:0 MD:Z:27G14C4T3 MR:f:0
=====================================
mapdamage/version.py
=====================================
@@ -2,4 +2,4 @@
try:
from mapdamage._version import __version__
except ImportError:
- __version__ = "2.1.0"
+ __version__ = "2.2.0"
=====================================
setup.py
=====================================
@@ -4,8 +4,13 @@
from distutils.core import setup
from distutils.command.install import install as DistutilsInstall
import os
+import sys
import subprocess
+if sys.version_info < (3, 5):
+ print("At least Python 3.5 is required.\n", file=sys.stderr)
+ exit(1)
+
def compile_seqtk():
"""Compiling the seqtk toolkit"""
old_wd = os.getcwd()
@@ -55,11 +60,11 @@ class compileInstall(DistutilsInstall):
setup(
cmdclass={'install': compileInstall},
name='mapdamage',
- version='2.1.0',
- author='Aurélien Ginolhac, Mikkel Schubert, Hákon Jónsson',
+ version='2.2.0',
+ author='Aurélien Ginolhac, Mikkel Schubert, Ãkon Jónsson',
author_email='MSchubert at snm.ku.dk, jonsson.hakon at gmail.com',
packages=['mapdamage'],
- package_data={'mapdamage': ['Rscripts/*.R','Rscripts/stats/*.R','seqtk/seqtk']},
+ package_data={'mapdamage': ['Rscripts/*.R','Rscripts/stats/*.R','tests/*','seqtk/seqtk']},
scripts=['bin/mapDamage'],
url='https://github.com/ginolhac/mapDamage',
license='LICENSE.txt',
View it on GitLab: https://salsa.debian.org/med-team/mapdamage/compare/ac991c66c1eb17abed708778741089fed0ddd3f5...90380206100e9dd6c74af3b0deae952371ccbcbd
--
View it on GitLab: https://salsa.debian.org/med-team/mapdamage/compare/ac991c66c1eb17abed708778741089fed0ddd3f5...90380206100e9dd6c74af3b0deae952371ccbcbd
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/20191213/3c243d23/attachment-0001.html>
More information about the debian-med-commit
mailing list