We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent ec3661c commit d1e5bd5Copy full SHA for d1e5bd5
src/junction.h
@@ -506,12 +506,14 @@ namespace torali
506
507
TSvtSRBamRecord altSR(2 * DELLY_SVT_TRANS, TSRBamRecord());
508
TConfig altConfig(c);
509
+ altConfig.hasExcludeFile = false;
510
+ altConfig.exclude.clear();
511
altConfig.genome = genome[i];
512
altConfig.files.clear();
513
altConfig.files.push_back(align[i]);
514
if (isBamCram(align[i].string())) {
515
// Alternate alignment in BAM/CRAM format
- samFile* samfile = sam_open(c.files[0].string().c_str(), "r");
516
+ samFile* samfile = sam_open(altConfig.files[0].string().c_str(), "r");
517
bam_hdr_t* hdr = sam_hdr_read(samfile);
518
TValidRegions vR;
519
_parseExcludeIntervals(altConfig, hdr, vR);
0 commit comments