View Single Post
Old 12-12-2016, 05:20 PM   #2
Brian Bushnell
Super Moderator
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707

There are three primary ways. I'll describe how to do them using the BBMap package.

1) Via mapping, which requires a reference: in1=r1.fastq in2=r2.fastq ref=ref.fasta ihist=ihist.txt reads=2m pairlen=2000

2) Via overlap, which requires overlapping reads (they probably overlap given you ran at 2x250): in1=r1.fastq in2=r2.fastq ihist=ihist.txt reads=2m

3) Via assembly, which requires sufficient read depth and memory to assemble the genome: in1=r1.fastq in2=r2.fastq ihist=ihist.txt extend2=200

2 is the fastest. The best choice and best settings depend on your data, though. Can you describe the organism, experiment, and target insert size?
Brian Bushnell is offline   Reply With Quote