#!/usr/bin/env perl

use Test::Most;

use autodie;
use feature qw(say);

use Path::Class qw(file);

use Bio::MUST::Core;
use Bio::MUST::Core::Utils qw(:tests);

my $class = 'Bio::MUST::Core::SeqMask::Pmsf';

{
    my $s_infile = file('test', 'test-pmsf-archaea.sitefreq');
    my $s_rates = $class->load($s_infile);
    my $o_infile = file('test', 'test-pmsf-euka.sitefreq');
    my $o_rates = $class->load($o_infile);

    my $delta_rates = $s_rates->chi_square_stats($o_rates);
    cmp_store(
        obj    => $delta_rates,
        method => 'store',
        file   => 'test-pmsf-chi-square.stats',
        test   => 'wrote expected chi-square file',
    );
}

done_testing;