View Single Post
Old 04-06-2017, 10:19 AM   #451
Brian Bushnell
Super Moderator
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707

Normally, I would just map the reads with the flag "maxindel=100k" (or whatever is appropriate, given the length of the genome) and then call variants with "" using the flag "rarity=0.01" to find low-frequency events. But alternatively, you can separate the reads with multiple BBMap passes, like this: in=reads.fq outm=mapped.fq maxindel=100k in=mapped.fq outm=normal.fq outu=longdeletion.fq maxindel=100k dellenfilter=20 in=longdeletion.fq out=mapped_longdeletion.sam maxindel=100k
Brian Bushnell is offline   Reply With Quote