%read two sequences
%%getgenbank('X00166', 'TOFILE', 'HGENBANKX00166.GBK');
%%getgenbank('V01153', 'TOFILE', 'HGENBANKV01153.GBK');

seq1 = genbankread('HGENBANKX00166.GBK');
seq2 = genbankread('HGENBANKV01153.GBK');

window=11; num=7;
seqdotplot(seq1,seq2,window,num)
xlabel('X00166');
ylabel('V01153');
title('Window 11 Num 7');

%%getgenpept('CAA24991', 'TOFILE', 'HGENBANKCAA24991.GBK');
%%getgenpept('CAA24470', 'TOFILE', 'HGENBANKCAA24470.GBK');

seq3 = genbankread('HGENBANKCAA24991.GBK');
seq4 = genbankread('HGENBANKCAA24470.GBK');

window=10; num=3;
seqdotplot(seq3,seq4,window,num)
xlabel('CAA24991');
ylabel('CAA24470');
title('Window 10 Num 3');

%%getgenpept('P01130', 'TOFILE', 'HGENBANKP01130.GBK');
seq5 = genbankread('HGENBANKP01130.GBK');

window=1; num=1;
seqdotplot(seq5,seq5,window,num)
xlabel('P01130 Human LDL receptor');
ylabel('P01130 Human LDL receptor');
title('Window 1 Num 1');

window=23; num=7;
seqdotplot(seq5,seq5,window,num)
xlabel('P01130 Human LDL receptor');
ylabel('P01130 Human LDL receptor');
title('Window 23 Num 7');