1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
   |  
int get_trinucls(char* sequence, int position_start, int position_end, char* trinucl_sens, char* trinucl_anti, int cut_rank) {
	int j;
	char base_anti, base_sens;
	for (j=0; j<3; j++) {
		base_anti = sequence[position_end-cut_rank-j-1];
		base_sens = sequence[position_start+j+(ENZYME_LEN-cut_rank)-1];
		if( (base_anti == 'a' || base_anti == 'c' || base_anti == 'g' || base_anti == 't') && (base_sens == 'a' || base_sens == 'c' || base_sens == 'g' || base_sens == 't') ) {
			trinucl_sens[j] = base_sens;
// 			if (j == 0) printf("%c", trinucl_sens[j]);
			switch (base_anti) {
				case 'a' : base_anti = 't'; break;
				case 'c' : base_anti = 'g'; break;
				case 'g' : base_anti = 'c'; break;
				case 't' : base_anti = 'a'; break;
			}
			trinucl_anti[j] = base_anti;
		}
		else {
			return 0;
		}
	}
	trinucl_sens[3]='\0';
	trinucl_anti[3]='\0';
// 	printf("\n");
//      et là j'ai exactement les memes caractères !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	printf("%c\t%c\n", trinucl_anti[4], trinucl_sens[1]);
// 	printf("%s\t%u\t%s\t%u\t%s\n", psub->name, position_start, trinucl_sens, position_end, trinucl_anti);
	return 1;
} | 
Partager