| 52 | map<int, READ_INT_TYPE>::iterator iter; |
| 53 | |
| 54 | void init(const char* imdName, const char* alignF) { |
| 55 | parser = new SamParser(alignF, aux, transcripts, imdName); |
| 56 | |
| 57 | memset(cat, 0, sizeof(cat)); |
| 58 | memset(readOutFs, 0, sizeof(readOutFs)); |
| 59 | |
| 60 | int tmp_n_os = -1; |
| 61 | |
| 62 | for (int i = 0; i < 3; i++) { |
| 63 | genReadFileNames(imdName, i, read_type, n_os, readOutFs[i]); |
| 64 | |
| 65 | assert(tmp_n_os < 0 || tmp_n_os == n_os); tmp_n_os = n_os; |
| 66 | |
| 67 | for (int j = 0; j < n_os; j++) |
| 68 | cat[i][j] = new ofstream(readOutFs[i][j]); |
| 69 | } |
| 70 | |
| 71 | counter.clear(); |
| 72 | } |
| 73 | |
| 74 | //Do not allow duplicate for unalignable reads and supressed reads in SAM input |
| 75 | template<class ReadType, class HitType> |
no test coverage detected