I am trying to extract long reads from the Blasr .sam output. However, I want to collect the entire raw reads that produced alignments, and not just the aligned segments.
Is there an easy way to do this?
Here is the blasr command I am using:
blasr <reads.fasta> <reference.fasta> -noSplitSubreads -clipping none -sam -out output.sam