Question: simulate reads with zero mutations
1
gravatar for jyu429
4.3 years ago by
jyu429120
United States
jyu429120 wrote:

Hi,

I'm interested in simulating reads with a mutation rate of zero. So basically I just want to generate reads of a certain length without any mutations. Is there a way to do this with wgsim or perhaps Art? Thanks!

 

wgsim • 1.4k views
ADD COMMENTlink modified 4.3 years ago • written 4.3 years ago by jyu429120
3
gravatar for Pierre Lindenbaum
4.3 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum121k wrote:
$ wgsim

Program: wgsim (short read simulator)
Version: 0.3.0
Contact: Heng Li <lh3@sanger.ac.uk>

Usage:   wgsim [options] <in.ref.fa> <out.read1.fq> <out.read2.fq>

Options: -e FLOAT      base error rate [0.020]
         -d INT        outer distance between the two ends [500]
         -s INT        standard deviation [50]
         -N INT        number of read pairs [1000000]
         -1 INT        length of the first read [70]
         -2 INT        length of the second read [70]
         -r FLOAT      rate of mutations [0.0010]
         -R FLOAT      fraction of indels [0.15]
         -X FLOAT      probability an indel is extended [0.30]
         -S INT        seed for random generator [-1]
         -h            haplotype mode

->

 

wgsim -e 0.0 -r 0.0 <in.ref.fa> <out.read1.fq> <out.read2.fq>
ADD COMMENTlink written 4.3 years ago by Pierre Lindenbaum121k
0
gravatar for jyu429
4.3 years ago by
jyu429120
United States
jyu429120 wrote:

Thanks! Is there also a parameter to play with read-depth coverage?

ADD COMMENTlink written 4.3 years ago by jyu429120

If you know the size of the genome, which you do, coverage is reads*readLength/geneomeLength, and the number of reads and their length are parameters.

ADD REPLYlink written 4.3 years ago by Brian Bushnell16k

see option '-N'

ADD REPLYlink written 4.3 years ago by Pierre Lindenbaum121k
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: 1665 users visited in the last hour