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

Not found: reference/mutation_type_dictionary.v4.txt #25

Open
ulloaa opened this issue Dec 19, 2023 · 8 comments
Open

Not found: reference/mutation_type_dictionary.v4.txt #25

ulloaa opened this issue Dec 19, 2023 · 8 comments

Comments

@ulloaa
Copy link

ulloaa commented Dec 19, 2023

Hello,

We have downloaded hg38 reference files from:

gs://getzlab-mutsig2cv-reference-files/hg38

and have followed instructions shown in:

https://github.com/getzlab/MutSig2CV/issues/10

However, we still get the error:

Not found: reference/mutation_type_dictionary.v4.txt
Could you please advise?

Thank you!

@julianhess
Copy link
Contributor

Hi,

The mutation type dictionary is the same between hg19 and hg38. Please use the file from the hg19 bundle.

—Julian

@ulloaa
Copy link
Author

ulloaa commented Dec 20, 2023

thanks so much for your prompt response; what about the other required files that are not included in the hg38 reference?

such as:
'reference/mutation_type_dictionary.v4.txt' 'reference/coverage_models.v5a.mat' 'reference/coverage_basewise.fwb' 'reference/target_list.hg19.v1a.txt' 'reference/context_and_effect.c65e29b.fwb' 'reference/context_and_effect.c65e29b.txt' 'reference/covariates_transformed.v5a.txt' 'reference/conservation46.fwb' 'reference/FixedWidthBinary.jar'

Thanks again!

@julianhess
Copy link
Contributor

Some of the files have been renamed between hg19 and hg38. Following the instructions in this post should give you everything you need.

@sabertadros
Copy link

Hi,
We ran standalone version following exactly the steps as illustrated in the homepage.

bin/MutSig2CV <path to mutations> <path to output directory> [params file]

The parameter file we used are adapted from issue 10, as followed:

skip_permutations true
mutation_type_dictionary_file reference/mutation_type_dictionary.v4.txt
coverage_models_mat_file       reference/coverage_models.agilent.mat
target_list_file        reference/target_list.gencode25.txt
context_and_effect_fwb_file     reference/context_and_effect_track.fwb
context_and_effect_fwi_file     reference/context_and_effect_track.fwi
context_and_effect_categs_file  reference/context_and_effect_dict.txt
covariates_file        reference/covariates.agilent.txt

We used hg38 files from GCP, for the rest necessary files we used those from hg19 as you instructed above. And when we ran the program, the error messages are shown as following:

2CV v3.11
Warning: "FixedWidthBinary.jar" is already specified on static java path.
> In javaclasspath>local_validate_dynamic_path at 293
  In javaclasspath>local_javapath at 182
  In javaclasspath at 102
  In MutSig_2CV_v3_11_core at 27
  In MutSig_2CV_v3_11_wrapper at 51 
LOADING DATA
Error using MutSig_2CV_v3_11_core (line 90)
Java exception occurred: 
java.lang.Exception: File and index are mismatched, no matter what width is tried
at org.broadinstitute.cga.tools.seq.FixedWidthBinary.open(FixedWidthBinary.java:170)
at org.broadinstitute.cga.tools.seq.FixedWidthBinary.open(FixedWidthBinary.java:148)
at org.broadinstitute.cga.tools.seq.FixedWidthBinary.open(FixedWidthBinary.java:144)
at org.broadinstitute.cga.tools.seq.FixedWidthBinary.<init>(FixedWidthBinary.java:113)

Error in MutSig_2CV_v3_11_wrapper (line 51)

We wonder if the error might be related to the rest of the reference files including:

conservation46.fwb 
context_and_effect_track.fwi 
coverage_basewise.fwb
coverage_basewise.fwi        
mutation_type_dictionary.v4.txt
conservation46.fwi 

What do you recommend? Thanks

@julianhess
Copy link
Contributor

Hi,

This is coming from an issue with coverage_basewise.fwb/coverage_basewise.fwi. Is it possible that these files are not being specified properly (or corrupted)?

—Julian

@sabertadros
Copy link

Thanks for your response.
We ran the program specifying these two files in the params.txt and it returned with the same exact errors as in the previous post. Next, we ran the program with hg19 built using the same files and they worked fine. So, it seems that neither path specification nor file corruption is the cause for the error.

skip_permutations true
mutation_type_dictionary_file reference/mutation_type_dictionary.v4.txt
coverage_models_mat_file       reference/coverage_models.agilent.mat
target_list_file        reference/target_list.gencode25.txt
context_and_effect_fwb_file     reference/context_and_effect_track.fwb
context_and_effect_fwi_file     reference/context_and_effect_track.fwi
context_and_effect_categs_file  reference/context_and_effect_dict.txt
basewise_coverage_fwb_file  reference/coverage_basewise.fwb
basewise_coverage_fwi_file  reference/coverage_basewise.fwi
covariates_file        reference/covariates.agilent.txt

Since we are using these 6 files from the hg19 built and comparable files for hg38 built are not included in the cloud link [as per issue 10], we believe that compatible reference files (coverage_basewise.fwb/coverage_basewise.fwi) for built hg38 if they are available might fix the problem.

What do you recommend?

@sunyl312
Copy link

coverage_basewise

I encounter the same problem. it worked well while using the hg19, but encounter errors in hg38. Have you solved the problem? how can you solve it? Thank you very much.

@sabertadros
Copy link

sabertadros commented Jun 20, 2024 via email

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

No branches or pull requests

4 participants