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.

Simulating reads to test methods

evolvedmicrobeevolvedmicrobe MGHPosts: 14Member
edited January 2013 in Ask the GATK team

I just wrote a walker to look for particular types of low frequency mutations, and I wanted to verify that the methods were working. I was hoping to simulate some illumina data with the variants and then run the methods against this data.

However, I don't know what a realistic error model is for common Illumina data and so am not sure how realistic my simulations are (Proportion of gaps, A->C versus A->G, etc.). Does the GATK include a read simulator? I saw one walker in the documentation but it seemed to rely on inputting settings that I didn't know about it and looked a bit out of date.

Any help appreciated.

Post edited by Geraldine_VdAuwera on

Best Answer


Sign In or Register to comment.