[Biojava-l] bug using getFeature() ?

Alberto Ambesi ambesi at tigem.it
Thu Jul 31 19:55:18 EDT 2003


my program fails to get a transcript feature on human chromosome 1  
within positions: 156947999-156952998. If the search is done via web at  
ensembl.org, I find ensembl transcript CASQ1, with an exon in position:  
156948909-156949040. In many other genomic regions the program works  
well. Why is this?
thank you in advance.
dr. Alberto Ambesi

--

import java.sql.*;
import java.util.*;
import org.ensembl.db.*;
import org.ensembl.db.sql.SQLDatabaseAdaptor;
import org.biojava.bio.*;
import org.biojava.bio.program.ensembl.*;
import org.biojava.bio.seq.*;
import org.biojava.bio.seq.db.*;
import org.biojava.bio.seq.io.*;
import org.biojava.bio.symbol.*;

public class bug {
     public static void main(String[] args) throws Exception {
         Driver driverName =  
(Driver)Class.forName("org.gjt.mm.mysql.Driver").newInstance();
         //Connection conn =  
DriverManager.getConnection("jdbc:mysql://kaka.sanger.ac.uk/ 
homo_sapiens_core_15_33", "anonymous", "");
         SQLDatabaseAdaptor sqlDBAdaptor =  
SQLDatabaseAdaptor.connectSQL("jdbc:mysql://kaka.sanger.ac.uk/ 
homo_sapiens_core_15_33", "anonymous", "", 9);
         Ensembl ens = new Ensembl(sqlDBAdaptor);
         SequenceDB chromosomes =ens.getChromosomes();
         Sequence chrSeq = chromosomes.getSequence("1");
         Sequence seq = new SubSequence(chrSeq, 156947999, 156952998);
         System.out.println(seq.seqString());
         for (Iterator fi = seq.features(); fi.hasNext(); ) {
             org.biojava.bio.seq.Feature f =  
(org.biojava.bio.seq.Feature) fi.next();
             String featureType = f.getType();
             System.out.println(featureType);
         }
     }
}
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: text/enriched
Size: 2985 bytes
Desc: not available
Url : http://pw600a.bioperl.org/pipermail/biojava-l/attachments/20030731/fbc795c8/attachment.bin


More information about the Biojava-l mailing list