how to import gene fasta file from NCBI using matlab 2016a

Question

Priyanka Roy 2017-2-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/327369-how-to-import-gene-fasta-file-from-ncbi-using-matlab-2016a

评论： Paola Favaretto 2017-3-1

I want to import gene fasta file from NCBI database using the Accession Number. I use the code Data = getgenbank('NP_752927.1'); but, getting this error : Error using getncbidata (line 191) The key NP_752927.1 was not found in the nucleotide database at this time. Please check that the input is a valid accession number or try again.

NOTE: This function is dependent on NCBI's Entrez tools and sequence databases. Changes to either may cause this function to break.

Error in getgenbank (line 70)

    [varargout{1:nargout}] = getncbidata(accessnum,'fileformat','GenBank','database','nucleotide',varargin{:});

How will i resolve the error?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Paola Favaretto 2017-2-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/327369-how-to-import-gene-fasta-file-from-ncbi-using-matlab-2016a#answer_256718

Are you sure it is a valid accession number? When I search the NCBI databases with the id you provided, I get 0 results.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Priyanka Roy 2017-3-1

Escherichia coli O157:H7 str. Sakai, complete genome. NCBI Reference Sequence: NC_002695.1 ACCESSION NC_002695

I can manually download the Fasta file from NCBI database, but getting the same error while using the code : Data = getgenbank('NC_002695');

what will be the code for importing the fasta file directly from NCBI database ?

请先登录，再进行评论。

Answer 2

Paola Favaretto 2017-3-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/327369-how-to-import-gene-fasta-file-from-ncbi-using-matlab-2016a#answer_256847

编辑：Paola Favaretto 2017-3-1

在 MATLAB Online 中打开

What version of Bioinformatics Toolbox are you using? I am able to download the sequence without issues.

You can get the sequence information by typing:

a = getgenbank('NC_002695', 'sequenceonly', true)

Or you can save the sequence in a FASTA formatted file by typing:

a = getgenbank('NC_002695', 'tofile', 'S:/myfile2.fa', 'fileformat', 'fasta')

This is a snippet of the output:

a = 
    struct with fields:
        Header: 'NC_002695.1 Escherichia coli O157:H7 str. Sakai, complete genome'
      Sequence: 'AGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTCTCTGACAGC ...'