<div dir="ltr"><p style="box-sizing: border-box; margin-bottom: 10px;"><span style="color: rgb(51, 51, 51); line-height: 9.285715103149414px; font-family: 'Helvetica Neue', Helvetica, Arial, sans-serif;">I have uploaded genomic fasta sequence and GFF3 file into local mysql database using "bp_seqfeature_load.pl". Now, I want to&nbsp;</span><br>extract sequence from position 8966215-8966961 on "+" strand of Chromosome I, and print only the mRNA segment between 8966215-8966961. If there are many transcripts/isoforms found between these position, print each isoform sequences with name and position as fasta format. Could you please tell me what is the best approach to do it for it. My code does not work well. Thanks<br><br><br>#!/user/bin/perl<br>use strict;<br>use warnings;<br>use Bio::Seq;<br>use Bio::SeqIO;<br>use Bio::SeqFeatureI;<br>use Bio::DB::SeqFeature::Store;<br>use Bio::DB::SeqFeature::Segment;<br><br>my @features;<br>my $db = Bio::DB::SeqFeature::Store-&gt;<br>    new(-adaptor=&gt; 'DBI::mysql',<br>        -dsn=&gt; '..',<br>        -user=&gt; '..',<br>        -pass=&gt; '..',<br>    );<br><br>my @foo = $db-&gt;features(-seq_id =&gt;'I',-start=&gt;'8966215',-end=&gt;'8966961', -strand =&gt;'+1',<br>                        -type =&gt; 'mRNA',<br>                            );<br><br>for my $f(@foo){<br><br>    print "&gt;",$-&gt;name,"_",$f-&gt;start,"_",$f-&gt;end,"\n";<br>    print $f-&gt;dna,"\n";<br><br>}<br><br>exit;<br></p></div>