void revcomp_dna( char *seq )
{
/* Martin A. Hansen, May 2008 */
/* Reverse complement a DNA sequence in place. */
complement_dna( seq );
reverse( seq );
}
void seqed_write_complement(char *sequence,
int pos,
int line_length,
char *line)
{
if (line_length <= 0)
return;
strncpy(line, sequence, line_length);
line[line_length] = '\0';
complement_dna(line, line_length);
}
void complement_nuc( char *seq )
{
/* Martin A. Hansen, May 2008 */
/* Complements a nucleotide sequence, */
/* after guess the type. */
if ( is_dna( seq ) ) {
complement_dna( seq );
} else if ( is_rna( seq ) ) {
complement_rna( seq );
} else {
abort();
}
}
