If you plan to use this tutorial in an interactive fashion, then you should download an example run folder for an Illumina sequencer. For example, you may download and extract this example run folder: MiSeqDemo
Please change the path below so that it points at the run folder you wish to use:
In [7]:
run_folder = r"D:\RTA.Data\InteropData\MiSeqDemo"
In [8]:
from interop import py_interop_run_metrics, py_interop_run
In [9]:
run_metrics = py_interop_run_metrics.run_metrics()
valid_to_load = py_interop_run.uchar_vector(py_interop_run.MetricCount, 0)
In [10]:
valid_to_load[py_interop_run.Extraction]=1
In [11]:
run_metrics.read(run_folder, valid_to_load)
In [12]:
extraction_metrics = run_metrics.extraction_metric_set()
In [13]:
print "Last Cycle: ", extraction_metrics.max_cycle()
In [14]:
lane, tile, cycle, channel = (1, 1101, 15, 0)
extraction_metric = extraction_metrics.get_metric(lane, tile, cycle)
intensity = extraction_metric.max_intensity (channel)
In [15]:
print "Intensity for channel", channel, " lane", lane, "tile", tile, "cycle", cycle, "is", intensity
In [16]:
extraction_metric = extraction_metrics.at(0)
intensity = extraction_metric.max_intensity (channel)
In [17]:
print "Intensity for channel", channel, " lane", extraction_metric.lane(), "tile", extraction_metric.tile(), "cycle", extraction_metric.cycle(), "is", intensity
In [ ]: