script for interrogating PacBio sequencing chemistry used

showChem.py

#!/usr/bin/env python
import sys
import h5py

fname = sys.argv[1]
f = h5py.File(fname, "r")

if fname.endswith("bax.h5"):
    ri = f["/ScanData/RunInfo"]
    try:
        chem = ri.attrs["SequencingChemistry"]
        print sys.argv[1], chem
    except:
        print sys.argv[1], "absent"
else:
    assert fname.endswith("cmp.h5")
    for (movieName, chemistry) in zip(f["/MovieInfo/Name"],
                                      f["/MovieInfo/SequencingChemistry"]):
        print movieName, chemistry