molecule m;
string filename;
int numstrand, i;
string seq[5000], strandname[5000], type[5000];

filename = "dna.pdb";

getseq_from_pdb( filename, numstrand, seq, strandname, type );

for( i=1; i <= numstrand; i++ ) {
    printf( "strand %s is %8s with sequence: %s\n", strandname[ i ],
                        type[ i ], seq[i] );
}