sam/bam文件转换成vcf

有没有了解samtools或者bcftools的，怎么用这些工具将sam/bam文件转换成vcf格式

bam2fasta的转变方式：
samtools view SL3003_SL3004_100122-hg18.bam | \
awk '{OFS="\t"; print ">"$1"\n"$10}' - > SL3003.fasta

sam2fasta的转变方式
cat *.sam | awk '{print ">"$1"\n"$10}' > *.fasta

1、sort sort对bam文件进行排序。 usage: samtools sort [-n] [-m ] -m 参数默认下是 500,000,000 即500m（不支持k，m，g等缩写）。对于处理大数据时，如果内存够用，则设置大点的值，以节约时间。 -n 设定排序方式按short reads的id排序。

如以上回答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！

点此我要举报以上问答信息