for ( racpi  = (*rmpi) -> best . begin( ); racpi != (*rmpi) -> best . end( ); racpi ++ ) { cout << "#\n"; for ( curraln  = (*racpi) -> head; curraln != NULL; curraln = curraln -> from ) cout << ' ' << curraln -> loR << '\t' << curraln -> hiR << endl; }
void PrintMapping (Mapping_t & mapping){  ReadMap_t temprm;  temprm . place = new ReadAlignChain_t( );  temprm . place -> head = new ReadAlign_t( );  map<string, Reference_t>::iterator rmi;  vector<ReadMap_t *>::iterator rmpi;  pair<    vector<ReadMap_t *>::iterator,    vector<ReadMap_t *>::iterator    > rmpip;  //-- For each reference  for ( rmi  = mapping . references . begin( );	rmi != mapping . references . end( ); rmi ++ )    {      //-- Find the reference range      temprm . place -> head -> ref = &(rmi -> second);      rmpip = equal_range (mapping . reads . begin( ), mapping . reads . end( ),			   &temprm, ReadReferenceCmp_t( ));      cout << '>' << rmi -> first << '\t'	   << (rmi -> second) . len	   << endl;      //-- For each read mapping to this reference      for ( rmpi = rmpip . first; rmpi != rmpip . second; rmpi ++ )	{	  cout << (*rmpi) -> id << '\t'	       << (*rmpi) -> place -> beg << '\t'	       << ((*rmpi) -> place -> head -> ori == FORWARD_CHAR ? '+' : '-');	  if ( (*rmpi) -> mate . read == NULL )	    cout << "\tNULL";	  else	    cout << '\t' << (*rmpi) -> mate . read -> id;	  cout << endl;	}    }  delete temprm . place -> head;  delete temprm . place;}