Mercurial > repos > portiahollyoak > temp

diff scripts/pickUniqPairFastq.pl @ 11:e19d9742c99b draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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];
 	}