comparison 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
comparison
equal deleted inserted replaced
10:3200a0299230 11:e19d9742c99b
19 my ($xt)=$f[11]=~/XT:A:(.)/; 19 my ($xt)=$f[11]=~/XT:A:(.)/;
20 20
21 ## revcom the read mapped to the reverse strand 21 ## revcom the read mapped to the reverse strand
22 if($f[1]=~/r/) 22 if($f[1]=~/r/)
23 { 23 {
24 my $seq=Bio::Seq->new(-seq=>$f[9]); 24 my $seq=Bio::Seq->new(-seq=>$f[9], -alphabet => 'dna');
25 $f[9]=$seq->revcom->seq; 25 $f[9]=$seq->revcom->seq;
26 $f[10]=reverse $f[10]; 26 $f[10]=reverse $f[10];
27 } 27 }
28 if (($rnum == 1) || ($rnum == 2)) 28 if (($rnum == 1) || ($rnum == 2))
29 { 29 {