华中农业大学<br>园艺植物生物学教育部重点实验室
IP属地:香港
请问这个脚本~/hychao/script/get_minimap_header_from_fai.py在哪里找呀
全长转录组更新########更新版 minimap2生成的sam文件是没有header的,那使用samtools转bam并排序的时候会报错,这一部分可以自己手动加上
########更新版 minimap2生成的sam文件是没有header的,那使用samtools转bam并排序的时候会报错,这一部分可以自己手动加上
本文中vcf文件的获取主要是使用gatk的HaplotypeCaller模块,该模块的原理大概如下: 定义active regions 沿着参考基因组以一定的窗口滑动,统计比...