Skip to content
This repository was archived by the owner on Feb 26, 2025. It is now read-only.
This repository was archived by the owner on Feb 26, 2025. It is now read-only.

tpmean gives somewhat unintuitive results when coverage is low #45

@wwood

Description

@wwood

Since there's so many 0 values in the pileups, given a true coverage of 0.11 and an assumption that no read overlaps another, the tpmean coverage will be 0.01 / (0.8*length_of_contig), which is <<0.11 or even <<0.1 depending on contig length.

This could be fixed by, instead of chopping the top and bottom 10% off, by setting the top 10% of bases to the 90th percentile, and the bottom 10% to the 10th percentile.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions