[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