Mercurial > repos > mcharles > rapsodyn
view rapsodyn/filtersam_mapped_and_unique.pl @ 28:fefe6d89dbf1 draft
Uploaded
author | mcharles |
---|---|
date | Wed, 10 Sep 2014 08:10:38 -0400 |
parents | f487491e1091 |
children | 7b8646f46010 |
line wrap: on
line source
#!/usr/bin/perl use strict; use warnings; open(IN, $ARGV[0]) or die ("Can't open $ARGV[0]\n"); while (my $line=<IN>){ if ($line =~ /^\@/){ #Header conservation print $line; } else { #Optionnal flag verification if (($line =~ /XT\:A\:U/)&&($line =~ /X0\:i\:1/)&&($line =~ /X1\:i\:0\s/)){ my @fields = split (/\s+/,$line); if (($fields[1]==83)||($fields[1]==163)||($fields[1]==147)||($fields[1]==99)){ print $line; } } } } close (IN);