[ensembl-dev] Problem with biomaRt::getSequence

Mohammad Tanvir Ahamed mashranga at yahoo.com
Wed May 8 12:16:37 BST 2013


Hi,
I can run the code some days ago . But cant run now. 

Problem 1: Output is ok
ensembl = useDataset("hsapiens_gene_ensembl",mart=ensembl)
utr5 = getSequence(chromosome=3, start=185514033, end=185535839, type="entrezgene",seqType="5utr", mart=ensembl) 
Output : 
                                                                                              5utr  entrezgene
                                                                             Sequence unavailable      10644
                                             GGAGCGCCGGGTACCGGGCCGGGGGAGCCGCGGGCTCTCGGGGAAGAGACGG      10644
3 GGGGGGCGGAGGAGGAGGAGAGACGAGGGCAGCGGAGGAGGCGAGGAGCGCCGGGTACCGGGCCGGGGGAGCCGCGGGCTCTCGGGGAAGAGACGG      10644
                                 CGGAGGAGGCGAGGAGCGCCGGGTACCGGGCCGGGGGAGCCGCGGGCTCTCGGGGAAGAGACGG      10644
                                                         No UTR is annotated for this transcript      10644
 
Problem 2:Problem is here

protein = getSequence(id=c(100, 5728),type="entrezgene",seqType="peptide", mart=ensembl)

Error in getBM(c(seqType, type), filters = type, values = id, mart = mart,  : 
  Query ERROR: caught BioMart::Exception::Database: Error during query execution: Can't create/write to file '/mnt/ephemeral0/mysqltmp/#sql_40a_0.MYI' (Errcode: 2)

I need help please.
 
/.......Tanvir Ahamed
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20130508/de2c792a/attachment.html>


More information about the Dev mailing list