Question: A Problem About Simulate Colour Reads With Wgwin Tool
0
gravatar for User 0473
6.5 years ago by
User 04730
User 04730 wrote:

I use wgsim with -c parameter to simulate solid reads. Why does it always produce base reads?

could anyone help me to resolve this problem?

thanks

simulation solid • 1.1k views
ADD COMMENTlink modified 3.6 years ago by Biostar ♦♦ 20 • written 6.5 years ago by User 04730
1
gravatar for Ben Lange
6.2 years ago by
Ben Lange180
Minneapolis, MN
Ben Lange180 wrote:

This may be from the wrong version of wgsim, but I'm not seeing a -c command line argument. https://github.com/lh3/wgsim/blob/master/wgsim.c

fprintf(stderr, "\n");
fprintf(stderr, "Program: wgsim (short read simulator)\n");
fprintf(stderr, "Version: %s\n", PACKAGE_VERSION);
fprintf(stderr, "Contact: Heng Li <lh3@sanger.ac.uk>\n\n");
fprintf(stderr, "Usage:   wgsim [options] <in.ref.fa> <out.read1.fq> <out.read2.fq>\n\n");
fprintf(stderr, "Options: -e FLOAT      base error rate [%.3f]\n", ERR_RATE);
fprintf(stderr, "         -d INT        outer distance between the two ends [500]\n");
fprintf(stderr, "         -s INT        standard deviation [50]\n");
fprintf(stderr, "         -N INT        number of read pairs [1000000]\n");
fprintf(stderr, "         -1 INT        length of the first read [70]\n");
fprintf(stderr, "         -2 INT        length of the second read [70]\n");
fprintf(stderr, "         -r FLOAT      rate of mutations [%.4f]\n", MUT_RATE);
fprintf(stderr, "         -R FLOAT      fraction of indels [%.2f]\n", INDEL_FRAC);
fprintf(stderr, "         -X FLOAT      probability an indel is extended [%.2f]\n", INDEL_EXTEND);
fprintf(stderr, "         -S INT        seed for random generator [-1]\n");
fprintf(stderr, "         -h            haplotype mode\n");
ADD COMMENTlink modified 3.6 years ago by Istvan Albert ♦♦ 71k • written 6.2 years ago by Ben Lange180
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1216 users visited in the last hour