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
	p=index(ed,mutb);
	if (p>0){
	    c=substr(ed,p-1,1);
	    if (0+c==c){
		pos=0;x=0;
		for (i=1;i<=p;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;}
			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==refb) {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        
    
 }