I'll be giving a course about NGS and I'd like the students hunting for some well defined non-synonymous variations in the domain of a given protein .
My 1st idea : duplicating one chromosome , inserting a mutation and let samtools wgsim generating the (heterozygous) mutations.
2nd idea: changing the source code of wgsim . However it's not clear to me where I should put my snippet of C code to insert the mutation; Do you have any idea ? Thanks.