我正在编写一个代码来使用 R 分析酵母测序数据。我能够使用 BWA 进行对齐，从而使用 samtools 获得 .bam 文件及其 .bai 索引。我的下一步是使用 bcftools 执行变体调用。这...