Hi Can, > Can I get some sample code using the Seq class? $seqo=new seq; // if we are on the - strand, reverse complement if ($strand=='minus') { $sequence=revcomp($sequence,'DNA'); } $seqo->attribute('sequence',$sequence); $seqo->attribute('seqlength',strlen($sequence)); $seqo->attribute('moltype','DNA'); $proto=new seq; $proto->sequence=$seqo->translate(0,1); $proto->attribute('moltype','PROTEIN'); > I wasn't sure how to use it since I couldn't find a constructor. I guess there is none at the moment... See the directory testfiles for examples of using the seqIOimport class (which gives you a parser that can return a seqobject). Best, Nico > Thanks, > Can > > > _______________________________________________ > Biophp-dev mailing list > Biophp-dev@bioinformatics.org > https://bioinformatics.org/mailman/listinfo/biophp-dev >