- while((pt = getline(fin, line)) && line[0] != '>') {
- gseq += line;
- }
-
- size_t len = gseq.length();
- assert(len > 0);
- for (size_t j = 0; j < len; j++) gseq[j] = check(gseq[j]);
-
- iter = sn2tr.find(seqname);
- if (iter == sn2tr.end()) continue;
-
- chrvec.push_back(ChrInfo(seqname, len));
-
- vector<int>& vec = iter->second;
- int s = vec.size();
- for (int j = 0; j < s; j++) {
- assert(vec[j] > 0 && vec[j] <= M);
- transcripts.getTranscriptAt(vec[j]).extractSeq(gseq, seqs[vec[j]]);
- }
+ while((getline(fin, line)) && (line[0] != '>')) {
+ gseq += line;
+ }
+
+ size_t len = gseq.length();
+ assert(len > 0);
+ for (size_t j = 0; j < len; j++) gseq[j] = check(gseq[j]);
+
+ iter = sn2tr.find(seqname);
+ if (iter == sn2tr.end()) continue;
+
+ chrvec.push_back(ChrInfo(seqname, len));
+
+ vector<int>& vec = iter->second;
+ int s = vec.size();
+ for (int j = 0; j < s; j++) {
+ assert(vec[j] > 0 && vec[j] <= M);
+ transcripts.getTranscriptAt(vec[j]).extractSeq(gseq, seqs[vec[j]]);
+ }