Skip to content

Commit

Permalink
Validate median insert size
Browse files Browse the repository at this point in the history
  • Loading branch information
seallard committed Dec 12, 2023
1 parent 418837d commit 4cabeca
Show file tree
Hide file tree
Showing 3 changed files with 10 additions and 0 deletions.
1 change: 1 addition & 0 deletions cg/constants/constants.py
Original file line number Diff line number Diff line change
Expand Up @@ -214,6 +214,7 @@ class MicrosaltQC:
QC_PERCENT_THRESHOLD_MWX: float = 0.1
COVERAGE_10X_THRESHOLD: float = 0.75
DUPLICATION_RATE_THRESHOLD: float = 0.8
INSERT_SIZE_THRESHOLD: int = 100
MAPPED_RATE_THRESHOLD: float = 0.3
NEGATIVE_CONTROL_READS_THRESHOLD: float = 0.2
TARGET_READS: int = 6000000
Expand Down
5 changes: 5 additions & 0 deletions cg/meta/workflow/microsalt/quality_checker.py
Original file line number Diff line number Diff line change
Expand Up @@ -8,6 +8,7 @@
from cg.meta.workflow.microsalt.utils import (
is_valid_duplication_rate,
is_valid_mapped_rate,
is_valid_median_insert_size,
is_valid_total_reads,
is_valid_total_reads_for_control,
parse_quality_metrics,
Expand Down Expand Up @@ -166,3 +167,7 @@ def is_valid_mapped_rate(self, metrics: SampleMetrics) -> bool:
def is_valid_duplication_rate(self, metrics: SampleMetrics) -> bool:
duplication_rate: float = metrics.picard_markduplicate.duplication_rate
return is_valid_duplication_rate(duplication_rate)

def is_valid_median_insert_size(self, metrics: SampleMetrics) -> bool:
insert_size: int = metrics.picard_markduplicate.insert_size
return is_valid_median_insert_size(insert_size)
4 changes: 4 additions & 0 deletions cg/meta/workflow/microsalt/utils.py
Original file line number Diff line number Diff line change
Expand Up @@ -21,6 +21,10 @@ def is_valid_duplication_rate(sample_duplication_rate: float) -> bool:
return sample_duplication_rate < MicrosaltQC.DUPLICATION_RATE_THRESHOLD


def is_valid_median_insert_size(sample_insert_size: int) -> bool:
return sample_insert_size > MicrosaltQC.INSERT_SIZE_THRESHOLD


def parse_quality_metrics(file_path: Path) -> QualityMetrics:
data = read_json(file_path)
return QualityMetrics.model_validate_json(data)

0 comments on commit 4cabeca

Please sign in to comment.