Skip to content
Snippets Groups Projects
Commit 6f2a249e authored by Thomas Faraut's avatar Thomas Faraut
Browse files

Merge branch 'newfilter' of forgemia.inra.fr:svdetection/svlib into newfilter

parents 85d17402 d26a53fe
No related branches found
No related tags found
No related merge requests found
...@@ -130,8 +130,8 @@ class AnnotatedRecordSample(object): ...@@ -130,8 +130,8 @@ class AnnotatedRecordSample(object):
class VCFReader(SVReader): class VCFReader(SVReader):
def __init__(self, file_name, sv_to_report=None): def __init__(self, file_name, sv_to_report=None):
<<<<<<< HEAD
super(VCFReader, self).__init__(file_name) super(VCFReader, self).__init__(file_name)
self.vcf_reader = VariantFile(file_name)
self.filename = file_name self.filename = file_name
self.sv_to_report = sv_to_report self.sv_to_report = sv_to_report
self.vcf_reader = VariantFile(file_name) self.vcf_reader = VariantFile(file_name)
...@@ -172,7 +172,6 @@ class VCFReader(SVReader): ...@@ -172,7 +172,6 @@ class VCFReader(SVReader):
self.addFilter("LOWSUPPORT", self.addFilter("LOWSUPPORT",
"total supp reads < 5 or supp samples < 2") "total supp reads < 5 or supp samples < 2")
def getOrderedSamples(self): def getOrderedSamples(self):
samples = self.vcf_reader.header.samples samples = self.vcf_reader.header.samples
sample_names = [sample.rsplit('.')[0] for sample in samples] sample_names = [sample.rsplit('.')[0] for sample in samples]
......
...@@ -396,6 +396,7 @@ class PindelReader(SVReader): ...@@ -396,6 +396,7 @@ class PindelReader(SVReader):
else: else:
return True return True
def remove_duplicate(self, records): def remove_duplicate(self, records):
""" """
returns a vector of records where duplicates were removed returns a vector of records where duplicates were removed
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment