提取相应的位置的碱基序列

从maf文件中提取突变位点附件的碱基序列

首先准备bed文件

 

 

bedtools getfasta

bedtools getfasta -fi $fa -bed test.bed -fo test.fa

结果既是附近的碱基