Skip to content

Instantly share code, notes, and snippets.

@avalind
Last active February 23, 2023 07:23
Show Gist options
  • Save avalind/f9d0930e8ca0aacbef745f905e4e1f7d to your computer and use it in GitHub Desktop.
Save avalind/f9d0930e8ca0aacbef745f905e4e1f7d to your computer and use it in GitHub Desktop.
MSMutect notes for running

MSMutect notes

I've run in to the following issues while trying to get MSMutect to work:

  1. the path to the data folder in shell_for_all.sh needs to be absolute.
  2. make a symlink with link name phobos in bin, pointing to the correct phobos binary.
  3. bowtie2 needs to have its indices specified by the -x flag, so this needs to be added to the awk-oneliner invoking bowtie2 for each repeat.
  4. Even when adding the -x flag bowtie2 complains about some missing reference files (using the Normal.bam) test data.
  5. on line 211 in Python_managing_file.py vec[:,fi[0]] would sometimes (also on the Normal.bam test data) raise an IndexError, when wrapping this in a try-catch clause the shell_for_all.sh would run without problems for the Normal.bam and Tumor.bam and ms_indel_calling.sh would also run without problems but the resulting maf_like-files would be all empty.

Internal UPPMAX notes

To setup a fitting environment on the milou cluster;

module load bioinfo-tools BEDTools/2.26.0 R/2.15.2 java/sun_jdk1.8.0_92 bowtie2/2.3.2 gcc/4.8 python/2.7

@guodudou
Copy link

Hello,

I am so excited to know that you have successfully run MSMuTect. I am trying to apply MSMuTect and MSMutSig to WES data. But I have several concerns and questions. Can I discuss with you? There is very few support/discussion online regarding these tools.

Thank you very much! I am looking forward to your response,
Wendy

@mattoslmp
Copy link

Dear,
I also have several doubts about these tools, and have no documentation. I would also like to discuss how to successfully run msmutsig. For example, I was able to run MSmutect, but I can not reproduce with successfully MSMutsig for WG data and others, the P-value values obtained weren't significant for samples that I previously know are MSI-High. Thanks to all who can help, I am looking forward to your response too,

@cboursnell
Copy link

What data did you put into MSMutSig? The R script doesn't have any readme or help. There are two blank variables in the file for the A and C motif data, but I don't know what files to put into here. I have tried the outputs from MSMutect but that just gives errors. I have emailed the author, but have had no reply.

@dpandey123
Copy link

Hi Anders,

I need your assistance to run MSMuTect can you please guide me on how to do it? Also, I wanted to know if this will do paired-end.

command used:
sh bin/shell_for_all.sh Normal /home/project/11003021/tools/MSMuTect1-master/data Normal.bam

error screenshot is attached below. KIndly see and please let me know how to resolve it.

image

Thanks

I'm eagerly awaiting your response.

@dpandey123
Copy link

Dear, I also have several doubts about these tools, and have no documentation. I would also like to discuss how to successfully run msmutsig. For example, I was able to run MSmutect, but I can not reproduce with successfully MSMutsig for WG data and others, the P-value values obtained weren't significant for samples that I previously know are MSI-High. Thanks to all who can help, I am looking forward to your response too,

Hi Mattos,

Can you assist me with how to use this tool? When I am using it it's throwing an error with a missing index file. Can you please tell me how you do the troubleshooting?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment