Hi, I have a simulation in .fasta format, I want to generate simulate read that has 40x coverage, I know the number of read I should generate, but wgsim do not seem to provide such option. Is there any tool I could use to generate such read? Thanks!
3 answers
The -N parameter let's you set the number of reads. Then, in addition, the -1 and -2 flags allow you to set the read length.
The above parameters together provide you with full control over the average coverage of the simulated reads.
ART can simulate outputs of many different sequencing platforms, and has an option to specify coverage.
randomreads.sh from BBMap suite can do coverage based read generation. You can do many other things like add mutations etc. Run the program by itself to see options.
Log in to answer this question.
-
FASTQ Simulation with known mutations
written by Raman2 •
-
What's the best genetic data simulation tool? GPOPSIM, GCTA or sim1000G
written by Friyunnwoo •
-
What is the best tool for generating Multiple sequence alignment (MSA) of fasta sequences ?
written by sunnykevin97
-
Simulation using sim1000G tool
written by katkota.a •
-
Simulation of single-end reads from reference genome
written by newscient •
-
wgsim no output
written by RJ •
-
tool simulate paired-end reads and report where each read come from
written by Chen Sun
-
Simulating Sequences From Simulated Svs At Different Cellular Frequencies
written by Aaronquinlan
-
What Ngs Read Simulators Are Available For Single-End 35Base Pair Length Data?
written by Chintzz •
-
Ngs Dna Read Simulator With Quality Scores Available?
written by Travis
