sequence.h

Go to the documentation of this file.

#ifndef VIENNA_RNA_PACKAGE_SEQUENCE_H
#define VIENNA_RNA_PACKAGE_SEQUENCE_H

typedef struct vrna_sequence_s vrna_seq_t;

typedef struct vrna_alignment_s vrna_msa_t;

#include <ViennaRNA/fold_compound.h>


#define VRNA_SEQUENCE_RNA       1U

#define VRNA_SEQUENCE_DNA       2U

typedef enum {
  VRNA_SEQ_UNKNOWN,   
  VRNA_SEQ_RNA,       
  VRNA_SEQ_DNA        
} vrna_seq_type_e;


struct vrna_sequence_s {
  vrna_seq_type_e type;       
  char            *name;
  char            *string;    
  short           *encoding;  
  short           *encoding5;
  short           *encoding3;
  unsigned int    length;     
};


struct vrna_alignment_s {
  unsigned int  n_seq;
  vrna_seq_t          *sequences;
  char                **gapfree_seq;
  unsigned int        *gapfree_size;  /* for MAF alignment coordinates */
  unsigned long long  *genome_size;     /* for MAF alignment coordinates */
  unsigned long long  *start;           /* for MAF alignment coordinates */
  unsigned char       *orientation;     /* for MAF alignment coordinates */
  unsigned int        **a2s;
};


vrna_seq_t *
vrna_sequence(const char    *string,
              unsigned int  options);


int
vrna_sequence_add(vrna_fold_compound_t  *fc,
                  const char            *string,
                  unsigned int          options);


int
vrna_sequence_remove(vrna_fold_compound_t *fc,
                     unsigned int         i);


void
vrna_sequence_remove_all(vrna_fold_compound_t *fc);


void
vrna_sequence_prepare(vrna_fold_compound_t *fc);


int
vrna_sequence_order_update(vrna_fold_compound_t *fc,
                           const unsigned int   *order);


int
vrna_msa_add( vrna_fold_compound_t      *fc,
              const char                **alignment,
              const char                **names,
              const unsigned char       *orientation,
              const unsigned long long  *start,
              const unsigned long long  *genome_size,
              unsigned int              options);


#endif