Skip to content
Snippets Groups Projects
Commit a21804fc authored by premanand17's avatar premanand17
Browse files

Added transl_table qualifier to genbank and embl flat file dumps - ENSCORESW-1971

parent a1098db4
No related branches found
No related tags found
2 merge requests!168Enscoresw 1971,!168Enscoresw 1971
......@@ -105,7 +105,20 @@ my $index_count_fh = sub {
is(scalar(@{$lines}), 1, 'Expect only 1 EMBL SQ line describing a sequence');
is($lines->[0], 'SQ Sequence 100001 BP; 24986 A; 24316 C; 24224 G; 26475 T; 0 other;', 'Formatting of SQ as expected');
}
# check if transl_table is included
$sd = Bio::EnsEMBL::Utils::SeqDumper->new();
$sd->{feature_types}->{$_} = 0 for keys %{$sd->{feature_types}};
$sd->{feature_types}->{'gene'} = 1;
{
my $mt_slice = $slice_adaptor->fetch_by_region('chromosome', 'MT_NC_001807', 10060, 10405);
my $fh = IO::String->new();
$sd->dump_embl($mt_slice, $fh);
my $lines = $index_fh->($fh, 'FT ');
like( $lines->[9], qr/FT\s+\/transl_table=2/, "Expected transl_table line at FT CDS");
}
{
my $fh = IO::String->new();
$sd->dump_genbank($slice, $fh);
......
......@@ -168,6 +168,6 @@
162030 21739 1
162031 21739 4
162032 21739 3
162033 21740 0
162033 21740 1
162034 21741 1
162035 21741 2
\ No newline at end of file
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment