Bug Bulletin: The recent 3.2 release fixes many issues. If you run into a problem, please try the latest version before posting a bug report, as your problem may already have been solved.

Mouse VQSR

lacasabonitalacasabonita Posts: 10Member

I was wondering if anyone has used VQSR for a mouse related genome project. I am working with mm10 dbsnp and DNA-seq short insert data for multiple homozygous mouse samples. I have obtained decent results so far using the mm10 dbsnp as the training set, but was curious to see if anyone had any recommendations as to what settings to use. Any input is appreciated. I also have a lot of RNA-seq data, but that will come at a much later point in time. Thanks!

Tagged:

Comments

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 5,869Administrator, GATK Developer admin

    Hi there,

    We do not have any mouse-specific recommendations, but I know we have other users who work on mouse genomes. Hopefully someone will volunteer to share their experience with you.

    In relation to this, we are considering setting up subcategories in the forum to consolidate organism-specific discussions. Is that something you would want to use to connect with other researchers who also apply GATK tools to the organism you study?

    Geraldine Van der Auwera, PhD

  • lacasabonitalacasabonita Posts: 10Member

    Yeah I think that would be very helpful. I have another question though. What are the key differences between SNP mode and INDEL mode when running VQSR? I have run in both with similar settings, and the number of variants filtered are almost identical. I have also run in BOTH mode and a significant amount of variants were left unfiltered compared to the other two modes.

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 5,869Administrator, GATK Developer admin

    The main thing you should consider when deciding which mode to use for VQSR is how your variants were called in the first place. If you called them together using the same caller, applying the same statistical model, then it is okay to filter them together too. But if you called them separately then merged the results, you shouldn't filter using BOTH, because some of the metrics will not be comparable directly. When in doubt I would recommend performing the procedure separately.

    Geraldine Van der Auwera, PhD

  • lacasabonitalacasabonita Posts: 10Member
Sign In or Register to comment.