Fixed a minor bug which only affects paired-end reads for reporting how many alignmen...

[rsem.git] / Transcript.h

diff --git a/Transcript.h b/Transcript.h
index 8d7be7ae68afae09c773144e45db4ba79b81b36f..90f6e6797b17f949b38e4776fcff7eacb8ace2ce 100644 (file)
--- a/Transcript.h
+++ b/Transcript.h
@@ -113,7 +113,10 @@ void Transcript::read(std::ifstream& fin) {
        int s;
        std::string tmp;
 
-       fin>>transcript_id>>gene_id>>seqname>>tmp>>length;
+       getline(fin, transcript_id);
+       getline(fin, gene_id);
+       getline(fin, seqname);
+       fin>>tmp>>length;
        assert(tmp.length() == 1 && (tmp[0] == '+' || tmp[0] == '-'));
        strand = tmp[0];
        structure.clear();
@@ -130,7 +133,10 @@ void Transcript::read(std::ifstream& fin) {
 void Transcript::write(std::ofstream& fout) {
        int s = structure.size();
 
-       fout<<transcript_id<<" "<<gene_id<<" "<<seqname<<" "<<strand<<" "<<length<<" ";
+       fout<<transcript_id<<std::endl;
+       fout<<gene_id<<std::endl;
+       fout<<seqname<<std::endl;
+       fout<<strand<<" "<<length<<std::endl;
        fout<<s;
        for (int i = 0; i < s; i++) fout<<" "<<structure[i].start<<" "<<structure[i].end;
        fout<<std::endl;