Skip to content

Commit

Permalink
Merge pull request #107 from BioInf-Wuerzburg/hotfix/ccseq-phred-offset
Browse files Browse the repository at this point in the history
ccseq: assume phred-offset 33 for read if guessing fails. Thanks to David Mathog
  • Loading branch information
thackl authored Sep 19, 2017
2 parents 4c17c5e + da6cf11 commit 21fecc0
Show file tree
Hide file tree
Showing 2 changed files with 5 additions and 1 deletion.
4 changes: 4 additions & 0 deletions bin/ccseq
Original file line number Diff line number Diff line change
Expand Up @@ -198,6 +198,10 @@ $L->info("Reading STDIN") unless $opt{reads};
my $rpr = Fastq::Parser->new($opt{reads} ? (file => $opt{reads}) : ())->check_format;
if ($rpr) {
my $po = $rpr->guess_phred_offset();
if(!defined($po)){
$po=33;
$L->info("Couldn't guess phred offset, but assuming 33 for PacBio");
}
$L->logdie("Expected phred offset 33") if $po && $po != 33;
$rpr->phred_offset($po); # should be 33 for PacBio
} else {
Expand Down
2 changes: 1 addition & 1 deletion bin/proovread
Original file line number Diff line number Diff line change
Expand Up @@ -41,7 +41,7 @@ use Time::HiRes qw(usleep);
use Data::Dumper;
$Data::Dumper::Sortkeys = 1;

use version 0.77; our $VERSION = qv("2.14.0");
use version 0.77; our $VERSION = qv("2.14.1");

my $RealLib = "$RealBin/../lib/";
my $RealPerl = "perl -I${RealLib} $RealBin";
Expand Down

0 comments on commit 21fecc0

Please sign in to comment.