Generate Uniform, Perfect Short Reads From Genome

1

Entering edit mode

10.3 years ago

Raygozak ★ 1.4k

HI, i would like to generate short reads that cover all the genome and that have no mutations, indels, or other artifacts from a fasta file genome.

I tried wgsim with the following command line and it still generated reads with some indels

wgsim -e 0 -d 800 -N 700000 -1 270 -2 270 -r 0 -R 0 -X 0 alignments/SS52_FINAL.fasta ss52_1.fastq ss52_2.fastq

Regards

• 2.2k views

ADD COMMENT • link 10.3 years ago by Raygozak ★ 1.4k

0

Entering edit mode

why don't you try metasim? It's for metagenomics but can also generate reads from a single genome. And you can choose option that applies non errors, and uniform length

ADD REPLY • link 10.3 years ago by matija.sosic ▴ 110

0

Entering edit mode

I suggest you send a report to the developer (maybe a bug?). From what I see, seems like you correctly set to zero all the possible parameters introducing SNPs or indels. I don't know metasim, suggested by @matija.sosic, but it might be worth a try. Also, is it possible that for some strange reason some reads that could be mapped perfectly to a position are erroneously mapped to a different position with an indel? It shouldn't be possible, but it also depends on the aligner...