SEQanswers (
-   Bioinformatics (
-   -   dwgsim to simulate Illumina reads (

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 ( 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 02:04 AM.

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