Added error detection for cases such as a read's two mates having different names...

[rsem.git] / Transcript.h

diff --git a/Transcript.h b/Transcript.h
index 90f6e6797b17f949b38e4776fcff7eacb8ace2ce..f7afa2b20542ebdb2f4883d1cc12631eb2b1176e 100644 (file)
--- a/Transcript.h
+++ b/Transcript.h
@@ -10,6 +10,10 @@
 
 #include "utils.h"
 
+/**
+   If no genome is provided, seqname field is used to store the allele name.
+ */
+
 struct Interval {
        int start, end;
 
@@ -49,7 +53,7 @@ public:
        }
 
        bool operator< (const Transcript& o) const {
-         return gene_id < o.gene_id || (gene_id == o.gene_id && transcript_id < o.transcript_id);
+         return gene_id < o.gene_id || (gene_id == o.gene_id && transcript_id < o.transcript_id) || (gene_id == o.gene_id && transcript_id == o.transcript_id && seqname < o.seqname);
        }
 
        const std::string& getTranscriptID() const { return transcript_id; }