How can I read a DNA sequence (Fasta file ) from NCBI and store it in array or in a specific variable

7 views (last 30 days)
I want to click on command in my window to enter the Accession Number of the sequence, then the output must be for example S=ATCCTGGAA....
i am try this code, but does not work
%% Read sequence data from FASTA files
s = fastaread('AC168908.txt')

Answers (1)

Smit
Smit on 10 Feb 2023
Hi,
I understand you want to read the DNA sequence from a FASTA file.
The function “fastaread” will read the FASTA file and store the data in a structure. This structure also contains the DNA sequence. So, all you need to do extract the sequence access the “Sequence” field in the structure in the following manner
sequence = s.Sequence;
For more information on “fastaread” you can visit the documentation page.

Categories

Find more on Genomics and Next Generation Sequencing in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!