view release on metacpan or  search on metacpan

lib/App/Fasops/Command/stat.pm  view on Meta::CPAN

                next if length $info_of->{ $full_names[0] }{seq} < $opt->{length};
            }

            # outgroup
            my $out_seq;
            if ( $opt->{outgroup} ) {
                $out_seq = pop @{$seq_refs};
            }

            my $first_name  = $full_names[0];
            my $result      = App::Fasops::Common::multi_seq_stat($seq_refs);
            my $indel_sites = App::Fasops::Common::get_indels($seq_refs);

            $csv->print( $out_fh, [ $first_name, @{$result}, scalar( @{$indel_sites} ) ] );
        }
        else {
            $content .= $line;
        }
    }

    close $out_fh;

t/90-common.t  view on Meta::CPAN

        [   [   qw{ GATTATCATCACCCCAGCCACATW
                    GATTTT--TCACTCCATTCGCATA }
            ],
            [ 24, 21, 16, 5, 2, 1, 0, 0.2381, ],
        ],
    );

    for my $i ( 0 .. $#data ) {
        my ( $seq_refs, $except_ref ) = @{ $data[$i] };

        my $result_ref = App::Fasops::Common::multi_seq_stat($seq_refs);
        Test::Number::Delta::delta_ok( $result_ref, $except_ref, "stat $i" );
    }
}

{
    print "#ref_pair_D\n";

    my @data = (

        #