How do I download a Reference Genome
4 ビュー (過去 30 日間)
古いコメントを表示
I am going walking through a MATLAB example for Exploring Genome Wide Methylation... The instructions go like this:
This example also assumes that you downloaded the reference human genome (GRCh37.p5). You can use the bowtie-inspect command to reconstruct the human reference directly from the bowtie indices. Or you may download the reference from the NCBI repository by uncommenting the following line:
% getgenbank('NC_000009','FileFormat','fasta','tofile','hsch9.fasta');
How ever this line of code produces an error as shown below:
Error using getncbidata (line 191) The key NC_000009 was not found in the nucleotide database at this time. Check that the input is a valid Accession Number.
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{:});
The accession number seem okay, as I have checked the NCBI website.
Does anyone know how to do this right?
Cheers.
1 件のコメント
Akhilesh Thakur
2017 年 7 月 31 日
https://www.mathworks.com/help/bioinfo/examples/working-with-whole-genome-data.html https://www.mathworks.com/help/bioinfo/examples/comparing-whole-genomes.html#zmw57dd0e13135 Did you check this links?
回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Bioinformatics Toolbox についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!