Skip to content

Commit

Permalink
Merge tag '2.0.5' into dev
Browse files Browse the repository at this point in the history
Bugfix - output filename for RSeQC module GeneCoverage has changed - updating process_qcstats gene_coverage subroutine to reflect this
  • Loading branch information
am26 committed Feb 23, 2016
2 parents ffe341f + 7499ab4 commit 502a916
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions perl/lib/Sanger/CGP/RnaQC/Implement.pm
Original file line number Diff line number Diff line change
Expand Up @@ -52,7 +52,7 @@ sub check_input {
# Check all the required input files exist
PCAP::Cli::file_for_reading('bam-stats', File::Spec->catfile($inputdir, "$sample.bas"));
PCAP::Cli::file_for_reading('bam-stats', File::Spec->catfile($inputdir, "$sample.transcriptome.bas"));
PCAP::Cli::file_for_reading('gene-coverage', File::Spec->catfile($inputdir, "$sample.coverage.geneBodyCoverage.r"));
PCAP::Cli::file_for_reading('gene-coverage', File::Spec->catfile($inputdir, "$sample.geneBodyCoverage.r"));
PCAP::Cli::file_for_reading('rrna-percentage', File::Spec->catfile($inputdir, "$sample.rrna.txt"));
PCAP::Cli::file_for_reading('read-distribution', File::Spec->catfile($inputdir, "$sample.read_dist.txt"));

Expand All @@ -64,8 +64,8 @@ sub gene_coverage {

my $inputdir = $options->{'indir'};
my $sample = $options->{'sample'};
my $gene_coverage_r_file = File::Spec->catfile($inputdir, "$sample.coverage.geneBodyCoverage.r");
my $updated_coverage_r_file = File::Spec->catfile($inputdir, "$sample.coverage.geneBodyCoverage_UPDATED.r");
my $gene_coverage_r_file = File::Spec->catfile($inputdir, "$sample.geneBodyCoverage.r");
my $updated_coverage_r_file = File::Spec->catfile($inputdir, "$sample.geneBodyCoverage_UPDATED.r");

# Modify the R script so that it runs the shapiro-wilk normality test and outputs that test statistic plus p-value to a text file. Also draw the gene body coverage graph with a more distinctly coloured line
open (my $ifh, $gene_coverage_r_file) or die "Could not open file $gene_coverage_r_file $!";
Expand Down

0 comments on commit 502a916

Please sign in to comment.