Mercurial > repos > portiahollyoak > temp

diff scripts/pickUniqPairFastq_MEM.pl @ 21:9672fe07a232 draft default tip

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 0fea84d05f8976b8360a8b4943ecb01b87e3ade0-dirty
author mvdbeek
date Mon, 05 Dec 2016 09:58:47 -0500
parents ca36262102d8
children
line wrap: on
line diff
--- a/scripts/pickUniqPairFastq_MEM.pl	Wed Oct 26 07:24:45 2016 -0400
+++ b/scripts/pickUniqPairFastq_MEM.pl	Mon Dec 05 09:58:47 2016 -0500
@@ -13,7 +13,9 @@
 	chomp;
 	my @f=split/\t/,$_,12;
 	## read number 1 or 2
-	my ($rnum)=$f[1]=~/(\d)$/;
+	#my ($rnum)=$f[1]=~/(\d)$/;
+        my $rnum=1;
+        if (($f[1] & 128) == 128) {$rnum=2;}
 
 	## XT:A:* 
 	my $xt="";
@@ -34,7 +36,7 @@
 	}
 
 	## revcom the read mapped to the reverse strand
-	if($f[1]=~/r/)
+	if (($f[1] & 16) == 16)
 	{
 		my $seq=Bio::Seq->new(-seq=>$f[9], -alphabet => 'dna');
 		$f[9]=$seq->revcom->seq;