use strict;
use warnings;
my %ASSEMBLY_TESTS;

BEGIN {
    use lib '.';
    use Bio::Root::Test;

    test_begin( -tests => 796,
                -requires_modules => [
                    'DB_File',
                    'Bio::Tools::Run::Samtools',
                   ],
               );

    use_ok('Bio::Seq');
    use_ok('Bio::LocatableSeq');
    use_ok('Bio::Seq::Quality');
    use_ok('Bio::Assembly::IO');
    use_ok('Bio::Assembly::Singlet');

}

use Bio::Root::IO;

my ($aio, $assembly, @contig_seq_ids, @singlet_ids, @contig_ids, @all_seq_ids);

my $file = 'test.bowtie';
my $refdb = 'test.ref.fas';
ok $aio = Bio::Assembly::IO->new( -file => test_input_file($file),
                  -index => test_input_file($refdb),
                  -format => 'bowtie' ), "init bowtie IO object";
isa_ok($aio, 'Bio::Assembly::IO');
$aio->_current_refseq_id( ($aio->sam->seq_ids)[0] ); # kludge

while (my $contig = $aio->next_contig) {
isa_ok($contig, 'Bio::Assembly::Contig');
}
ok $aio = Bio::Assembly::IO->new( -file => test_input_file($file),
                  -index => test_input_file($refdb),
                  -format => 'bowtie' ),"reopen";
ok $assembly = $aio->next_assembly, "get sam assy";
is( $assembly->get_nof_contigs, 23, "got all contigs");

ok(@contig_seq_ids = $assembly->get_contig_seq_ids, "get_contig_seq_ids");
is(@contig_seq_ids, 312);
for my $contig_seq_id (@contig_seq_ids) {
ok ($contig_seq_id =~ m/^SRR/i);
}
ok(@contig_ids = $assembly->get_contig_ids, "get_contig_ids");
is(@contig_ids, 23);
for my $contig_id (@contig_ids) {
ok ($contig_id =~ m/sam_assy/i);
}
ok(@singlet_ids = $assembly->get_singlet_ids, "get_singlet_ids");
is(@singlet_ids, 36);
for my $singlet_id (@singlet_ids) {
ok ($singlet_id =~ m/^sam_assy/i);
}
ok(@all_seq_ids = $assembly->get_all_seq_ids, "get_all_seq_ids");
for my $seq_id (@all_seq_ids) {
ok ($seq_id =~ m/^SRR/i);
}
is(@all_seq_ids, 348);

exit;