#include <string>\r
#include <map>\r
#include <utility> // for pair\r
+#include <sstream>\r
\r
using namespace std;\r
\r
// index mapping reference/position pairings to bamreaders and their alignments\r
typedef multimap<pair<int, int>, pair<BamReader*, BamAlignment*> > AlignmentIndex;\r
\r
+\r
class BamMultiReader {\r
\r
// constructor / destructor\r