小编给大家分享一下如何提取bam/sam文件指定区域reads，相信大部分人都还不怎么了解，因此分享这篇文章给大家参考一下，希望大家阅读完这篇文章后大有收获，下面让我们一起去了解一下吧！若想要从sam或bam文件中提取指定区域内的reads，可以使用samtools和bedtools来实现。首先准备一个区域信息文件。region.bed #第一列为染色体ID，第二三列分别为起始终止位置 例： 12 2免费云主机域名1100 41200 #取12号染色体21100-41200区间的序列之后如果是sam文件,先转成bam文件：samtools view -Sb reads.sam > reads.bam然后就可以使用bedtools来提取reads啦，命令如下：bedtools intersect -a reads.bam -b region.bed > target_reads.bam
以上是“如何提取bam/sam文件指定区域reads”这篇文章的所有内容，感谢各位的阅读！相信大家都有了一定的了解，希望分享的内容对大家有所帮助，如果还想学习更多知识，欢迎关注百云主机行业资讯频道！
这篇文章将为大家详细讲解有关微信小程序swiper怎么用，小编觉得挺实用的，因此分享给大家做个参考，希望大家阅读完这篇文章后可以有所收获。 滑块视图容器。 从公共库v1.4.0开始，change事件返回detail中包含一个source字段，表示导致变更的原因…
免责声明：本站发布的图片视频文字，以转载和分享为主，文章观点不代表本站立场，本站不承担相关法律责任；如果涉及侵权请联系邮箱：360163164@qq.com举报，并提供相关证据，经查实将立刻删除涉嫌侵权内容。
微信扫一扫 
