Merge pull request #100 from J35P312/master

TIDDIT 3.5.0
SciLifeLab · Jan 2, 2023 · 46f81fa · 46f81fa
2 parents 2bd9492 + 8a10667
commit 46f81fa
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 3 deletions.
diff --git a/setup.py b/setup.py
@@ -20,7 +20,7 @@
 
 setup(
     name = 'tiddit',
-    version = '3.4.0',
+    version = '3.5.0',
 
     url = "https://github.com/SciLifeLab/TIDDIT",
     author = "Jesper Eisfeldt",

diff --git a/tiddit/__main__.py b/tiddit/__main__.py
@@ -17,7 +17,7 @@
 import tiddit.tiddit_contig_analysis as tiddit_contig_analysis
 
 def main():
-	version="3.4.0"
+	version="3.5.0"
 	parser = argparse.ArgumentParser("""tiddit-{}""".format(version),add_help=False)
 	parser.add_argument("--sv"	 , help="call structural variation", required=False, action="store_true")
 	parser.add_argument("--cov"        , help="generate a coverage bed file", required=False, action="store_true")

diff --git a/tiddit/tiddit_coverage_analysis.pyx b/tiddit/tiddit_coverage_analysis.pyx
@@ -1,4 +1,5 @@
 import numpy
+import re
 cimport numpy
 import pysam
 import gzip
@@ -19,7 +20,9 @@ def get_gc(str reference_fasta,bam_header,bin_size):
 		with gzip.open(reference_fasta, 'r') as f:
 			sequence=f.read()
 
-	split_reference=sequence.split(">")
+
+	#split_reference=sequence.split(">")
+	split_reference=re.split("\n>|^>", sequence)
 	del sequence
 	del split_reference[0]