BEGIN{FS=OFS="\t";}{
    if (substr($1,1,1)==">"){
	rd=0;
	st=$3;refs=$4;refe=$5;rds=$6;rde=$7;sc=$9;ed=$10;rd=substr($1,2);ri=$2;shr=$0;
	if (st=="+"){
	    refb="T";mutb="C";
	}else{
	    refb="A";mutb="G";
	}
	#loop over multiple Ts/As
	el=length(ed);nrefT=0;
	for (i=1;i<=el;i++){
	    c=substr(ed,i,1);
	    if (c==refb){
		c2=substr(ed,i-1,1);
		if (0+c2==c2){
		    refTpos[++nrefT]=i;
		}
	    }
	}
#	p=index(ed,refb);
#	    c=substr(ed,p-1,1);
#	    if (0+c==c){
	if (nrefT>0){
		pos=0;x=0;
		for (i=1;i<=el;i++){
		    c=substr(ed,i,1);
		    if (0+c==c){r=1;
			x=x*10+c;#print i" "c" "x"n";
		    }else{
			if(r==1) {pos+=x;}			
			for (pi=1;pi<=nrefT;pi++){
			    if (refTpos[pi]==i) {refTpos2[pi]=pos;}
			}
			if (i<p){
			    if ((c!="(")&&(c!=")")&&(c!="x")){pos++;}
			}
#			if (c=="-"){pos++;}
#			if ((c=="x")||(c=="X")){pos++;}
#			print i" "c" "pos" c";			
			r=0;x=0;
		    }
		}
		getline;
		n=split($1,a," ");
		mutpos=substr(a[3],pos+rds,1);
		if (mutpos==mutb) {rd=1;
		    print mutpos" pos"pos" "pos+rds" "shr;
		    print $0;
		    ntc++;
		}else{nntc++;}
	    }
	}
    }else{
	if(rd==1){print $0}
    }
}END{
    print "nTC "ntc"  nonTC "nntc"  TCrate "100.0*ntc/(ntc+nntc)
#- G AGCCCTCCCGGCCTACAGCGAGCTCTCGCCAATAGACACTGTGTCTCATC 6x5x5x1C5x1x2x2A3A3 pos27 41
#G:       2772    AGCCCTCCCGGCCTACAGCGAGCTCTCGCCAATAGACACTGTGTCTCATC    2737      
#                              ||||||X||||X||||X |||||XX|X| ||| |||
#T:               -------------TACAGCgAGCTcTCGCcCATAGAcaCtGAGTCACAT-
#R:         14   T00233021113031311231223202233201332203131221211131    49        
#+ C TGCGGATGGACACCTTGGCCCATGT--TCTCTACTATCCTCAGAAACCCC 1x9x5(CG)3x2T5 pos20 31
#G:       2700    TGCGGATGGACACCTTGGCCCATGT--TCTCTACTATCCTCAGAAACCCC    2725      
#                           |X||||||||X||||  |||X| |||||
#T:               ----------CaCCTTGGCCcATGTCGTCTcTTCTATC------------
#R:         11   T31013132112310201030313112312232022332212221111320    38        
    
 }