0
点赞
收藏
分享

微信扫一扫

使用bwa 批量比对整个文件夹的fasta文件

书坊尚 2021-09-19 阅读 67
日记本

1、准备fasta文件,这里使用的是13份野生稻

2 、准备日本晴基因组文件

3、bwa 生成索引文件  

bwa index /public/home/fengting/database/rice/ribenqing/IRGSP1.0_genome.fasta

4、生成sam文件:

for i in $(ls /public/home/fengting/task/4.25/raw_wide);

do bwa mem -t 10 /public/home/fengting/database/rice/ribenqing/IRGSP-1.0_genome.fasta /public/home/fengting/task/4.25/raw_wide/$i > /public/home/fengting/task/4.25/raw_bam/$i.sam;

done

举报

相关推荐

0 条评论