X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bp_bin%2Fbowtie_seq;h=07374216593d0a6f24d7fd858f6de3b21c6377ec;hb=5de6112b70b59420b245ce636a8b2e3c90acbe00;hp=57b673117c427524b8df7b8ca61aaf3506a9511f;hpb=12bb05ba3197c8e0ccb478d4b634f5555da14752;p=biopieces.git diff --git a/bp_bin/bowtie_seq b/bp_bin/bowtie_seq index 57b6731..0737421 100755 --- a/bp_bin/bowtie_seq +++ b/bp_bin/bowtie_seq @@ -80,7 +80,7 @@ while ( $record = Maasha::Biopieces::get_record( $in ) ) if ( $entry = Maasha::Fastq::biopiece2fastq( $record ) ) { Maasha::Common::error( "Mixed FASTA and FASTQ entries in stream" ) if defined $type and $type ne "FASTQ"; - Maasha::Common::error( "Sequence longer than 1024 not allowed") if length( $entry[ SEQ ] ) > 1024; + Maasha::Common::error( "Sequence longer than 1024 not allowed") if length( $entry->[ SEQ ] ) > 1024; Maasha::Fastq::put_entry( $entry, $fh_out ); $type = "FASTQ"; @@ -88,7 +88,7 @@ while ( $record = Maasha::Biopieces::get_record( $in ) ) elsif ( $entry = Maasha::Fasta::biopiece2fasta( $record ) ) { Maasha::Common::error( "Mixed FASTA and FASTQ entries in stream" ) if defined $type and $type ne "FASTA"; - Maasha::Common::error( "Sequence longer than 1024 not allowed") if length( $entry[ SEQ ] ) > 1024; + Maasha::Common::error( "Sequence longer than 1024 not allowed") if length( $entry->[ SEQ ] ) > 1024; Maasha::Fasta::put_entry( $entry, $fh_out ); $type = "FASTA"; @@ -103,7 +103,7 @@ push @args, "-n $options->{ 'mismatches' }"; push @args, "-v $options->{ 'mismatches' }"; # DANGER: using seed mismatches as alignment mismatches - may work, may not! push @args, "-f" if $type eq "FASTA"; push @args, "-p $options->{ 'cpus' }"; -push @args, "--phred64-quals" unless $type eq "FASTA"; +push @args, "--phred33-quals" unless $type eq "FASTA"; if ( defined $options->{ 'max_hits' } ) { push @args, "-k $options->{ 'max_hits' }";