I''m trying to write a program that can analyze protein relationship. I'm trying to make the user input the accession number of the protein, then save it in a variable, then load the sequence using the getgenpept function, but it doesn't seem to work! (Input AAA69808, NP_989750)
protein1=input('NCBI accession number of the first protein: ', 's');protein2=input('NCBI accession number of the second protein: ', 's');seq1 = getgenpept('protein1','SequenceOnly',true);seq2 = getgenpept('protein2','SequenceOnly',true);
However, doing the folloiwng in MATLAB works fine:
human = getgenpept('AAA69808','SequenceOnly',true);chicken = getgenpept('NP_989750','SequenceOnly',true);
Best Answer