[extract] Add bam flags

tanghaibao · Jul 6, 2018 · 8a8d3a5 · 8a8d3a5
1 parent e576250
commit 8a8d3a5
Showing 1 changed file with 14 additions and 1 deletion.
diff --git a/extract.go b/extract.go
@@ -589,8 +589,21 @@ func (r *Extracter) ExtractInterContigLinks() {
 			}
 			break
 		}
+		// Filtering: Unmapped | Secondary | QCFail | Duplicate | Supplementary
+		if rec.MapQ == 0 || rec.Flags&3844 != 0 {
+			continue
+		}
+
+		// Make sure we have these contig ids
 		at, bt := rec.Ref.Name(), rec.MateRef.Name()
-		ai, bi := r.contigToIdx[at], r.contigToIdx[bt]
+		ai, ok := r.contigToIdx[at]
+		if !ok {
+			continue
+		}
+		bi, ok := r.contigToIdx[bt]
+		if !ok {
+			continue
+		}
 
 		//         read1                                               read2
 		//     ---a-- X|----- dist = a2 ----|         |--- dist = b ---|X ------ b2 ------