fix bug in viral profiling, causin dublicated reads in viral file

smangul1 · Mar 15, 2018 · d4bbd4e · d4bbd4e
1 parent 2ea0b6a
commit d4bbd4e
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/rop_functions.py b/rop_functions.py
@@ -228,7 +228,7 @@ def nMicrobialReads(inFile_name, outFile_name,flag):
 
         if  alignmentLength >= 0.8*readLength and identity >= 0.9:
             readsMicrobiome.add(read.query_name)
-        set2file(readsMicrobiome,outFile_name)
+    set2file(readsMicrobiome,outFile_name)
 
     samfile.close()