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】:
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。