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

fixing bug in lumpy parser

parent f2a499b0
No related branches found
No related tags found
No related merge requests found
...@@ -213,7 +213,7 @@ class SVReader(object): ...@@ -213,7 +213,7 @@ class SVReader(object):
if (sv[4] != "." if (sv[4] != "."
and sv[4] != "-1" and sv[4] != "-1"
and float(sv[-1]) < cutoff_proximity): and float(sv[-1]) < cutoff_proximity):
print("-".join(sv)) #print("-".join(sv))
excluded[sv[3]] = 1 excluded[sv[3]] = 1
vcf_records_filtered = [] vcf_records_filtered = []
......
...@@ -184,7 +184,7 @@ class LumpyReader(SVReader): ...@@ -184,7 +184,7 @@ class LumpyReader(SVReader):
# Filtering criteria more than 4 PE # Filtering criteria more than 4 PE
# see # see
# http://bcb.io/2014/08/12/validated-whole-genome-structural-variation-detection-using-multiple-callers # http://bcb.io/2014/08/12/validated-whole-genome-structural-variation-detection-using-multiple-callers
return record.MaxIndSupportingRP() > 4 and record.sr >0 return record.MaxIndSupportingRP() > 4 and record.sr[0] >0
class LumpyWriter(SVWriter): class LumpyWriter(SVWriter):
......
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