We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent e623bc5 commit 73c38a0Copy full SHA for 73c38a0
1 file changed
TCSfilter.py
@@ -35,7 +35,7 @@ def ListParser(infile_name, minqual, mincov):
35
covs = sorted(list(map(int, covs)))
36
if tcov <= mincov: # Discard positions with less than mincov
37
pass
38
- elif int(lines.split()[12]) >= tcov//2: # Discard pos with many gaps
+ elif int(lines.split()[12]) >= tcov//2 or lines.split()[2] == "-": # Discard pos with many gaps
39
40
elif covs[-2] <= (ceil(tcov * 0.2)): # Discard low freq second variant
41
0 commit comments