Skip to content

Commit

Permalink
Added a test for specifying the reference via samtools -T. This
Browse files Browse the repository at this point in the history
appeared to work in our previous tests, but did not due to the code
sneakily reverting to the @sq UR: field instead.
  • Loading branch information
jkbonfield authored and jmarshall committed Oct 30, 2014
1 parent 01465e8 commit 6c17071
Showing 1 changed file with 16 additions and 0 deletions.
16 changes: 16 additions & 0 deletions test/test.pl
Original file line number Diff line number Diff line change
Expand Up @@ -1825,6 +1825,22 @@ sub test_view
}
}

# Test CRAM with explicit -T
my $cram_no_ur = sprintf("%s.test%03d.cram", $out, $test);
run_view_test($opts,
msg => "$test: Make CRAM with no UR field",
args => ['-C', $sam_no_ur],
ref_path => "$$opts{path}/dat/cram_md5",
out => $cram_no_ur);
$test++;

run_view_test($opts,
msg => "$test: Decoding CRAM with no UR field via -T",
args => ['-T', $ref_file, $cram_no_ur],
out => sprintf("%s.test%03d.sam", $out, $test),
compare => $sam_no_ur);
$test++;


# CIGAR B-operator removal tests.
my $b_op_sam = "$$opts{path}/dat/view.003.sam";
Expand Down

0 comments on commit 6c17071

Please sign in to comment.