Mercurial > repos > portiahollyoak > temp
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 { |