We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent a5219e2 commit 983e83eCopy full SHA for 983e83e
lib/BtIO.py
@@ -107,7 +107,7 @@ def readBam(infile, set_of_blobs):
107
cigar_match_re = re.compile(r"(\d+)M") # only gets digits before M's
108
# execute samtools to get only mapped reads
109
#command = "samtools view -F 12 -F 256" + infile
110
- command = "samtools view -F 4" + infile
+ command = "samtools view -F 4 -F 256 " + infile
111
# ADD flag picard -F 1028 to not consider optical duplicates
112
#command = "samtools view -F 1028 " + infile
113
# only one counter since only yields mapped reads
0 commit comments