Question: How to Specify Genome Build in MutSigCV Input ?
0
gravatar for H.Hasani
5.3 years ago by
H.Hasani980
Freiburg, Germany
H.Hasani980 wrote:

Hello,

I wonder if there is a way to to specify the genome assembly given to MutSigCV? I would like to use GRCh37.75 as I'm getting the following error :

Error using MutSigCV>MutSig_preprocess (line 542)
probable build mismatch between mutation_file and chr_files

For I've been using the reference files coming with the tool.

I'm very grateful for any quick help

 

assembly snp genome gene • 2.2k views
ADD COMMENTlink modified 3.3 years ago by jz3140 • written 5.3 years ago by H.Hasani980

Hi,

Were you able to resolve this issue? I am running into the same problem now.

Thanks!

ADD REPLYlink written 4.9 years ago by ssekar0

Unfortunately no! I would love to know how!

ADD REPLYlink written 4.9 years ago by H.Hasani980
0
gravatar for jz314
3.3 years ago by
jz3140
jz3140 wrote:

I had this issue as well using the reference files that came with the tool, and it turned out I had a 1-off error- on the Start_position in the input mutation file. It should be 1-based relative to the Reference_Allele.
The way I figured this out was to check the base at specific positions in the provided reference files with commands such as head -c <position> chr1.txt|tail -c 1 and make sure the output matched the Reference_Allele base(s) in my input maf.
If you wanted to use a custom genome build instead of the hg19 provided, your format would probably need to be similar to the provided files, 1 file per chromosome with no linebreaks, so you'd have to convert a typical fasta accordingly (something like cat <your_reference_fasta_chromosome>|grep -v ">"|awk '{printf $1}' > <out_chr.txt>

ADD COMMENTlink modified 3.3 years ago • written 3.3 years ago by jz3140
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2243 users visited in the last hour
_