X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=sequence.hpp;h=6a50cb0b8aa31e4a3a3b6c9c97e77416f53eeeb1;hb=8dd3c225255d7084e3aff8740aa4f1f1cabb367a;hp=03cbab7ffddce5531e272cf6116069edfa0ebade;hpb=c196b6b4768ccb84955d773ff0f22e4994d1ba7b;p=mothur.git

diff --git a/sequence.hpp b/sequence.hpp
index 03cbab7..6a50cb0 100644
--- a/sequence.hpp
+++ b/sequence.hpp
@@ -8,11 +8,17 @@
  *  Created by Pat Schloss on 12/15/08.
  *  Copyright 2008 Patrick D. Schloss. All rights reserved.
  *
+ *	A sequence object has three components: i) an accession number / name, ii) the unaligned primary sequence, iii) a
+ *	pairwise aligned sequence, and iv) a sequence that is aligned to a reference alignment.  This class has methods
+ *	to set and get these values for the other classes where they are needed.  *
+ *
  */
-using namespace std;
 
-#include "mothur.h"
 
+//Data Structure for a fasta file.
+
+#include "mothur.h"
+#include "mothurout.h"
 
 /**************************************************************************************************/
 
@@ -20,28 +26,59 @@ class Sequence {
 public:
 	Sequence();
 	Sequence(string, string);
+	Sequence(ifstream&);
+	Sequence(istringstream&);
+	Sequence(const Sequence& se) : name(se.name), unaligned(se.unaligned), aligned(se.aligned), pairwise(se.pairwise), numBases(se.numBases), startPos(se.startPos), endPos(se.endPos),
+									alignmentLength(se.alignmentLength), isAligned(se.isAligned), longHomoPolymer(se.longHomoPolymer), ambigBases(se.ambigBases) { m = MothurOut::getInstance(); }
+	
+	//these constructors just set the unaligned string to save space
+	Sequence(string, string, string);  
+	Sequence(ifstream&, string);
+	Sequence(istringstream&, string);
 	
 	void setName(string);
 	void setUnaligned(string);
 	void setPairwise(string);
 	void setAligned(string);
 	void setLength();
+	void reverseComplement();
+	void trim(int);
 	
 	string convert2ints();
 	string getName();
 	string getAligned();
 	string getPairwise();
 	string getUnaligned();
-	int getLength();
+	string getInlineSeq();
+	int getNumBases();
+	int getStartPos();
+	int getEndPos();
+	void padToPos(int);
+	void padFromPos(int);
+	int getAlignLength();
+	int getAmbigBases();
+	void removeAmbigBases();
+	int getLongHomoPolymer();
+	bool getIsAligned();
 	void printSequence(ostream&);
 	
 private:
+	MothurOut* m;
+	void initialize();
+	string getSequenceString(ifstream&, int&);
+	string getCommentString(ifstream&);
+	string getSequenceString(istringstream&, int&);
+	string getCommentString(istringstream&);
 	string name;
 	string unaligned;
 	string aligned;
 	string pairwise;
-	int length;
-	int lengthAligned;
+	int numBases;
+	int alignmentLength;
+	bool isAligned;
+	int longHomoPolymer;
+	int ambigBases;
+	int startPos, endPos;
 };
 
 /**************************************************************************************************/