标签：average 测序 基因组 混样 hisat2 depth bam ID
WGS和宏基因组混样测序-判断和host基因组的乘深和覆盖度
数据准备
在前期已经做好map到host基因组的步骤以后 （使用了hisat2）得到的 ID.hisat2.bam文件
第一步:sort
因为我们下一步要用的samtools depth的功能需要sorted的bam文件，操作如下：
samtools sort ID.hisat2.bam ID.hisat2.sorted.bam
第二步：判断depth
简单粗暴，操作如下：
samtools depth ID.hisat2.sorted.bam > ID_depth_result
第三步：求average depth
average_depth=$(awk ‘{sum+=$3}END{print sum/NR}’ ID_depth_result);echo ${average_depth}
第四步：求average coverage
name=“ID”;satisfied=$(awk '{if(KaTeX parse error: Expected 'EOF', got '}' at position 20: …erage1){sum+=1}}̲' average1={average_depth} ID_depth_result）;all=$(wc -l ID_depth_result|awk ‘{print $1}’);echo ${statisfied};echo ${all}
再用statisfiel/all 就是我们需要的average coverage啦
EffeAuuuu
发布了1 篇原创文章 · 获赞 0 · 访问量 8
私信
关注
标签：average,测序,基因组,混样,hisat2,depth,bam,ID 来源： https://blog.csdn.net/EffeAuuuu/article/details/103983360
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。