基因数据处理14之BWA三种方式bwa、BWA_SW、BWA_MEM使用

发布时间：2020-12-14 01:59:52 所属栏目：大数据来源：网络整理

导读：1.构建索引： bwa index ref.fa 或者从ftp下载，请参考【1】 2.BWA： bwa.sh为脚本文件 hadoop@Master:~/cloud/adam/xubo/data/test20160310/bwa3sh$ cat bwa.sh #!/bin/bashbwa aln ../GCA_000001405.15_GRCh38/GCA_000001405.15_GRCh38_full_analysis_set.

1.构建索引：

bwa index ref.fa

或者从ftp下载，请参考【1】

2.BWA：

bwa.sh为脚本文件

hadoop@Master:~/cloud/adam/xubo/data/test20160310/bwa3sh$ cat bwa.sh 
#!/bin/bash
bwa aln ../GCA_000001405.15_GRCh38/GCA_000001405.15_GRCh38_full_analysis_set.fna ../SRR003161h20.fastq > SRR003161h20t1.sai
bwa samse ../GCA_000001405.15_GRCh38/GCA_000001405.15_GRCh38_full_analysis_set.fna SRR003161h20t1.sai ../SRR003161h20.fastq >SRR003161h20t1.sam

文件的获取请参考【1】和【2】

3.BWA_SW:

bwasw.sh也是脚本文件

hadoop@Master:~/cloud/adam/xubo/data/test20160310/bwasw$ cat bwasw.sh 
#!/bin/bash
bwa bwasw ../GCA_000001405.15_GRCh38/GCA_000001405.15_GRCh38_full_analysis_set.fna ../SRR003161h20.fastq >SRR003161h20t1.sam

4.BWA_MEM:

bwa mem GCA_000001405.15_GRCh38/GCA_000001405.15_GRCh38_full_analysis_set.fna SRR003161h20.fastq >SRR003161h20.sam

具体的格式请参考【3】

参考：

【1】http://blog.csdn.net/xubo245/article/details/50858334

【2】?http://blog.csdn.net/xubo245/article/details/50879044

【3】：

（编辑：李大同）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!