[Biopython] gff3 problem
Michal
mictadlo at gmail.com
Fri Apr 8 06:35:33 UTC 2011
On Thu, Apr 7, 2011 at 10:54 PM, Michal <mictadlo at gmail.com> wrote:
> On 04/05/2011 11:22 PM, Brad Chapman wrote
>
>> in_handle = open(in_file)
>> for rec in GFF.parse(in_handle):
>> for feature in rec.features:
>> print feature.type, feature.location
>> print feature.qualifiers
>> for sub_feature in feature.sub_features:
>> print " ", sub_feature.type, sub_feature.location
>> in_handle.close()
>>
>> This will print out details of each feature. For instance, here is
>> a gene with exon sub_features:
>>
>> gene [2234:3344]
>> {'Note': ['Elongation factor P (EF-P) family protein n:2 Tax:Arabidopsis
>> RepID:D7L774_ARALY'],
>> 'source': ['x'], 'ID': ['BC-x.1'], 'Name': ['BC-x.1']}
>> exon [2234:2279]
>> exon [2422:2535]
>> exon [2609:2691]
>> exon [2762:2864]
>> exon [2971:3049]
>> exon [3125:3251]
>> exon [3320:3344]
>>
>> Hope this helps,
>> Brad
>>
>> Thank you it works.
>
How could I get also the strand position?
More information about the Biopython
mailing list