root/branches/stable_5.0/ARB_GDE/GDE_FileIO.cxx

#include <stdio.h>
#include <stdlib.h>
// #include <malloc.h>
#include <string.h>

#include <arbdb.h>
#include <arbdbt.h>
#include <aw_root.hxx>
#include <aw_device.hxx>
#include <aw_window.hxx>
#include <awt.hxx>

#include "gde.hxx"
#include "GDE_def.h"
#include "GDE_menu.h"
#include "GDE_extglob.h"

int MAX(int a,int b)
{
    if(a>b) return a;
    return b;
}

int MIN(int a,int b)
{
    if(a<b) return a;
    return b;
}

void Regroup(NA_Alignment *alignment)
{
    size_t j;
    size_t group;
    int last;

    for(j=0;j<alignment->numelements;j++)
    {
        alignment->element[j].groupf = NULL;
        alignment->element[j].groupb = NULL;
    }

    for (group = 1;group <= alignment->numgroups;group++)
    {
        last = -1;
        for(j=0;j<alignment->numelements;j++)
            if(alignment->element[j].groupid == group)
            {
                if(last != -1)
                {
                    alignment->element[j].groupb =
                        &(alignment->element[last]);
                    alignment->element[last].groupf =
                        &(alignment->element[j]);
                }
                last = j;
            }
    }
    return;
}


/*
 *   Print error message, and die
 */
void ErrorOut5(int code,const char *string)
{
    if (code == 0)
    {
        fprintf(stderr,"Error:%s\n",string);
        exit(1);
    }
    return;
}


/*
 *   More robust memory management routines
 */
char *Calloc(int count,int size)
{
    char *temp;
    size *= count;
#ifdef SeeAlloc
    extern int TotalCalloc;
    TotalCalloc += count*size;
    fprintf(stderr,"Calloc %d %d\n",count*size,TotalCalloc);
#endif
    temp = (char *)malloc(size);
    ErrorOut5(0!= temp,"Cannot allocate memory");
    memset(temp,0,size);
    return(temp);
}

char *Realloc(char *block,int size)
{
    char       *temp;
#ifdef SeeAlloc
    extern int  TotalRealloc;
    TotalRealloc += size;
    fprintf(stderr,"Realloc %d\n",TotalRealloc);
#endif
    temp          = (char *)realloc(block,size);
    ErrorOut5(0   != temp,"Cannot change memory size");

    return(temp);
}

void Cfree(char *block)
{
    if (block)
    {
        /*if(cfree(block) == 0)
          Warning("Error in Cfree...");*/
        free(block);
    }
    else
        Warning("Error in Cfree, NULL block");
    return;
}


static void ReadNA_Flat(char *filename,char *dataset,int type)
{
    size_t j;
    int i, jj, c, curelem=0,offset;
    char buffer[GBUFSIZ];
    char in_line[GBUFSIZ];
    char curname[GBUFSIZ];
    i=0;c=0;type=0;

    NA_Sequence *this_elem;
    NA_Alignment *data;

    FILE *file;

    curname[0] = '\0';
    data = (NA_Alignment*)dataset;

    file = fopen(filename,"r");
    if(file == NULL)
    {
        fprintf(stderr,"Cannot open %s.\n",filename);
        return;
    }
    for(;fgets(in_line,GBUFSIZ,file) !=0;)
    {
        if (in_line[0] == '#' ||
            in_line[0] == '%' ||
            in_line[0] == '"' ||
            in_line[0] == '@')
        {
            offset = 0;
            for(j=0;j<strlen(in_line);j++)
            {
                if(in_line[j] == '(')
                {
                    sscanf((char*)
                           &(in_line[j+1]),"%d",&offset);
                    in_line[j] = '\0';
                }
            }

            curelem = data->numelements++;
            if( curelem == 0 )
            {
                data->element=(NA_Sequence*)
                    Calloc(5,sizeof(NA_Sequence));
                data->maxnumelements = 5;
            }
            else if (curelem==data->maxnumelements)
            {
                (data->maxnumelements) *= 2;
                data->element=
                    (NA_Sequence*)Realloc((char*)data->element
                                          ,data->maxnumelements*sizeof(NA_Sequence));
            }

            InitNASeq(&(data->element[curelem]),
                      in_line[0] == '#'?DNA:
                      in_line[0] == '%'?PROTEIN:
                      in_line[0] == '"'?TEXT:
                      in_line[0] == '@'?MASK:TEXT);
            this_elem= &(data->element[curelem]);
            if(in_line[strlen(in_line)-1] == '\n')
                in_line[strlen(in_line)-1] = '\0';
            strncpy(this_elem->short_name,(char*)&(in_line[1]),31);
            this_elem->offset = offset;
        }
        else if(in_line[0] != '\n')
        {
            size_t strl = strlen(in_line);
            for(j=0,jj=0;j<strl;j++)
                if(in_line[j] != ' ' && in_line[j] != '\n' &&
                   in_line[j] != '\t')
                    buffer[jj++] = in_line[j];

            if(data->element[curelem].rmatrix)
                Ascii2NA(buffer,jj,data->element[curelem].rmatrix);
            AppendNA((NA_Base*)buffer,jj,&(data->element[curelem]));
        }
    }

    for(j=0;j<data->numelements;j++)
        data->maxlen = MAX(data->maxlen,data->element[j].seqlen +
                           data->element[j].offset);

    for(j=0;j<data->numelements;j++)
        if(data->element[j].seqlen==0)
            data->element[j].protect =
                PROT_BASE_CHANGES+ PROT_GREY_SPACE+
                PROT_WHITE_SPACE+ PROT_TRANSLATION;

    NormalizeOffset(data);
    Regroup(data);
    return;
}

/*
  LoadFile():
  Load the given filename into the given dataset.  Handle any
  type conversion needed to get the data into the specified data type.
  This routine is used in situations where the format and datatype is known.

  Copyright (c) 1989-1990, University of Illinois board of trustees.  All
  rights reserved.  Written by Steven Smith at the Center for Prokaryote Genome
  Analysis.  Design and implementation guidance by Dr. Gary Olsen and Dr.
  Carl Woese.

  Copyright (c) 1990,1991,1992 Steven Smith at the Harvard Genome Laboratory.
  All rights reserved.
*/

static void LoadFile(char *filename,NA_Alignment *dataset,int type,int format)
{

    if (DataType != type)
        fprintf(stderr,"Warning, datatypes do not match.\n");
    /*
      Handle the overwrite/create/merge dialog here.
    */
    switch(format)
    {
        case NA_FLAT:
            ReadNA_Flat(filename,(char*)dataset,type);
            ((NA_Alignment*)dataset)->format = GDE;
            break;

        case GENBANK:
            ReadGen(filename,dataset,type);
            ((NA_Alignment*)dataset)->format = GENBANK;
            break;

        case ARBDB:
            ReadArbdb_plain(filename,dataset,type);
            ((NA_Alignment*)dataset)->format = ARBDB;
            break;

        case GDE:
            ReadGDE(filename,dataset,type);
            ((NA_Alignment*)dataset)->format = GDE;
            break;
        case COLORMASK:
            ReadCMask(filename);

        default:
            break;
    }
    return;
}

static int FindType(char *name,int *dtype,int *ftype)
{
    FILE *file;
    char in_line[GBUFSIZ];

    file = fopen(name,"r");
    *dtype=0;
    *ftype=0;

    if (file == NULL)
        return(1);

    /*
     *   Is this a flat file?
     *   Get the first non blank line, see if a type marker shows up.
     */
    if (fgets(in_line,GBUFSIZ,file) == 0) {
        return 1;
    }
    for(;strlen(in_line)<2 && fgets(in_line,GBUFSIZ,file) != NULL;) ;

    if (in_line[0] == '#' || in_line[0] == '%' ||
       in_line[0]  == '"' || in_line[0] == '@' )
    {
        *dtype=NASEQ_ALIGN;
        *ftype=NA_FLAT;
    }

    /*
     *   Else, try genbank
     */
    else
    {
        fclose(file);
        file = fopen(name,"r");
        *dtype=0;
        *ftype=0;

        if (file == NULL)
            return(1);

        for(;fgets(in_line,GBUFSIZ,file) != NULL;)
            if(Find(in_line,"LOCUS"))
            {
                *dtype=NASEQ_ALIGN;
                *ftype=GENBANK;
                fclose(file);
                return(0);
            }
        /*
         *   and last, try GDE
         */
            else if(Find(in_line,"sequence"))
            {
                *dtype = NASEQ_ALIGN;
                *ftype = GDE;
                fclose(file);
                return(0);
            }
            else if(Find(in_line,"start:"))
            {
                *dtype = NASEQ_ALIGN;
                *ftype = COLORMASK;
                fclose(file);
                return(0);
            }
    }

    fclose(file);
    return(0);
}

/*
  LoadData():
  Load a data set from the command line argument.

  Copyright (c) 1989, University of Illinois board of trustees.  All rights
  reserved.  Written by Steven Smith at the Center for Prokaryote Genome
  Analysis.  Design and implementation guidance by Dr. Gary Olsen and Dr.
  Carl Woese.

  Copyright (c) 1990,1991,1992 Steven Smith at the Harvard Genome Laboratory.
  All rights reserved.

*/

void LoadData(char *filen)
{

    FILE *file;
    NA_Alignment *DataNaAln;
    char temp[1024];
    /*
     *   Get file name, determine the file type, and away we go..
     */
    if(Find2(filen,"gde")!=0)
        strcpy(FileName,filen);

    if (strstr(filen,".arb") || strchr(filen, ':')) {   /* ARBDB TYPE */
        if (DataSet == NULL) {
            DataSet = (NA_Alignment *) Calloc(1,
                                              sizeof(NA_Alignment));
            DataNaAln = (NA_Alignment *) DataSet;
            DataSet->rel_offset = 0;
        } else{
            DataNaAln = (NA_Alignment *) DataSet;
        }
        DataType = NASEQ_ALIGN;
        FileFormat = ARBDB;
        LoadFile(filen, DataNaAln,
                 DataType, FileFormat);

        sprintf(temp,"Remote ARBDB access (%s)",filen);
        return;
    }


    if( (file=fopen(filen,"r"))!=0 )
    {
        FindType(filen,&DataType,&FileFormat);
        switch(DataType)
        {
            case NASEQ_ALIGN:
                if(DataSet == NULL)
                {
                    DataSet = (NA_Alignment*)Calloc(1,
                                                    sizeof(NA_Alignment));
                    DataNaAln =(NA_Alignment*)DataSet;
                    DataSet->rel_offset = 0;
                }else{
                    DataNaAln = (NA_Alignment*)DataSet;
                }

                LoadFile(filen,DataNaAln,
                         DataType,FileFormat);

                break;
            default:
                aw_message(GBS_global_string("Internal error: unknown file type of file %s",filen));
                break;
        }
        fclose(file);
    }
    sprintf(temp,"Genetic Data Environment 2.2 (%s)",FileName);
    return;
}


void AppendNA(NA_Base *buffer,int len,NA_Sequence *seq)
{
    int curlen=0,j;
    NA_Base *temp;
    temp=0;
    if(seq->seqlen+len >= seq->seqmaxlen)
    {
        if(seq->seqlen>0)
            seq->sequence = (NA_Base*)Realloc((char*)seq->sequence,
                                              (seq->seqlen + len+GBUFSIZ) * sizeof(NA_Base));
        else
            seq->sequence = (NA_Base*)Calloc(1,(seq->seqlen +
                                                len+GBUFSIZ) * sizeof(NA_Base));
        seq->seqmaxlen = seq->seqlen + len+GBUFSIZ;
    }
    /*
     *   seqlen is the length, and the index of the next free
     *   base
     */
    curlen = seq->seqlen + seq->offset;
    for(j=0;j<len;j++)
        putelem(seq,j+curlen,buffer[j]);

    seq->seqlen += len;
    return;
}

void Ascii2NA(char *buffer,int len,int matrix[16])
{
    /*
     *   if the translation matrix exists, use it to
     *   encode the buffer.
     */
    int i;
    if(matrix != NULL) {
        for(i=0;i<len;i++) {
            buffer[i] = matrix[(unsigned char)buffer[i]];
        }
    }
    return;
}

int WriteNA_Flat(NA_Alignment *aln,char *filename,int method,int maskable)
{
    size_t j;
    int kk,mask = -1,k,offset;
    char offset_str[100],buf[100];
    NA_Sequence *seqs;
    FILE *file;
    if(aln == NULL)
        return(1);
    if(aln->numelements == 0)
        return(1);
    seqs = aln->element;

    file = fopen(filename,"w");
    if(file == NULL)
    {
        Warning("Cannot open file for output");
        return(1);
    }
    if(maskable && (method != SELECT_REGION))
    {
        for(j=0;j<aln->numelements;j++)
            if(seqs[j].elementtype == MASK &&
               seqs[j].selected)
                mask = j;
    }
    /* Removed by OLIVER
       for(j=0;j<aln->numelements;j++)
       {
       SeqNorm(&(seqs[j]));
       }
    */

    for(j=0;j<aln->numelements;j++)
    {
        if (method != SELECT_REGION) {
            offset = seqs[j].offset;
        }
        else {
            for(offset=seqs[j].offset; aln->selection_mask[offset] == '0'; offset++) ;
        }

        if(offset+aln->rel_offset != 0)
            sprintf(offset_str,"(%d)",offset+aln->rel_offset);
        else
            offset_str[0] = '\0';

        if((((int)j!=mask) && (seqs[j].selected) && method != SELECT_REGION)
           || (method == SELECT_REGION && seqs[j].subselected)
           || method == ALL)
        {
            fprintf(file,"%c%s%s\n",
                    seqs[j].elementtype == DNA?'#':
                    seqs[j].elementtype == RNA?'#':
                    seqs[j].elementtype == PROTEIN?'%':
                    seqs[j].elementtype == TEXT?'"':
                    seqs[j].elementtype == MASK?'@':'"',
                    seqs[j].short_name,
                    (offset+aln->rel_offset  == 0)? "":offset_str);
            if(seqs[j].tmatrix)
            {
                if(mask == -1)
                    for(k=0,kk=0;kk<seqs[j].seqlen;kk++)
                    {
                        if((k)%60 == 0 && k>0)
                        {
                            buf[60] = '\0';
                            fputs(buf,file);
                            putc('\n',file);
                        }
                        if(method == SELECT_REGION)
                        {
                            if(aln->selection_mask[kk+offset]=='1')
                            {
                                buf[k%60] =((char)seqs[j].tmatrix[
                                                                  (int)getelem( &(seqs[j]),kk+offset) ]);
                                k++;
                            }
                        }
                        else
                        {
                            buf[k%60] =((char)seqs[j].tmatrix[
                                                              (int)getelem( &(seqs[j]),kk+offset) ]);
                            k++;
                        }
                    }
                else
                    for(k=0,kk=0;kk<seqs[j].seqlen;kk++)
                    {
                        if(getelem(&(seqs[mask]),kk+seqs[mask].offset) != '0'
                           && (getelem(&(seqs[mask]),kk+seqs[mask].offset)
                               != '-'))
                        {
                            if((k++)%60 == 0 && k>1)
                            {
                                buf[60] = '\0';
                                fputs(buf,file);
                                putc('\n',file);
                            }
                            buf[k%60] = ((char)seqs[j].tmatrix
                                         [getelem(&(seqs[j]),kk+offset)]);
                        }
                    }
            }
            else
            {
                if(mask == -1)
                    for(k=0,kk=0;kk<seqs[j].seqlen;kk++)
                    {
                        if((k)%60 == 0 && k>0)
                        {
                            buf[60] = '\0';
                            fputs(buf,file);
                            putc('\n',file);
                        }
                        if(method == SELECT_REGION)
                        {
                            if(aln->selection_mask[kk+offset]=='1')
                            {
                                buf[k%60] =(getelem( &(seqs[j]),kk+offset));
                                k++;
                            }
                        }
                        else
                        {
                            buf[k%60] =( getelem( &(seqs[j]),kk+offset) );
                            k++;
                        }
                    }
                else
                    for(k=0,kk=0;kk<seqs[j].seqlen;kk++)
                    {
                        if(getelem(&(seqs[mask]),kk+offset) == '1')
                        {
                            if((k++)%60 == 0 && k>1)
                            {
                                buf[60] = '\0';
                                fputs(buf,file);
                                putc('\n',file);
                            }
                            buf[k%60] =((char)getelem(&(seqs[j]),
                                                      kk+offset));
                        }
                    }
            }
            buf[(k%60)>0 ? (k%60):60] = '\0';
            fputs(buf,file);
            putc('\n',file);
        }
    }
    fclose(file);
    return(0);
}


void Warning(const char *s)
{
    /*extern Frame frame;
      extern Panel_item left_foot,right_foot;
      Beep();
      xv_set(frame,FRAME_RIGHT_FOOTER,s,0);
      xv_set(right_foot,PANEL_LABEL_STRING,s,0);*/
    aw_message(s);
}


void InitNASeq(NA_Sequence *seq,int type)
{


    SetTime(&(seq->t_stamp.origin));
    SetTime(&(seq->t_stamp.modify));
    strncpy(seq->id,uniqueID(),79);
    seq->seq_name[0] = '\0';
    seq->barcode[0] = '\0';
    seq->contig[0] = '\0';
    seq->membrane[0] = '\0';
    seq->authority[0] = '\0';
    seq->short_name[0] = '\0';
    seq->sequence = NULL;
    seq->offset = 0;
    seq->baggage = NULL;
    seq->baggage_len = 0;
    seq->baggage_maxlen = 0;
    seq->comments = NULL;
    seq->comments_len = 0;
    seq->comments_maxlen = 0;
    seq->description[0] = '\0';
    seq->mask = NULL;
    seq->seqlen = 0;
    seq->seqmaxlen = 0;
    seq->protect = PROT_WHITE_SPACE + PROT_TRANSLATION;
#ifdef HGL
    seq->attr = 0;
#else
    seq->attr = IS_5_TO_3 + IS_PRIMARY;
#endif
    seq->elementtype = type;
    seq->groupid = 0;
    seq->groupb = NULL;
    seq->groupf = NULL;
    seq->cmask = NULL;
    seq->selected = 0;
    seq->subselected = 0;

    switch (type)
    {
        case DNA:
            seq->tmatrix = Default_DNA_Trans;
            seq->rmatrix = Default_NA_RTrans;
            seq->col_lut = Default_NAColor_LKUP;
            break;
        case RNA:
            seq->tmatrix = Default_RNA_Trans;
            seq->rmatrix = Default_NA_RTrans;
            seq->col_lut = Default_NAColor_LKUP;
            break;
        case PROTEIN:
            seq->tmatrix = NULL;
            seq->rmatrix = NULL;
            seq->col_lut = Default_PROColor_LKUP;
            break;
        case MASK:
        case TEXT:
        default:
            seq->tmatrix = NULL;
            seq->rmatrix = NULL;
            seq->col_lut = NULL;
            break;
    }
    return;
}


void ReadCMask(const char *filename)
{

    char in_line[GBUFSIZ];
    char head[GBUFSIZ];
    char curname[GBUFSIZ];
    char temp[GBUFSIZ];
    bool IGNORE_DASH = false;
    int  offset;
    
    /*NA_DisplayData *NAdd;*/
    NA_Alignment *aln;

    size_t  j;
    size_t  curlen = 0;
    int    *colors = 0,orig_ctype,jj,indx = 0;
    FILE   *file;

    if(DataSet == NULL) return;

    /*NAdd = (NA_DisplayData*)((NA_Alignment*)DataSet)->na_ddata;

    if(NAdd == NULL)
    return;
    */
    aln = (NA_Alignment*)DataSet;

    curname[0] = '\0';
    orig_ctype = COLOR_MONO;
    file = fopen(filename,"r");
    if(file == NULL)
    {
        Warning("File not found");
        Warning(filename);
        return;
    }

    /*NAdd->color_type = COLOR_ALN_MASK;*/
    for(;fgets(in_line,GBUFSIZ,file) !=0;)
    {
        if(Find(in_line,"offset:"))
        {
            crop(in_line,head,temp);
            sscanf(temp,"%d",&(aln->cmask_offset));
        }
        else if(Find(in_line,"nodash:"))
            IGNORE_DASH = true;
        else if(Find(in_line,"dash:"))
            IGNORE_DASH = true;
        else if(Find(in_line,"name:"))
        {
            crop(in_line,head,curname);
            curname[strlen(curname)-1] = '\0';
            for(j=0;j<strlen(curname);j++)
                if(curname[j] == '(')
                    curname[j] = '\0';
        }
        else if(Find(in_line,"length:"))
        {
            crop(in_line,head,temp);
            sscanf(temp,"%zu",&curlen);
        }
        else if(Find(in_line,"start:"))
        {
            indx = -1;
            if(curlen == 0)
            {
                Warning("illegal format in colormask");
                /*NAdd->color_type = orig_ctype;*/
                return;
            }
            if(strlen(curname) != 0)
            {
                indx = -1;
                for(j=0;j<aln->numelements;j++)
                    if(Find(aln->element[j].short_name,curname)
                       || Find(aln->element[j].id,curname))
                    {
                        if(aln->element[j].cmask != NULL)
                            Cfree((char*)aln -> element[j].cmask);
                        colors=(int*)Calloc(aln->element[j]
                                            .seqmaxlen+1+aln->element[j].offset
                                            ,sizeof(int));
                        aln->element[j].cmask = colors;
                        /*NAdd->color_type = COLOR_SEQ_MASK;*/
                        indx = j;
                        j = aln->numelements;
                    }
                if(indx == -1)
                    colors=NULL;
            }
            else
            {
                if(aln->cmask != NULL) Cfree((char*)aln->cmask);
                colors=(int*)Calloc(curlen,sizeof(int));
                aln->cmask = colors;
                aln->cmask_len = curlen;
                /*NAdd->color_type = COLOR_ALN_MASK;*/
                for(j=0;j<curlen;j++)
                    colors[j] = 12;
            }

            if(IGNORE_DASH && (indx != -1))
            {
                for(jj=0,j=0;(j<curlen) &&
                        (jj<aln->element[indx].seqlen);j++,jj++)
                {
                    offset = aln->element[indx].offset;
                    if(fgets(in_line,GBUFSIZ,file)==NULL)
                    {
                        Warning
                            ("illegal format in colormask");
                        /*NAdd->color_type = orig_ctype;*/
                        return;
                    }
                    /*
                     *  Fixed so that the keyword nodash causes the colormask to be mapped
                     *  to the sequence, not the alignment.
                     *
                     *  The allocated space is equal the seqlen of the matched sequence.
                     *
                     */
                    if(aln->element[indx].tmatrix)
                        for(;(getelem(&(aln->element[indx]),jj
                                      +offset)
                              ==(aln->element[indx].tmatrix['-'])
                              || (getelem(&(aln->element[indx]),jj
                                          +offset)
                                  ==aln->element[indx].tmatrix['~']))
                                && jj < aln->element[indx].seqlen;)
                            colors[jj++] = 12;
                    else
                        for(;getelem(&(aln->element[indx]),jj
                                     +offset)
                                =='-' && jj < aln->element[indx].seqlen;)
                            colors[jj++] = 12;

                    sscanf(in_line,"%d",&(colors[jj]));
                }
            }
            else if((indx == -1)  && (strlen(curname) != 0))
                for(j=0;j<curlen;j++)
                    fgets(in_line,GBUFSIZ,file);
            else
                for(j=0;j<curlen;j++)
                {
                    if(fgets(in_line,GBUFSIZ,file)==NULL)
                    {
                        Warning
                            ("illegal format in colormask");
                        /*NAdd->color_type = orig_ctype;*/
                        return;
                    }
                    sscanf(in_line,"%d",&(colors[j]));
                }
            IGNORE_DASH = false;
            curname[0] = '\0';
        }

    }
    /*RepaintAll(true);*/
    return;
}


int WriteStatus(NA_Alignment *aln,char *filename,int method)
{
    //  extern int EditMode;
    //  NA_DisplayData *NAdd;
    NA_Sequence *this_seq;
    int j;
    FILE *file;
    method=0;filename=0;

    if(DataSet == NULL)
        return(1);

    /*
      NAdd = (NA_DisplayData*)((NA_Alignment*)DataSet)->na_ddata;
      if(NAdd == NULL)
      return(1);
    */

    file = fopen(filename,"w");
    if (file == NULL)
    {
        Warning("Cannot open status file.");
        return(1);
    }
    fprintf(file,"File_format: %s\n",FileFormat==GENBANK?"genbank":"flat");
    /*
      fprintf(file,"EditMode: %s\n",EditMode==INSERT?"insert":
      "check");
    */

    this_seq = &(aln->element[1]); /* Nadd->cursor !? */
    if(this_seq->id != NULL)
        fprintf(file,"sequence-ID %s\n",this_seq->id);
    fprintf(file,"Column: %d\nPos:%d\n",1,1);/*NAdd->cursor_x,NAdd->position*/
    switch(this_seq->elementtype)
    {
        case DNA:
        case RNA:
            fprintf(file,"#%s\n",
                    this_seq->short_name);
            break;
        case PROTEIN:
            fprintf(file,"%%%s\n",
                    this_seq->short_name);
            break;
        case MASK:
            fprintf(file,"@%s\n",
                    this_seq->short_name);
            break;
        case TEXT:
            fprintf(file,"%c%s\n",'"',
                    this_seq->short_name);
            break;
        default:
            break;
    }
    if(this_seq->tmatrix)
        for(j=0;j<this_seq->seqlen;j++)
            putc(this_seq->tmatrix[getelem(this_seq,j)],file);
    else
        for(j=0;j<this_seq->seqlen;j++)
            putc(getelem(this_seq,j),file);

    fclose(file);
    return(0);
}

void ReadStatus(char *filename)
{
    filename=0;
    /*
      int i,j;
      FILE *file;
      filename=0;

      char in_line[GBUFSIZ],head[GBUFSIZ];
      file = fopen(filename,"r");
      for(;!DONE;)
      {
      fgets(in_line,GBUFSIZ,file);
      if(strlen(in_line) == 0)
      DONE = true;
      else
      {
      sscanf(in_line,"%s",head);
      if(strncmp(head,"Col",3) != 0)
      {
      sscanf(in_line,"%*s %d",head,&(DataSet->nadd->
      cursor_x),&(DataSet->nadd->cursory);
      }
      else if(strncmp(head,"Pos",3) != 0)
      {
      }
      }
      }

    */
}


void NormalizeOffset(NA_Alignment *aln)
{
    int i;
    size_t j;
    int offset = 99999999;
    i=0;

    for(j=0;j<aln->numelements;j++)
        offset = MIN(offset,aln->element[j].offset);

    for(j=0;j<aln->numelements;j++)
        aln->element[j].offset -= offset;

    aln->maxlen = -999999999;
    for(j=0;j<aln->numelements;j++)
        aln->maxlen = MAX(aln->element[j].seqlen+aln->element[j].offset,
                          aln->maxlen);

    aln->rel_offset += offset;

    if(aln->numelements == 0)
        aln->rel_offset = 0;

    return;
}

int WriteCMask(NA_Alignment *aln,char *filename,int method,int maskable)
{
    size_t j;
    int kk,mask = -1,k,offset;
    char offset_str[100];
    int *buf;
    NA_Sequence *seqs;
    FILE *file;
    if(aln == NULL)
        return(1);

    if(aln->numelements == 0)
        return(1);
    seqs = aln->element;

    file = fopen(filename,"w");
    if(file == NULL)
    {
        Warning("Cannot open file for output");
        return(1);
    }
    if(maskable && (method != SELECT_REGION))
    {
        for(j=0;j<aln->numelements;j++)
            if(seqs[j].elementtype == MASK &&
               seqs[j].selected)
                mask = j;
    }
    for(j=0;j<aln->numelements;j++)
    {
        SeqNorm(&(seqs[j]));
    }

    for(j=0;j<aln->numelements;j++)
    {
        if(method != SELECT_REGION) {
            offset = seqs[j].offset;
        }
        else {
            for(offset=seqs[j].offset; aln->selection_mask[offset] == '0'; offset++) ;
        }

        if(offset+aln->rel_offset != 0) {
            sprintf(offset_str,"(%d)",offset+aln->rel_offset);
        }
        else {
            offset_str[0] = '\0';
        }

        if((((int)j!=mask) && (seqs[j].selected) && method != SELECT_REGION)
           || (method == SELECT_REGION && seqs[j].subselected)
           || method == ALL)
        {
            fprintf(file,"%c%s%s\n",
                    seqs[j].elementtype == DNA?'#':
                    seqs[j].elementtype == RNA?'#':
                    seqs[j].elementtype == PROTEIN?'%':
                    seqs[j].elementtype == TEXT?'"':
                    seqs[j].elementtype == MASK?'@':'"',
                    seqs[j].short_name,
                    (offset+aln->rel_offset  == 0)? "":offset_str);

            if(seqs[j].cmask != NULL)
            {

                buf =(int*) Calloc(seqs[j].seqlen,sizeof(int) );

                if(mask == -1)
                {
                    for(k=0,kk=0;kk<seqs[j].seqlen;kk++)
                    {
                        if(method == SELECT_REGION)
                        {
                            if(aln->selection_mask[kk+offset]=='1')
                                buf[k++] = (getcmask( &(seqs[j]),kk+offset));
                        }

                        else
                            buf[k++] =( getcmask( &(seqs[j]),kk+offset) );
                    }
                }
                else
                {
                    for(k=0,kk=0;kk<seqs[j].seqlen;kk++)
                        if(getelem(&(seqs[mask]),kk+offset) == '1')
                            buf[k++] =(getcmask(&(seqs[j]), kk+offset));
                    /*
                     *          Looks like k might be one behind?
                     */
                }
                fprintf(file,"name:%s\noffset:%d\nlength:%d\nstart:\n",
                        seqs[j].short_name,seqs[j].offset,k);

                for(kk = 0; kk < k; kk++)
                    fprintf(file, "%d\n", buf[kk]);

                Cfree((char*)buf);
            }
        }
    }
    fclose(file);
    return(0);
}