SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   dwgsim to simulate Illumina reads (http://seqanswers.com/forums/showthread.php?t=12558)

gene coder 07-06-2011 11:57 PM

dwgsim to simulate Illumina reads
 
Hello, I want to use dwgsim to simulate Illumina reads from the 1000 Genomes Project. I am particularly interested in individual NA12878.

Does anybody know what settings are (approximately or exactly) correct when it comes to the flags -e, -E, -r -R, -X and -Y (http://sourceforge.net/apps/mediawik...ome_Simulation)? Do you know where I can find technical details (does Illumina publish that?) so I can guess myself.

Can anybody explain those flags to a novice? Sequencing is not my domain.

nilshomer 07-07-2011 09:51 AM

Do you have someone who has an Illumina machine that you could ask for typical settings for -e and -E (base error rate)? For human, the rate of mutation (SNPs) is 1/1000, while for indels it is ~1/10000.

If you are using this tool, sequencing is your domain. -e/-E are the per base mean error rate, -r is the rate of mutation, -R is the fraction of mutations that are indels, and -X controls the indel length distribution (probability of extending an indel).


All times are GMT -8. The time now is 07:07 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.