Merge branch 'add_pair_read_name'

arq5x · Apr 16, 2019 · 8f93887 · 8f93887
2 parents ca61d4e + 95ffc02
commit 8f93887
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 0 deletions.
diff --git a/src/lumpy/SV_Pair.cpp b/src/lumpy/SV_Pair.cpp
@@ -52,6 +52,8 @@ SV_Pair(const BamAlignment &bam_a,
     else
         min_mapping_quality = bam_b.MapQuality;
 
+    read_id = bam_a.Name;
+
     struct interval tmp_a, tmp_b;
     tmp_a.start = bam_a.Position;
     tmp_a.end = bam_a.GetEndPosition(false, false) - 1;
@@ -341,6 +343,7 @@ void
 SV_Pair::
 print_evidence()
 {
+    cout << read_id << "\t";
     print_bedpe(0);
 }
 //}}}

diff --git a/src/lumpy/SV_Pair.h b/src/lumpy/SV_Pair.h
@@ -61,6 +61,7 @@ class SV_Pair: public SV_Evidence
 		struct interval read_r;
 		bool read_l_is_split, read_r_is_split;
 		SV_PairReader *reader;
+		std::string read_id;
 
 
 		SV_Pair(const BamAlignment &bam_a,