%matrix1=zeros(62,1,844); rotalign=rot90(align); l=length(rotalign); for j=1:74 seq1=(rotalign(j).Sequence); len=length(seq1); for i=1:len if seq1(i)=='R' hydro1([i])=[-7.5]; isoel1([i])=[10.8]; size1([i])=[174]; end if seq1(i)=='K' hydro1([i])=[-4.6]; isoel1([i])=[9.7]; size1([i])=[146]; end if seq1(i)=='D' hydro1([i])=[-3.0]; isoel1([i])=[2.8]; size1([i])=[133]; end if seq1(i)=='Q' hydro1([i])=[-2.9]; isoel1([i])=[5.7]; size1([i])=[146]; end if seq1(i)=='N' hydro1([i])=[-2.7]; isoel1([i])=[5.4]; size1([i])=[132]; end if seq1(i)=='E' hydro1([i])=[-2.6]; isoel1([i])=[3.2]; size1([i])=[147]; end if seq1(i)=='H' hydro1([i])=[-1.7]; isoel1([i])=[7.6]; size1([i])=[155]; end if seq1(i)=='S' hydro1([i])=[-1.1]; isoel1([i])=[5.7]; size1([i])=[105]; end if seq1(i)=='T' hydro1([i])=[-0.8]; isoel1([i])=[5.9]; size1([i])=[119]; end if seq1(i)=='P' hydro1([i])=[-0.3]; isoel1([i])=[6.5]; size1([i])=[115]; end if seq1(i)=='Y' hydro1([i])=[0.1]; isoel1([i])=[5.7]; size1([i])=[181]; end if seq1(i)=='C' hydro1([i])=[0.2]; isoel1([i])=[5.1]; size1([i])=[121]; end if seq1(i)=='G' hydro1([i])=[0.7]; isoel1([i])=[6.0]; size1([i])=[75]; end if seq1(i)=='A' hydro1([i])=[1.0]; isoel1([i])=[6.0]; size1([i])=[89]; end if seq1(i)=='M' hydro1([i])=[1.1]; isoel1([i])=[5.7]; size1([i])=[149]; end if seq1(i)=='W' hydro1([i])=[1.5]; isoel1([i])=[6.0]; size1([i])=[204]; end if seq1(i)=='L' hydro1([i])=[2.2]; isoel1([i])=[6.0]; size1([i])=[131]; end if seq1(i)=='V' hydro1([i])=[2.3]; isoel1([i])=[6.0]; size1([i])=[117]; end if seq1(i)=='F' hydro1([i])=[2.5]; isoel1([i])=[5.5]; size1([i])=[165]; end if seq1(i)=='I' hydro1([i])=[-0.3]; isoel1([i])=[5.9]; size1([i])=[126.9]; end if seq1(i)=='-' hydro1([i])=[0]; isoel1([i])=[0]; size1([i])=[0]; end % prot=vertcat(hydro1, isoel1, size1); % prot=rot90(hydro1); % p=rot90(prot) hydro(j,i)=hydro1(:,i); isoel(j,i)=isoel1(:,i); size(j,i)=size1(:,i); end %matrix1(j,:,:)=[hydro1]; hydrocut=hydro(:,100:200); isoelcut=isoel(:,70:505); sizecut=size(:,70:505); mash=horzcat(hydrocut, isoelcut, sizecut); end