X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=sequence.hpp;h=95c3e51839029393798dc2f57848d5e325c6c2df;hb=b8f3e71c8a0310e25261464d27e6fcfe891483b0;hp=03cbab7ffddce5531e272cf6116069edfa0ebade;hpb=d5d2761f88b41f1006d0b700e0ab51e2ce48b875;p=mothur.git

diff --git a/sequence.hpp b/sequence.hpp
index 03cbab7..95c3e51 100644
--- a/sequence.hpp
+++ b/sequence.hpp
@@ -8,11 +8,14 @@
  *  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"
-
+#include "mothurout.h"
 
 /**************************************************************************************************/
 
@@ -20,28 +23,52 @@ class Sequence {
 public:
 	Sequence();
 	Sequence(string, string);
+	Sequence(ifstream&);
+	Sequence(istringstream&);
+	
+	//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();
 	
 	string convert2ints();
 	string getName();
 	string getAligned();
 	string getPairwise();
 	string getUnaligned();
-	int getLength();
+	int getNumBases();
+	int getStartPos();
+	int getEndPos();
+	int getAlignLength();
+	int getAmbigBases();
+	int getLongHomoPolymer();
+	bool getIsAligned();
 	void printSequence(ostream&);
 	
 private:
+	MothurOut* m;
+	void initialize();
+	string getSequenceString(ifstream&);
+	string getCommentString(ifstream&);
+	string getSequenceString(istringstream&);
+	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;
 };
 
 /**************************************************************************************************/