Overlapping paired end read simulator
I am trying to create some synthetic datasets. The reads in the datasets need to be 150bp and paired end. I want the read pairs to overlap to some extent. I have been using MetaSim but it does not seem to be able to make overlapping reads. Does anyone have any suggestions? NeSSM appears to allow for a gap but not an overlap.
• 2,329 views
•
link
updated
by
Ram
•
written
by
eheizer31 •
•
written
by
eheizer31 •
1 answer
BBMap's RandomReads tool [randomreads.sh] will do that. For example:
randomreads.sh ref=ecoli.fasta out=synth.fq.gz paired interleaved reads=100k length=150 mininsert=200 maxinsert=400 adderrors gaussian
That will result in insert sizes from 200bp to 400bp, with a Gaussian distribution, so about 50% of the reads will overlap.
• 0 views
•
link
updated
by
Ram
•
written
by
Brian Bushnell
•
written
by
Brian Bushnell
Log in to answer this question.
More posts like this
-
Merge overlapping paired end reads from BAM file.
written by chrisk •
-
Can ART reads generator create NextSeq 550 synthetic reads?
written by Tawny
-
Rsamtools mpileup discrepancy
written by Noushin N
-
How would you filter PCR duplicates for merged paired-end reads
written by hohoku •
-
How to quantify the overlapping reads in paired-end DNA sequencing to check the sequencing efficien…
written by cjgunase •
-
bed file intersection overlapping regions in at least x files
written by User6891
-
How to count the reads that overlap with 3' splice sites or 5' splice sites?
written by vaslanzadeh •
-
Has any one created a MetaSIM error profile file (.mconf) for Illumina reads >100bp?
written by Tawny
-
Assembly with Abyss question
written by khmkhm87 •
-
How Reliable Are Statistics Derived From "Subsampling" Of Paired-End Illumina Reads ?
written by thecuriousbiologist
