--- /dev/null
+/* -*- mode: c; tab-width: 4; c-basic-offset: 4; indent-tabs-mode: nil -*- */
+
+/*********************************************************************
+ * Clustal Omega - Multiple sequence alignment
+ *
+ * Copyright (C) 2010 University College Dublin
+ *
+ * Clustal-Omega is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU General Public License as
+ * published by the Free Software Foundation; either version 2 of the
+ * License, or (at your option) any later version.
+ *
+ * This file is part of Clustal-Omega.
+ *
+ ********************************************************************/
+
+/*
+ * RCS $Id: hhfullalignment-C.h 243 2011-05-31 13:49:19Z fabian $
+ */
+
+// hhfullalignment.C
+
+#ifndef MAIN
+#define MAIN
+#include <iostream> // cin, cout, cerr
+#include <fstream> // ofstream, ifstream
+#include <stdio.h> // printf
+#include <stdlib.h> // exit
+#include <string> // strcmp, strstr
+#include <math.h> // sqrt, pow
+#include <limits.h> // INT_MIN
+#include <float.h> // FLT_MIN
+#include <time.h> // clock
+#include <ctype.h> // islower, isdigit etc
+using std::ios;
+using std::ifstream;
+using std::ofstream;
+using std::cout;
+using std::cerr;
+using std::endl;
+#include "util-C.h" // imax, fmax, iround, iceil, ifloor, strint, strscn, strcut, substr, uprstr, uprchr, Basename etc.
+#include "list.h" // list data structure
+#include "hash.h" // hash data structure
+#include "hhdecl-C.h" // constants, class
+#include "hhutil-C.h" // imax, fmax, iround, iceil, ifloor, strint, strscn, strcut, substr, uprstr, uprchr, Basename etc.
+#include "hhhmm.h" // class HMM
+#include "hhalignment.h" // class Alignment
+#include "hhhit.h"
+#include "hhhalfalignment.h"
+#endif
+
+
+
+//////////////////////////////////////////////////////////////////////////////
+//////////////////////////////////////////////////////////////////////////////
+// Methods of class FullAlignment
+//////////////////////////////////////////////////////////////////////////////
+//////////////////////////////////////////////////////////////////////////////
+
+
+//////////////////////////////////////////////////////////////////////////////
+// Output Results: use classes HalfAlignment and FullAlignment
+//
+// Example:
+// Each column list contains at least a match state
+// Insert states between the match states are omitted
+//
+// step 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1
+// i 0 0 1 2 3 4 5 6 7 8 9 9 9 9 10 11 12 13 14
+//
+// Q ~ ~ X X X X X X X X X ~ ~ ~ X X X X X
+// T Y Y Y Y Y Y Y Y ~ Y Y Y Y Y Y Y Y ~ ~
+//
+// j 7 8 9 10 11 12 13 14 14 15 16 17 18 19 20 21 22 22 22
+// state IM IM MM MM MM MM MM MM DG MM MM GD GD GD MM MM MM MI MI
+//
+// nsteps=19
+//
+
+
+//////////////////////////////////////////////////////////////////////////////
+// Constructor
+FullAlignment::FullAlignment(int maxseqdis)
+{
+ qa = new HalfAlignment(maxseqdis);
+ ta = new HalfAlignment(maxseqdis);
+ for (int h=0; h<LINELEN-1; h++) symbol[h]=' ';
+}
+
+//////////////////////////////////////////////////////////////////////////////
+// Destructor
+FullAlignment::~FullAlignment()
+{
+ delete qa; qa = NULL;
+ delete ta; ta = NULL;
+}
+
+//////////////////////////////////////////////////////////////////////////////
+// Free memory of arrays s[][], l[][], and m[][] in HalfAlignment
+void FullAlignment::FreeMemory()
+{
+ qa->Unset();
+ ta->Unset();
+}
+
+//////////////////////////////////////////////////////////////////////////////
+/**
+ * @brief Add columns for match (and delete) states.
+ */
+void
+FullAlignment::AddColumns(int i, int j, char prev_state, char state, float S)
+{
+ switch(state)
+ {
+ case MM: //MM pair state (both query and template in Match state)
+ AddGaps(); //fill up gaps until query and template parts have same length
+ symbol[qa->pos] =ScoreChr(S);
+ qa->AddColumn(i);
+ ta->AddColumn(j);
+ qa->AddInsertsAndFillUpGaps(i);
+ ta->AddInsertsAndFillUpGaps(j);
+ break;
+
+ case GD: //-D state
+ if (prev_state==DG) AddGaps();
+ symbol[ta->pos]='Q';
+ ta->AddColumn(j); //query has gap -> add nothing
+ ta->AddInsertsAndFillUpGaps(j);
+ break;
+ case IM: //IM state
+ if (prev_state==MI) AddGaps();
+ symbol[ta->pos]='Q';
+ ta->AddColumn(j); //query has gap -> add nothing
+ ta->AddInsertsAndFillUpGaps(j);
+ break;
+
+ case DG: //D- state
+ if (prev_state==GD) AddGaps();
+ symbol[qa->pos]='T';
+ qa->AddColumn(i);//template has gap -> add nothing
+ qa->AddInsertsAndFillUpGaps(i);
+ break;
+ case MI: //MI state
+ if (prev_state==IM) AddGaps();
+ symbol[qa->pos]='T';
+ qa->AddColumn(i);//template has gap -> add nothing
+ qa->AddInsertsAndFillUpGaps(i);
+ break;
+ }
+}
+
+/////////////////////////////////////////////////////////////////////////////////////
+/**
+ * @brief Fill up gaps until query and template parts have same length
+ */
+void
+FullAlignment::AddGaps()
+{
+ while (qa->pos<ta->pos) qa->AddChar('.');
+ while (ta->pos<qa->pos) ta->AddChar('.');
+}
+
+
+//////////////////////////////////////////////////////////////////////////////
+/**
+ * @brief Build full alignment -> qa->s[k][h] and ta->s[k][h]
+ */
+int
+FullAlignment::Build(HMM& q, Hit& hit)
+{
+ int step;
+ char prev_state=MM, state=MM;
+ int n;
+ int hh;
+ int k;
+ identities=0; // number of identical residues in query and template sequence
+ score_sim=0.0f;// substitution matrix similarity score between query and template
+
+ ClearSymbols();
+
+ // Set up half-alignments
+ // n is the sequence index up to which sequences are prepared for display
+ n = imin( q.n_display,par.nseqdis+( q.nss_dssp>=0)+( q.nsa_dssp>=0)+( q.nss_pred>=0)+( q.nss_conf>=0)+( q.ncons>=0));
+ qa->Set( q.name, q.seq, q.sname, n, q.L, q.nss_dssp , q.nss_pred , q.nss_conf, q.nsa_dssp, q.ncons, hit.L/*<--FS*/);
+ n = imin(hit.n_display,par.nseqdis+(hit.nss_dssp>=0)+(hit.nsa_dssp>=0)+(hit.nss_pred>=0)+(hit.nss_conf>=0)+(hit.ncons>=0));
+ ta->Set(hit.name,hit.seq,hit.sname, n,hit.L,hit.nss_dssp,hit.nss_pred,hit.nss_conf,hit.nsa_dssp, hit.ncons, q.L/*<--FS*/);
+
+
+// printf("HMM: %s\nstep nst i j state hq ht\n",hit.name);
+
+#ifdef CLUSTALO
+ if ((1 == hit.i1) && (1 == hit.j1)){
+ /* do nothing */
+ }
+ else if ((1 == hit.i1) && (1 != hit.j1)){
+ for (int j = 1; j < hit.j1; j++){
+ AddColumns(0, j, MM, IM, 0.0);
+ }
+ if (rLog.iLogLevelEnabled <= LOG_DEBUG){
+ int iK, iL;
+ printf("%d: j1=%d -> temp has leading gaps\n", __LINE__, hit.j1);
+ for (iK = 0; iK < ta->n; iK++){
+ //printf("T%d: %s\n", iK, ta->s[iK]);
+ for (iL = 0; iL < hit.j1; iL++){
+ ta->s[iK][iL] = tolower(ta->s[iK][iL]);
+ }
+ }
+ }
+ }
+ else if ((1 != hit.i1) && (1 == hit.j1)){
+ for (int i = 1; i < hit.i1; i++){
+ AddColumns(i, 0, MM, MI, 0.0);
+ }
+ if (rLog.iLogLevelEnabled <= LOG_DEBUG){
+ FILE *fp = rLog.prFP[LOG_DEBUG];
+
+ fprintf(fp, "%d: i1=%d -> query has leading gaps\n", __LINE__, hit.i1);
+ int iK, iL;
+ for (iK = 0; iK < qa->n; iK++){
+ //printf("Q%d: %s\n", iK, qa->s[iK]);
+ for (iL = 0; iL < hit.i1; iL++){
+ qa->s[iK][iL] = tolower(qa->s[iK][iL]);
+ }
+ }
+ }
+ }
+ else {
+ fprintf(stderr,
+ "+-------------------------------+\n"
+ "| both sequences truncated left |\n"
+ "+-------------------------------+\n"
+ "i1 = %d, j1 = %d\n", hit.i1, hit.j1);
+ return FAILURE; /* FS, r241 -> r243 */
+ }
+#endif
+
+ for (step=hit.nsteps; step>=1; step--)
+ {
+ prev_state = state;
+ state = hit.states[step];
+
+ // Add column to alignment and compute identities and sequence-sequence similarity score
+ AddColumns(hit.i[step],hit.j[step],prev_state,state,hit.S[step]);
+ if (state==MM)
+ {
+ char qc=qa->seq[ q.nfirst][ qa->m[ q.nfirst][hit.i[step]] ];
+ char tc=ta->seq[hit.nfirst][ ta->m[hit.nfirst][hit.j[step]] ];
+ if (qc==tc) identities++; // count identical amino acids
+ score_sim += S[(int)aa2i(qc)][(int)aa2i(tc)];
+ // fprintf(stderr,"%3i %3i %3i %3i %3i %1c %1c %6.2f %6.2f %6.2f %6.2f \n",step,hit.nsteps,hit.i[step],hit.j[step],int(state),qc,tc,S[(int)aa2i(qc)][(int)aa2i(tc)],score_sim,hit.P_posterior[step],hit.sum_of_probs); //DEBUG
+ }
+ }
+
+#ifdef CLUSTALO
+ if ((qa->L == hit.i2) && (ta->L == hit.j2)){
+ /* do nothing */
+ }
+ else if ((qa->L == hit.i2) && (ta->L != hit.j2)){
+ for (int j = hit.j2+1; j <= ta->L; j++){
+ AddColumns(0, j, MM, IM, 0.0);
+ }
+ if (rLog.iLogLevelEnabled <= LOG_DEBUG){
+ int iK;
+ unsigned int uiL;
+ FILE *fp = rLog.prFP[LOG_DEBUG];
+
+ fprintf(fp, "%d: j2=%d (%d) -> temp has trailing gaps\n", __LINE__, hit.j2, ta->L);
+ for (iK = 0; iK < ta->n; iK++){
+ //printf("T%d: %s\n", iK, ta->s[iK]+strlen(ta->s[iK])-(ta->L-hit.j2));
+ for (uiL = strlen(ta->s[iK])-(ta->L-hit.j2); uiL < strlen(ta->s[iK]); uiL++){
+ ta->s[iK][uiL] = tolower(ta->s[iK][uiL]);
+ }
+ }
+ }
+ }
+ else if ((qa->L != hit.i2) && (ta->L == hit.j2)){
+ for (int i = hit.i2+1; i <= qa->L; i++){
+ AddColumns(i, 0, MM, MI, 0.0);
+ }
+ if (rLog.iLogLevelEnabled <= LOG_DEBUG){
+ int iK;
+ unsigned int uiL;
+ printf("%d: i2=%d (%d)-> query has trailing gaps\n", __LINE__, hit.i2, qa->L);
+ for (iK = 0; iK < qa->n; iK++){
+ //printf("Q%d: %s\n", iK, qa->s[iK]+strlen(qa->s[iK])-(qa->L-hit.i2));
+ for (uiL = strlen(qa->s[iK])-(qa->L-hit.i2); uiL < strlen(qa->s[iK]); uiL++){
+ qa->s[iK][uiL] = tolower(qa->s[iK][uiL]);
+ }
+ }
+ }
+ }
+ else{
+ fprintf(stderr,
+ "+--------------------------------+\n"
+ "| both sequences truncated right |\n"
+ "+--------------------------------+\n"
+ "i2 = %d != %d = qa->L, j2 = %d != %d = ta->L\n",
+ hit.i2, qa->L, hit.j2, ta->L);
+ return FAILURE; /* FS, r241 -> r243 */
+ }
+
+#endif
+
+ AddGaps(); //fill up gaps until query and template parts have same length
+ qa->AddChar('\0');
+ ta->AddChar('\0');
+
+ // Change gap symbol '.' (gap aligned to insert) to '~' if one HMM has gap with respect to other HMM
+ for (hh=1; hh<qa->pos; hh++)
+ {
+ if (symbol[hh]=='Q')
+ {
+ // Gap in query (IM or GD state)
+ symbol[hh]=' ';
+ for (k=0; k<qa->n; k++) if (qa->s[k][hh]=='.') qa->s[k][hh]='-';
+ }
+ else if (symbol[hh]=='T')
+ {
+ // Gap in target (MI or DG state)
+ symbol[hh]=' ';
+ for (k=0; k<ta->n; k++) if (ta->s[k][hh]=='.') ta->s[k][hh]='-';
+ }
+ }
+ return OK;
+} /* this is the end of FullAlignment::Build() */
+
+
+//////////////////////////////////////////////////////////////////////////////
+/**
+ * @brief Print out header before full alignment
+ */
+void
+FullAlignment::PrintHeader(FILE* outf, HMM& q, Hit& hit)
+{
+ fprintf(outf,">%s\n",hit.longname);
+ fprintf(outf,"Probab=%-.2f E-value=%-.2g Score=%-.2f Aligned_cols=%i Identities=%i%% Similarity=%-.3f Sum_probs=%.1f\n\n",
+ hit.Probab,hit.Eval,hit.score,hit.matched_cols,iround(100.0*identities/hit.matched_cols),score_sim/hit.matched_cols,hit.sum_of_probs);
+}
+
+//////////////////////////////////////////////////////////////////////////////
+/**
+ * @brief Print out full alignment in HHR format
+ */
+void
+FullAlignment::PrintHHR(FILE* outf, Hit& hit)
+{
+ const int NLEN=14; //Length of name field in front of multiple alignment
+ int h=0; //counts position (column) in alignment
+ int hh=0; //points to column at start of present output block of alignment
+ int k; //counts sequences in query and template
+ //short unsigned int lq[MAXSEQ]; // lq[k] counts index of residue from query sequence k to be printed next;
+ short unsigned int *lq = NULL; // lq[k] counts index of residue from query sequence k to be printed next;
+ //short unsigned int lt[MAXSEQ]; // lt[k] counts index of residue from template sequence k to be printed next;
+ short unsigned int *lt = NULL; // lt[k] counts index of residue from template sequence k to be printed next;
+ char namestr[NAMELEN]; //name of sequence
+ int iq=hit.i1; // match state counter for query HMM (displayed in consensus line)
+ int jt=hit.j1; // match state counter for template HMM (displayed in consensus line)
+
+ lq = new(short unsigned int[qa->n+2]);
+ lt = new(short unsigned int[ta->n+2]);
+
+ for (k=0; k<qa->n; k++) lq[k]=qa->l[k][hit.i1];
+ for (k=0; k<ta->n; k++) lt[k]=ta->l[k][hit.j1];
+
+ while (hh<ta->pos-1) // print alignment block
+ {
+ // Print query secondary structure sequences
+ for (k=0; k<qa->n; k++)
+ {
+ if (k==qa->nsa_dssp) continue;
+ if (!(k==qa->nss_dssp || k==qa->nsa_dssp || k==qa->nss_pred || k==qa->nss_conf)) continue;
+ if (k==qa->nss_dssp && !par.showdssp) continue;
+ if ((k==qa->nss_pred || k==qa->nss_conf) && !par.showpred) continue;
+ strncpy(namestr,qa->sname[k],NAMELEN-2);
+ namestr[NAMELEN-1]='\0';
+ strcut(namestr);
+ fprintf(outf,"Q %-*.*s ",NLEN,NLEN,namestr);
+ for (h=hh; h<imin(hh+par.aliwidth,qa->pos-1); h++) fprintf(outf,"%1c",qa->s[k][h]);
+ fprintf(outf,"\n");
+ }
+
+ // Print query sequences
+ for (k=0; k<qa->n; k++)
+ {
+ if (k==qa->nss_dssp || k==qa->nsa_dssp || k==qa->nss_pred || k==qa->nss_conf || k==qa->ncons) continue;
+ strncpy(namestr,qa->sname[k],NAMELEN-2);
+ namestr[NAMELEN-1]='\0';
+ strcut(namestr);
+ fprintf(outf,"Q %-*.*s %4i ",NLEN,NLEN,namestr,lq[k]);
+ for (h=hh; h<imin(hh+par.aliwidth,qa->pos-1); h++)
+ {fprintf(outf,"%1c",qa->s[k][h]); lq[k]+=WordChr(qa->s[k][h]);} //WordChr() returns 1 if a-z or A-Z; 0 otherwise
+ fprintf(outf," %4i (%i)\n",lq[k]-1,qa->l[k][qa->L+1]);
+ }
+
+ // Print query consensus sequence
+ if (par.showcons && qa->ncons>=0)
+ {
+ k=qa->ncons;
+ strncpy(namestr,qa->sname[k],NAMELEN-2);
+ namestr[NAMELEN-1]='\0';
+ strcut(namestr);
+ fprintf(outf,"Q %-*.*s %4i ",NLEN,NLEN,namestr,iq);
+ for (h=hh; h<imin(hh+par.aliwidth,qa->pos-1); h++)
+ {
+ if (qa->s[k][h]=='x') qa->s[k][h]='~';
+ if (qa->s[k][h]!='-' && qa->s[k][h]!='.') iq++;
+ fprintf(outf,"%1c",qa->s[k][h]);
+ }
+ fprintf(outf," %4i (%i)\n",iq-1,qa->L);
+ }
+
+
+ // Print symbols representing the score
+ fprintf(outf," %*.*s ",NLEN,NLEN," ");
+ for (h=hh; h<imin(hh+par.aliwidth,qa->pos-1); h++) fprintf(outf,"%1c",symbol[h]);
+ fprintf(outf,"\n");
+
+ // Print template consensus sequence
+ if (par.showcons && ta->ncons>=0)
+ {
+ k=ta->ncons;
+ strncpy(namestr,ta->sname[k],NAMELEN-2);
+ namestr[NAMELEN-1]='\0';
+ strcut(namestr);
+ fprintf(outf,"T %-*.*s %4i ",NLEN,NLEN,namestr,jt);
+ for (h=hh; h<imin(hh+par.aliwidth,ta->pos-1); h++)
+ {
+ if (ta->s[k][h]=='x') ta->s[k][h]='~';
+ if (ta->s[k][h]!='-' && ta->s[k][h]!='.') jt++;
+ fprintf(outf,"%1c",ta->s[k][h]);
+ }
+ fprintf(outf," %4i (%i)\n",jt-1,ta->L);
+ }
+ // Print template sequences
+ for (k=0; k<ta->n; k++)
+ {
+ if (k==ta->nss_dssp || k==ta->nsa_dssp || k==ta->nss_pred || k==ta->nss_conf || k==ta->ncons) continue;
+ strncpy(namestr,ta->sname[k],NAMELEN-2);
+ namestr[NAMELEN-1]='\0';
+ strcut(namestr);
+ fprintf(outf,"T %-*.*s %4i ",NLEN,NLEN,namestr,lt[k]);
+ for (h=hh; h<imin(hh+par.aliwidth,ta->pos-1); h++)
+ {fprintf(outf,"%1c",ta->s[k][h]); lt[k]+=WordChr(ta->s[k][h]);} //WordChr() returns 1 if a-z or A-Z; 0 otherwise
+ fprintf(outf," %4i (%i)\n",lt[k]-1,ta->l[k][ta->L+1]);
+ }
+
+ // Print template secondary structure sequences
+ for (k=0; k<ta->n; k++)
+ {
+ if (k==ta->nsa_dssp) continue;
+ if (!(k==ta->nss_dssp || k==ta->nss_pred || k==ta->nss_conf)) continue;
+ if (k==ta->nss_dssp && !par.showdssp) continue;
+ if ((k==ta->nss_pred || k==ta->nss_conf)&& !par.showpred) continue;
+ strncpy(namestr,ta->sname[k],NAMELEN-2);
+ namestr[NAMELEN-1]='\0';
+ strcut(namestr);
+ fprintf(outf,"T %-*.*s ",NLEN,NLEN,namestr);
+ for (h=hh; h<imin(hh+par.aliwidth,ta->pos-1); h++) fprintf(outf,"%1c",ta->s[k][h]);
+ fprintf(outf,"\n");
+ }
+
+ hh=h;
+ fprintf(outf,"\n\n");
+ }
+
+ delete[](lq); lq = NULL;
+ delete[](lt); lt = NULL;
+
+} /* this is the rnd of PrintHHR() */
+
+//////////////////////////////////////////////////////////////////////////////
+// Print out full alignment in A2M format
+//////////////////////////////////////////////////////////////////////////////
+#define TELOMER_PRINT 0
+void FullAlignment::PrintA2M(FILE* outf, Hit& hit)
+{
+ int k; //counts sequences in query and template
+ int h,hh;
+#if TELOMER_PRINT
+ int iTelomerLeft = hit.i1 - hit.j1;
+ int iTelomerRght = (qa->L - hit.i2) - (ta->L - hit.j2);
+#endif
+
+ // Print query sequences
+ for (k=0; k<qa->n; k++)
+ {
+ if (k==qa->nsa_dssp) continue;
+ if (k==qa->nss_dssp && !par.showdssp) continue;
+ if ((k==qa->nss_pred || k==qa->nss_conf) && !par.showpred) continue;
+ if (k==qa->ncons && !par.showcons) continue;
+ fprintf(outf,">%s\n",qa->sname[k]);
+#if TELOMER_PRINT
+ /* @<print left cut-off@> */
+ if (iTelomerLeft > 0){
+ for (int iI = 1; iI <= iTelomerLeft; iI++){
+ fprintf(outf, "%1c", qa->seq[k][iI]);
+ }
+ /* this may be unnecessary */
+ for (int iI = iTelomerLeft+1; iI < hit.i1; iI++){
+ fprintf(outf, "%1c", qa->seq[k][iI]);
+ }
+ }
+ else {
+ for (int iI = iTelomerLeft; iI < 0; iI++){
+ fprintf(outf, "%1c", '-');
+ }
+ /* this may be unnecessary */
+ /* think about it */
+ }
+ /* @<print consensus@> */
+ for (h = 0; qa->s[k][h] > 0; h++){
+ fprintf(outf, "%1c", qa->s[k][h]);
+ }
+ /* @<print out rght cut-off@> */
+ for (int iI = hit.i2+1; iI <= qa->L; iI++){
+ fprintf(outf, "%1c", qa->seq[k][iI]);
+ }
+ for (int iI = 0; iI > iTelomerRght; iI--){
+ fprintf(outf, "%1c", '-');
+ }
+#else
+ for (h=0,hh=-par.aliwidth; qa->s[k][h]>0; h++,hh++)
+ {
+ if (!hh) {fprintf(outf,"\n"); hh-=par.aliwidth;}
+ fprintf(outf,"%1c",qa->s[k][h]);
+ }
+#endif
+ fprintf(outf,"\n");
+ }
+
+ // Print template sequences
+ for (k=0; k<ta->n; k++)
+ {
+ if (k==ta->nsa_dssp) continue;
+ if (k==ta->nss_dssp && !par.showdssp) continue;
+ if ((k==ta->nss_pred || k==ta->nss_conf) && !par.showpred) continue;
+ if (k==ta->ncons && !par.showcons) continue;
+ fprintf(outf,">%s\n",ta->sname[k]);
+#if TELOMER_PRINT
+ /* @<print left cut-off@> */
+ if (iTelomerLeft > 0){
+ for (int iI = 1; iI <= iTelomerLeft; iI++){
+ fprintf(outf, "%1c", '-');
+ }
+ /* this may be unnecessary */
+ for (int iI = iTelomerLeft+1; iI < hit.j1; iI++){
+ fprintf(outf, "%1c", ta->seq[k][iI]);
+ }
+ }
+ else{
+ for (int iI = 1; iI <= -iTelomerLeft; iI++){
+ fprintf(outf, "%1c", ta->seq[k][iI]);
+ }
+ /* this may be unnecessary */
+ /* think about it */
+ }
+ /* @<print consensus@> */
+ for (h = 0; ta->s[k][h] > 0; h++){
+ fprintf(outf, "%1c", ta->s[k][h]);
+ }
+
+ /* @<print out rght cut-off@> */
+ for (int iI = hit.j2+1; iI <= ta->L; iI++){
+ fprintf(outf, "%1c", ta->seq[k][iI]);
+ }
+ for (int iI = 0; iI < iTelomerRght; iI++){
+ fprintf(outf, "%1c", '-');
+ }
+
+#else
+ for (h=0,hh=-par.aliwidth; ta->s[k][h]>0; h++,hh++)
+ {
+ if (!hh) {fprintf(outf,"\n"); hh-=par.aliwidth;}
+ fprintf(outf,"%1c",ta->s[k][h]);
+ }
+#endif
+ fprintf(outf,"\n");
+ }
+ fprintf(outf,"\n");
+
+} /* this is the end of PrintA2M() */
+
+////////////////////////////////////////////////////////////////////////////
+/**
+ * @brief Print out full alignment in A2M format
+ */
+void
+FullAlignment::PrintFASTA(FILE* outf, Hit& hit)
+{
+ // Transform sequences to uppercase and '.' to '-'
+ qa->ToFASTA();
+ ta->ToFASTA();
+ PrintA2M(outf,hit);
+}
+
+//////////////////////////////////////////////////////////////////////////////
+/**
+ * @brief Print out full alignment in A2M format
+ */
+void
+FullAlignment::PrintA3M(FILE* outf, Hit& hit)
+{
+ // Remove all '.' from sequences
+ qa->RemoveChars('.');
+ ta->RemoveChars('.');
+ PrintA2M(outf,hit);
+}
+
+
+/*
+ * @* OverWriteSeqs()
+ */
+void
+FullAlignment::OverWriteSeqs(char **ppcFirstProf, char **ppcSecndProf){
+
+ int iS, iR;
+ char cRes;
+
+ for (iS = 0; iS < qa->n; iS++){
+ for (iR = 0; iR < qa->pos; iR++){
+ cRes = qa->s[iS][iR];
+ ppcFirstProf[iS][iR] = ('.' == cRes) ? '-' : cRes;
+
+ } /* 0 <= iR < qa.L */
+ ppcFirstProf[iS][iR] = '\0';
+ } /* 0 <= iS < qa.n */
+
+ for (iS = 0; iS < ta->n; iS++){
+ for (iR = 0; iR < ta->pos; iR++){
+ cRes = ta->s[iS][iR];
+ ppcSecndProf[iS][iR] = ('.' == cRes) ? '-' : cRes;
+
+ } /* 0 <= iR < ta.L */
+ ppcSecndProf[iS][iR] = '\0';
+ } /* 0 <= iS < ta.n */
+
+} /* this is the end of OverWriteSeqs() */
+
+