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

corrected pindel filter

parent b689588e
No related branches found
No related tags found
No related merge requests found
...@@ -377,7 +377,7 @@ class PindelReader(SVReader): ...@@ -377,7 +377,7 @@ class PindelReader(SVReader):
# return (record.length() > 60) # return (record.length() > 60)
def SpecificFilterPass(self, record): def SpecificFilterPass(self, record):
# fILTER # fILTER
if (abs(record.start-record.end+1) >= 2000 and if (abs(record.start-record.end+1) >= 2000 or
record.MaxIndSupportingRP() <= 3): record.MaxIndSupportingRP() <= 3):
return False return False
else: else:
......
...@@ -80,6 +80,9 @@ class VCFReader(SVReader): ...@@ -80,6 +80,9 @@ class VCFReader(SVReader):
def getHeader(self): def getHeader(self):
return self.vcf_reader.header return self.vcf_reader.header
def removeInfoHeaderLine(self, key):
self.vcf_reader.header.info.remove_header(key)
def addInfo(self, name, number, type, description): def addInfo(self, name, number, type, description):
self.vcf_reader.header.info.add(id=name, self.vcf_reader.header.info.add(id=name,
number=number, number=number,
......
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