wwylab / MuSE

Somatic point mutation caller
GNU General Public License v2.0
18 stars 6 forks source link

Minimal output from MuSE call #6

Closed deniando closed 1 year ago

deniando commented 1 year ago

Hi,

I do not think I am getting the expected output from MuSE call (see attached). It runs without errors but the output contains minimal information. How important is it that we perform the exact pre-processing steps you have outlined? We have our own pre-processing pipeline and have not used all the same software for each step of the pre-processing pipeline that you have outlined. I assumed MuSE would still work with our processed BAMs but maybe that is not the case?

Thanks,

Denise

12-01_S3_L001.MuSE.txt

jiyunmaths commented 1 year ago

@deniando Sorry for addressing your issue late. I think the reason for the empty MuSE.txt file is that you use a compressed reference genome file "GRCh38.primary_assembly.genome.fa.gz". Please can you uncompress it (gzip -d GRCh38.primary_assembly.genome.fa.gz) and run MuSE call again? Let me know if you still have the problem. The preprocessing steps are important, please follow the tutorial. Thanks.

deniando commented 1 year ago

That was the issue @jiyunmaths, thanks very much for your help.