Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Rseq bamstat #155

Merged
merged 18 commits into from
Oct 26, 2024
Merged

Rseq bamstat #155

merged 18 commits into from
Oct 26, 2024

Conversation

emmarousseau
Copy link
Collaborator

Description

Add RSeQC bamstat component

Checklist before requesting a review

  • I have performed a self-review of my code

  • Conforms to the Contributing guidelines

  • Proposed changes are described in the CHANGELOG.md

  • I have tested my code with viash ns test --parallel -q <name or namespace>

  • Check the correct box. Does this PR contain:

    • Breaking changes
    • New functionality
    • Major changes
    • Minor changes
    • Documentation
    • Bug fixes

@emmarousseau emmarousseau marked this pull request as ready for review September 24, 2024 12:50
@KaiWaldrant KaiWaldrant self-requested a review September 27, 2024 08:59

engines:
- type: docker
image: ubuntu:22.04
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
image: ubuntu:22.04
image: python:3.10

Comment on lines 50 to 51
- type: apt
packages: [ python3-pip ]
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
- type: apt
packages: [ python3-pip ]

src/rseqc/rseqc_bamstat/config.vsh.yaml Show resolved Hide resolved
argument_groups:
- name: "Input"
arguments:
- name: "--input"
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Although the input parameter works I would use the parameter defined in the help page as the input parameter name.

Suggested change
- name: "--input"
- name: "--input_file"

type: file
required: true
description: Input alignment file in BAM or SAM format.
- name: "--map_qual"
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
- name: "--map_qual"
- name: "--mapq"


bam_stat.py \
--input "${par_input}" \
${par_map_qual:+--mapq "${par_map_qual}"} \
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
${par_map_qual:+--mapq "${par_map_qual}"} \
${par_mapq:+--mapq "${par_mapq}"} \

echo "> Running $meta_functionality_name."

"$meta_executable" \
--input "$meta_resources_dir/test_data/$input_bam" \
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
--input "$meta_resources_dir/test_data/$input_bam" \
--input_file "$meta_resources_dir/test_data/$input_bam" \

"$meta_executable" \
--input "$meta_resources_dir/test_data/$input_bam" \
--output "$output_summary" \
--map_qual 50
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
--map_qual 50
--mapq 50

echo "> Running $meta_functionality_name."

"$meta_executable" \
--input "$meta_resources_dir/test_data/$input_bam" \
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
--input "$meta_resources_dir/test_data/$input_bam" \
--input_file "$meta_resources_dir/test_data/$input_bam" \


echo ">>> Test 2: Test with non-default mapping quality threshold"

output_summary="mapping_quality_mapq_30.txt"
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
output_summary="mapping_quality_mapq_30.txt"
output_summary="mapping_quality_mapq_50.txt"

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Maybe we can subsample to half to get the file to about 10kb

@emmarousseau
Copy link
Collaborator Author

Thanks Kai, changes have been made!

@rcannood rcannood merged commit 2f8bf02 into viash-hub:main Oct 26, 2024
1 check passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants