Skip to content

Losing Precision between InterOp and SAV, Python docs #355

@mdooley93-lab

Description

@mdooley93-lab

I have been working to get PF Read information from InterOp and have had trouble getting results that match what SAV is reporting. For example, on Lane 1 of my most recent test of this, InterOp gives me 6372801024.0 while SAV is showing 6,372,801,242. I have had trouble tracking down more information on using InterOp, as well, as there doesn't seem to be Python documentation readily available here.

`# read in sequencing run
run_metrics = iop.read(run_dir)
print("Run data read in")

index_summary_lanes = iop.index_summary(run_metrics) # summary by lane read in as ndarray
indexing_out = iop.indexing(run_metrics) # sample indices as ndarray

print(f"Index summary by lane results:\n{index_summary_lanes}\n")`

Gives me:

Run data read in
Index summary by lane results:
[(1, 0.206 , 0.3533, 0., 95.0816, 6.3728010e+09, 7.6600443e+09)
(2, 0.2037, 0.2388, 0., 95.1361, 6.4046075e+09, 7.6600443e+09)
(3, 0.3738, 0.2415, 0., 94.9932, 6.3131494e+09, 7.6600443e+09)
(4, 0.4561, 0.2674, 0., 94.3576, 6.4385106e+09, 7.6600443e+09)]

My ultimate goal is to export the per lane PF reads along with the per sample % PF reads as they appear in SAV.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions