I wanted to simulate PE reads with wgsim but by using the default params, i didn't realise that I was creating reads with 500 bp insert. which doesn't go well with Bowtie's PE default of 250 bp which to cut a long story short, resulted in low pairing rates.
Only found this out by googling biostar of course there is the obligatory discussion of merits of using bowtie versus bwa or vice versa. and there are some nice graphs in there that illustrate the problem at hand.
Thank you SO much for posting this! I had exactly the same problem, and I would've wasted hours debugging it if you hadn't posted this.
ReplyDelete