售后

脚本路径：

/TJPROJ6/RNA_SH/personal_dir/wanghongyu/script/consensus/get_consensus_work.sh

参数设置：

#第一个参数：vcf文件
#第二个参数：fa文件
#第三个参数：gene文件,包含chr start end gene_name四列
#第四个参数：需要提取的样本名
#第五个参数：输出路径

结果展示

每个样本每个基因的fa

Scaffold00452:35365-41699

AAAAAACTGAGGTTGGTGACTGTGTTCACCAGCAACATCAAAGCACTCTCTCTCTCTCTC TCTCGTCTCTTGCTCAGCTTCTACCCAGTTCTCTCTCTGCATACAAACATATTAATTATC ATCATTCATTTCAAATTTTCAATTCACGCAATTGGATTTGAAATCTGAAATGATGCAGAG