A Problem About Simulate Colour Reads With Wgwin Tool
1
0
Entering edit mode
10.8 years ago
User 0473 • 0

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

solid simulation • 2.0k views
ADD COMMENT
1
Entering edit mode
10.6 years ago
Ben Lange ▴ 190

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 COMMENT

Login before adding your answer.

Traffic: 2602 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6