Подвыборка файла BAM для включения только чтений, длина которых превышает X

У меня есть несколько файлов .bam (выровненных с bwa), и я хотел бы сгенерировать файл подмножества .bam, который включает только чтения выше 50bp (длина чтения, а не выровненная длина). Есть ли способ сделать это из файла bam? Или мне следует вернуться и подгруппировать fastq, а затем перестроить?