Mercurial > repos > portiahollyoak > temp
comparison scripts/pickUniqPos_MEM.pl @ 12:ca36262102d8 draft
planemo upload for repository https://github.com/portiahollyoak/Tools commit 5d021f520b653582862ec98dd812a051b804aa50
author | portiahollyoak |
---|---|
date | Fri, 29 Apr 2016 05:47:54 -0400 |
parents | 28d1a6f8143f |
children | 9672fe07a232 |
comparison
equal
deleted
inserted
replaced
11:e19d9742c99b | 12:ca36262102d8 |
---|---|
25 elsif ($a[$i] =~ /^XS:i:/) { | 25 elsif ($a[$i] =~ /^XS:i:/) { |
26 $a[$i] =~ s/XS:i://; | 26 $a[$i] =~ s/XS:i://; |
27 $xs=$a[$i]; | 27 $xs=$a[$i]; |
28 } | 28 } |
29 if (($xs > 0) && ($as-$xs <= $ARGV[1])) {$xt="R";} | 29 if (($xs > 0) && ($as-$xs <= $ARGV[1])) {$xt="R";} |
30 else {$xt="U";} | 30 elsif ($as > 0) {$xt="U";} |
31 } | 31 } |
32 | 32 |
33 my $strand="+"; | 33 my $strand="+"; |
34 ## revcomp | 34 ## revcomp |
35 if($f[1]=~/r/) | 35 if($f[1]=~/r/) |
36 { | 36 { |
37 my $seq=Bio::Seq->new(-seq=>$f[9]); | 37 my $seq=Bio::Seq->new(-seq=>$f[9], -alphabet => 'dna'); |
38 $f[9]=$seq->revcom->seq; | 38 $f[9]=$seq->revcom->seq; |
39 $strand="-"; | 39 $strand="-"; |
40 } | 40 } |
41 | 41 |
42 ## parse CIGAR | 42 ## parse CIGAR |