#include "oechem.h"
#include <iostream>

using namespace OEChem;
using namespace std;

int main()
{
  OEUniMolecularRxn umr("[N:1]>>[Nh3+:1]");

  OEGraphMol mol;
  OEParseSmiles(mol,"NCCN");
  umr(mol);

  std::string smi;
  OECreateSmiString(smi,mol);

  cout << "smiles = " << smi << endl;
  
  return 0;
}