[Biojava-dev] Genbank Feature extraction question.

[michael walsh]

BioJava,

I am writing a program that extracts features from a Genbank file using
BioJava The program needs to extract the feature locations from the file.
This is easy enough to do but the location information returned by the
program does not specify whether or not the location is on the complementary
strand of DNA.

This is some of my code:
FeatureHolder fltrHold = seq.filter(codeFltr);//codeFltr is a filter that
retrieves coding features only.
        //iterate over the Features in fh
        for (Iterator i = fltrHold.features(); i.hasNext(); ){
            Feature f = (Feature)i.next();
            System.out.println(f.getLocation().toString());

An example of the output of the print command is:
join:[32775..32948,31801..32052]

However, the entry in the Genbank file reads:
complement(join(31801..32052,32775..32948))

Is there any way to get my program to output the fact that a features
location is on the complementary strand of DNA?

Any help that anyone can give me would be greatly appreciated.

Sincerely,

M Walsh