TBtools攻略试读

TBtools | 基因组重测序数据分析常用上游软件「三兄弟」,

2022-06-03  本文已影响0人  生信石头

写在前面

是的,最后,我选择了「一条路走到黑」~ 咱们先不管一些人如何看,东西先整出来。大体如下。


image.png

没错。从「原始测序数据」到「SNP Calling」,也就是 从 .fastq 到 .vcf,一套聚齐了。其中涉及到三个目前最常用的软件:

  1. bwa-mem2 -> 读段回帖
  2. samtools -> 重复标记排序
  3. bcftools -> 变异检测过滤

有了这三个插件,所有的「TBtools」用户,可以在 Windows 或者 MacOS 下,点点点,分析自己想分析的数据。当然,时间成本前面有提到,最好是按照至少一天跑一个样品来算。做做 BSAseq 或者少数两三个样品还是可以的。大的群体,那么就可能需要足够好配置的机器....这个暂时不是咱们要考虑的事情。
前两个插件都介绍过了。今天介绍「BCFtools GUI Wrapper」,整体如下。



一次包括两个功能,变异检测和变异过滤。

使用「TBtools」进行「SNP Calling」

从界面来看,比较简单。只需要将 MarkDup.PositionSorted.bam 都拖拽放置进去就可以了。



逻辑上需要等一些时间,可以得到初步的 bcf 文件。


使用「TBtools」进行「SNP Filtering」

「SNP Calling」之后,几乎所有情况下,都会做下一步过滤。具体过滤参数会对后续分析产生一定的影响。此处,默认给了一个相对严格的过滤参数,尽量保留「真阳性」SNP。具体使用方法如下。



简单来说...就是把文件放进去,然后设置一下输出目录,最后点击「Start」即可。


写在最后

基本完事。三个「玩具」,欢迎「试玩」。我觉得,BSAseq真的可以试试了。因为有了 VCF 之后,剩下的其实就是一个 R包 的事儿。这个就用 R-Plugin 了。而如果你的 PC 真的强,你甚至还可以用搞 GWAS,因为有了 VCF,剩下的不就是? R 包吗?

上一篇下一篇

猜你喜欢

热点阅读