bam文件截取指定长度的序列

截取长度在50-100的序列

samtools view -h test.bam | awk 'length($10) > 50 && length($10) < 100 || $1 ~ /^@/' | samtools view -bS - > new.bam

截取长度为50，51的序列

samtools view -h test.bam | awk 'length($10) ==50  || length($10) == 51 || $1 ~ /^@/' | samtools view -bS - > new.bam
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。