regex: samtools command to “refine” PacBio IsoSeq data?

The samtools view command you’ve provided just uses a filter expression to discard records that do not meet the specified expression criteria. In this case, the expression is looking for alignments where the value of the ‘rq’ (read quality) tag is greater or equal to 0.9. The isoseq tags documentation says that the read quality of the alignment is really the predicted accuracy of the polished isoform. Note that the command produces BAM output, so be sure to redirect stdout or write to a file for example using:

samtools view -b -e '[rq] >= 0.9' -o out.bam flnc.bam

Read more here: Source link