I have two gene sequences (1)GAGCAAGCGCCATACTCCTGTGGAGAAGTCTGCCGTTACTGCCCTGTGGGGCAAGGTGAACGTGGA (2)GAGCAAGCGCCATAGTCCTGTGGAGAAGTCTGCCGTTACTGCCCTGTGGGGCAAGGTGAACGTGGA   I want to mark their differences.  Now I marked the differences with lowcase(char) (see my code below). My question is how I can mark the difference with Red color. By the way, I appreciate if someone can optimize my code.   Thanks.       ***Code Start******************************************* data a; length f1 $ 200; input f1; datalines; GAGCAAGCGCCATACTCCTGTGGAGAAGTCTGCCGTTACTGCCCTGTGGGGCAAGGTGAACGTGGA GAGCAAGCGCCATAGTCCTGTGGAGAAGTCTGCCGTTACTGCCCTGTGGGGCAAGGTGAACGTGGA ; data b; set a; retain base len_1; if _n_=1 then do; base=f1; len_1=length(base); end; f2=f1; len_2=length(f2); x=min(len_1, len_2); do i=1 to x; substr_1=substr(base,i,1); substr_2=substr(f2,i,1); if substr_1 ^=substr_2 then substr(f2,i,1)=lowcase(substr_2); else; end; run;   proc print data=b ; var f1 f2; run;   ***Code end***********************************         ... View more

I put the below code. Why are there no data print out?     data ab; input a1-a5 b1-b5; datalines; 1 2 3 4 5 101 102 103 104 105 6 7 8 9 10 106 107 108 109 110 11 12 13 14 15 111 112 113 114 115 ; run; %MACRO render(template, dataset); DATA _NULL_; SET &dataset; FILE PRINT ODS=(TEMPLATE="&template"); PUT ODS; RUN; %MEND render; PROC TEMPLATE; DEFINE TABLE mytable; COLUMN a1-a5 b1-b5; END; RUN; %render(mytable,ab) ... View more