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

automate exome depth #506

Closed
wants to merge 17 commits into from
Closed

automate exome depth #506

wants to merge 17 commits into from

Conversation

aledj2
Copy link

@aledj2 aledj2 commented Sep 27, 2023

This change is Reviewable

@aledj2 aledj2 changed the title untested changes - still awaiting inputs automate exome depth Oct 6, 2023
Copy link
Contributor

@RachelDuffin RachelDuffin left a comment

Choose a reason for hiding this comment

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

Reviewed 2 of 3 files at r1, 1 of 1 files at r2, all commit messages.
Reviewable status: all files reviewed, 16 unresolved discussions (waiting on @aledj2)


automate_demultiplex_config.py line 11 at r2 (raw file):

# Set debug mode
testing = True

To change at end of testing


automate_demultiplex_config.py line 292 at r2 (raw file):

ED_VCP2_readcount_BEDfile_pannum = "Pan5188_exomedepth.bed" 
#VCP3 normals data file
ED_readcount_normals_VCP3_file=None#"Pan5149_normals_v1.0.0.RData"

Should this have a file?


automate_demultiplex_config.py line 493 at r2 (raw file):

    "Pan4146",  # VCP3 Viapath R81 (CM)
    "Pan4149",  # VCP2 Viapath R208 (BRCA)
    "Pan4150",  # VCP2 Viapath R207 (ovarian)

Was it intentional to remove Pan4127 and Pan4818?


automate_demultiplex_config.py line 609 at r2 (raw file):

    "Pan4146",
    "Pan4151",
    #"Pan4043",

Should this be commented out?


automate_demultiplex_config.py line 767 at r2 (raw file):

    #     "variant_calling_bedfile": "Pan5119data.bed",
    #     "sambamba_bedfile": "Pan5123dataSambamba.bed",
    #     "exome_depth_cnvcalling_BED":  "PanXXXX"

Should these two pan numbers be commented out?


automate_demultiplex_config.py line 820 at r2 (raw file):

    #     "congenica_credentials": "STG",
    #     "congenica_IR_template": "non-priority",
    #     "congenica_project": "4201",

Should these two pan numbers be commented out?


automate_demultiplex_config.py line 1249 at r2 (raw file):

        "sambamba_bedfile": "Pan4995dataSambamba.bed",
        "variant_calling_bedfile": "Pan4995data.bed",
        # "exome_depth_cnvcalling_BED": "Pan5179", bedfile not yet made

Replace with bedfile?


automate_demultiplex_config.py line 1440 at r2 (raw file):

        "variant_calling_bedfile": "Pan4995data.bed",
        "sambamba_bedfile": "Pan4995dataSambamba.bed",
        # "exome_depth_cnvcalling_BED": "Pan5174" BEDfile not yet available

Replace with bedfile?


automate_demultiplex_config.py line 1503 at r2 (raw file):

        "variant_calling_bedfile": "Pan4995data.bed",
        "sambamba_bedfile": "Pan4995dataSambamba.bed",
        #"exome_depth_cnvcalling_BED": "Pan5179" BEDfile not yet available

Replace with bedfile?


automate_demultiplex_config.py line 1597 at r2 (raw file):

        "sambamba_bedfile": "Pan5018dataSambamba.bed",
        "masked_reference": "project-ByfFPz00jy1fk6PjpZ95F27J:file-GF84GF00QfBfzV35Gf8Qg53q",  # hs37d5_Pan4967.bwa-index.tar.gz
        

Remove


automate_demultiplex_config.py line 1610 at r2 (raw file):

        "sambamba_bedfile": "Pan5018dataSambamba.bed",
        "masked_reference": "project-ByfFPz00jy1fk6PjpZ95F27J:file-GF84GF00QfBfzV35Gf8Qg53q",  # hs37d5_Pan4967.bwa-index.tar.gz
        

Remove


automate_demultiplex_config.py line 1780 at r2 (raw file):

    "tso_pannumbers": "-itso_pannumbers=Pan4969,Pan5085,Pan5114",
    "stg_pannumbers": (
        "-istg_pannumbers=Pan4821,Pan4822,Pan4823,Pan4824,Pan4825,Pan4816,Pan4817,Pan4819,Pan4820,"

Why have these pan numbers been removed?


upload_and_setoff_workflows.py line 122 at r2 (raw file):

            + "_TSO_post_run_commands.sh"
        )
        #TODO copy lines above to create separate dx run commands output script for TSO (to be run by duty binfx)

To do or delete


upload_and_setoff_workflows.py line 644 at r2 (raw file):

        take TSO500 samplesheet and split in to parts with <=16 samples/sheet
        write samplesheets to runfolder
        return list of samplesheet paths? or just names (if they're saved in the runfolder, 

Clarify / remove the '?'


upload_and_setoff_workflows.py line 1795 at r2 (raw file):

                command_list.append(self.build_ED_cnv_calling_cmd(panel))

        return command_list

Add a space after the function


upload_and_setoff_workflows.py line 2363 at r2 (raw file):

        return cleaned_list

    # TODO set this up so it only runs the RPKM app if there are enough samples files (minimum 3 required by the app)

To resolve?

Copy link
Author

@aledj2 aledj2 left a comment

Choose a reason for hiding this comment

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

Reviewable status: all files reviewed, 16 unresolved discussions (waiting on @RachelDuffin)


automate_demultiplex_config.py line 292 at r2 (raw file):

Previously, RachelDuffin wrote…

Should this have a file?

this hasn't been made yet (this is why exome depth doesn't support vcp3)


automate_demultiplex_config.py line 493 at r2 (raw file):

Previously, RachelDuffin wrote…

Was it intentional to remove Pan4127 and Pan4818?

yes, this test (R209) has been withdrawn (combined with R211)


automate_demultiplex_config.py line 609 at r2 (raw file):

Previously, RachelDuffin wrote…

Should this be commented out?

have removed


automate_demultiplex_config.py line 767 at r2 (raw file):

Previously, RachelDuffin wrote…

Should these two pan numbers be commented out?

they aren't needed. have removed


automate_demultiplex_config.py line 820 at r2 (raw file):

Previously, RachelDuffin wrote…

Should these two pan numbers be commented out?

they aren't needed. have removed


automate_demultiplex_config.py line 1249 at r2 (raw file):

Previously, RachelDuffin wrote…

Replace with bedfile?

if this variable is commented out exome depth won't be called - this is commented out to flag that it will be turned on in future but what other actions may be needed. For many VCP3 tests these have gone to WGS, and for those that haven't CNV analysis isn't going to be assessed routinely, so creation of these BED files has been low priority.


automate_demultiplex_config.py line 1503 at r2 (raw file):

Previously, RachelDuffin wrote…

Replace with bedfile?

if this variable is commented out exome depth won't be called - this is commented out to flag that it will be turned on in future but what other actions may be needed. For many VCP3 tests these have gone to WGS, and for those that haven't CNV analysis isn't going to be assessed routinely, so creation of these BED files has been low priority.


automate_demultiplex_config.py line 1597 at r2 (raw file):

Previously, RachelDuffin wrote…

Remove

Done.


automate_demultiplex_config.py line 1610 at r2 (raw file):

Previously, RachelDuffin wrote…

Remove

Done.


automate_demultiplex_config.py line 1780 at r2 (raw file):

Previously, RachelDuffin wrote…

Why have these pan numbers been removed?

They are not R number specific, which is needed to perform exome depth


upload_and_setoff_workflows.py line 122 at r2 (raw file):

Previously, RachelDuffin wrote…

To do or delete

Done.


upload_and_setoff_workflows.py line 644 at r2 (raw file):

Previously, RachelDuffin wrote…

Clarify / remove the '?'

Done.


upload_and_setoff_workflows.py line 1795 at r2 (raw file):

Previously, RachelDuffin wrote…

Add a space after the function

Done.


upload_and_setoff_workflows.py line 2363 at r2 (raw file):

Previously, RachelDuffin wrote…

To resolve?

Done.

Copy link
Contributor

@RachelDuffin RachelDuffin left a comment

Choose a reason for hiding this comment

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

Reviewed 2 of 2 files at r3, all commit messages.
Reviewable status: all files reviewed, 1 unresolved discussion (waiting on @aledj2)

Copy link
Contributor

@RachelDuffin RachelDuffin left a comment

Choose a reason for hiding this comment

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

Reviewed all commit messages.
Reviewable status: 2 of 3 files reviewed, 2 unresolved discussions (waiting on @aledj2)


upload_and_setoff_workflows.py line 1571 at r4 (raw file):

                # If panel is to be processed using mokapipe
                if self.panel_dictionary[panel]["mokapipe"]:
                    commands_list.append("#For each sample there are 5 lines of commands. The dx run command for the workflow. and the jobid is then added to two depends_on lists.")

Are these and other comments going to be removed?

Copy link
Contributor

@RachelDuffin RachelDuffin left a comment

Choose a reason for hiding this comment

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

Reviewable status: 2 of 3 files reviewed, 3 unresolved discussions (waiting on @aledj2)


upload_and_setoff_workflows.py line 689 at r4 (raw file):

        number_of_batches = len(batches)
        #capture path for samplesheet in runfolder
        runfolder_samplesheet_file = samplesheet_file = os.path.join(self.runfolder_obj.runfolderpath, self.runfolder_obj.runfolder_samplesheet_name)

Is this a mistake?

Copy link
Contributor

@RachelDuffin RachelDuffin left a comment

Choose a reason for hiding this comment

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

Reviewed 1 of 1 files at r4.
Reviewable status: all files reviewed, 3 unresolved discussions (waiting on @aledj2)

Copy link
Contributor

@rebeccahaines1 rebeccahaines1 left a comment

Choose a reason for hiding this comment

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

Reviewable status: 2 of 3 files reviewed, 3 unresolved discussions (waiting on @aledj2 and @RachelDuffin)


upload_and_setoff_workflows.py line 689 at r4 (raw file):

Previously, RachelDuffin wrote…

Is this a mistake?

yes, changed this now


upload_and_setoff_workflows.py line 1571 at r4 (raw file):

Previously, RachelDuffin wrote…

Are these and other comments going to be removed?

done

Copy link
Contributor

@RachelDuffin RachelDuffin left a comment

Choose a reason for hiding this comment

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

Reviewed 1 of 1 files at r5, all commit messages.
Reviewable status: all files reviewed, 1 unresolved discussion (waiting on @aledj2)

Copy link
Collaborator

@natashapinto natashapinto left a comment

Choose a reason for hiding this comment

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

Reviewed 1 of 1 files at r6, all commit messages.
Reviewable status: :shipit: complete! all files reviewed, all discussions resolved (waiting on @aledj2)

@rebeccahaines1 rebeccahaines1 marked this pull request as ready for review October 23, 2023 10:08
@rebeccahaines1
Copy link
Contributor

PR being closed without merge because set up incorrectly- should have been set to merge feature to development branch

@natashapinto natashapinto mentioned this pull request Oct 23, 2023
@natashapinto natashapinto deleted the feature/exomedepth branch October 23, 2023 10:20
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.

5 participants