dseq.h

Go to the documentation of this file.

#ifndef DNA_SEQUENCE
#define DNA_SEQUENCE
#include <iostream>
#include "gberr.h"
//file dseq.h
/* to handle the sequence part of genbank sequence */

class dseq
{
        public:
                /* ln may contain the first line or it will read from ins. 
                 * return 0 for normal reading, 1 for sequence composition error
                 */
                int input(istream &ins, char ln[]) throw(gberr); 

                void write(ostream &ous, const int bpl=70) const; //bpl base per line
                dseq();
                dseq(const dseq &ds);
                dseq &operator=(const dseq &ds);
                ~dseq();
                void subseq(char *&sb, int s, int f) const;
                //start index 1 as appose to 0 in C++

        private:
                char *seq;
                int len;             // the correct length of the DNA sequence
                int A, C, G, T, O;
                static const int LL=82;
};

#endif