Can you please guide me how to solve this problem? I'd like to use TCGA data as MuSiC input but don't have controlled-access at a moment. This means that I can't use BAM file as MuSiC input. Do you have any suggestion to circumvent this problem to run 'genome music play'? From the installation guide at Installation of the MuSiC suite on unsupported linux distributions, it suggest using wig file to replace BAM file but it didn't explain how to plug this replacement in 'genome music play'.
Do I have to create my own pipeline? Or I can just need to add a few more parameters to indicate wig file in 'music play'?
•
written
by
jessada115 •
I think you will need to use music bmr calc-wig-covg to generate the gene coverage files.
•
written
by
mdm-two •
This is correct. The "music play" script only works with "music bmr calc-covg" that requires BAM files. If you only have access to wig files, you will need to run "bmr calc-wig-covg", followed by "bmr calc-bmr", and then the "smg" test and other MuSiC modules.
Yes, calc-wig-covg generates a gene_covgs directory and a file named total_covgs which is all the information needed by calc-bmr... which is the necessary next step of either calc-covg or calc-wig-covg.
Log in to answer this question.
