marbl / Winnowmap

Long read / genome alignment software
Other
235 stars 22 forks source link

Way to filter winnowmap2 BAM files for read lengths #20

Open Musketeer-D opened 3 years ago

Musketeer-D commented 3 years ago

Could I use 'samtools view -h Winnowmap.bam | awk 'length($10) > 50000 || $1 ~ /^@/' | samtools view -bS - > 50k.bam' to get winnowmap2 BAM files filtered for read lengths greater than 50 kb?