Mercurial > repos > portiahollyoak > temp
diff scripts/pickUniqPairFastq.pl @ 11:e19d9742c99b draft
planemo upload for repository https://github.com/portiahollyoak/Tools commit bfcc36291d97d1c0986f034ca1a03433e8299ff2-dirty
author | portiahollyoak |
---|---|
date | Thu, 28 Apr 2016 07:57:12 -0400 |
parents | 28d1a6f8143f |
children | 9672fe07a232 |
line wrap: on
line diff
--- a/scripts/pickUniqPairFastq.pl Wed Apr 27 11:44:13 2016 -0400 +++ b/scripts/pickUniqPairFastq.pl Thu Apr 28 07:57:12 2016 -0400 @@ -21,7 +21,7 @@ ## revcom the read mapped to the reverse strand if($f[1]=~/r/) { - my $seq=Bio::Seq->new(-seq=>$f[9]); + my $seq=Bio::Seq->new(-seq=>$f[9], -alphabet => 'dna'); $f[9]=$seq->revcom->seq; $f[10]=reverse $f[10]; }