首页 大数据 基因数据处理20之BWASW算法ref分块建立索引然后比对

基因数据处理20之BWASW算法ref分块建立索引然后比对

SRR003161h20.fastq <span style="font-family: Arial,Helvetica,sans-serif; font-size: 12px; background-color: rgb(255,255,255);">三次测试:</span>

hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub$ bwa bwasw GRCH38chr1L3556522.fna SRR003161h20.fastq >SRR003161h20BWASWchr1.sam
[M::bwa_idx_load_from_disk] read 0 ALT contigs
[bsw2_aln] read 5 sequences/pairs (2691 bp) ...
[main] Version: 0.7.13-r1126
[main] CMD: bwa bwasw GRCH38chr1L3556522.fna SRR003161h20.fastq
[main] Real time: 0.315 sec; CPU: 0.257 sec
hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub$ bwa bwasw GRCH38chr1L3556522.fna SRR003161h20.fastq >SRR003161h20BWASWchr1.sam
[M::bwa_idx_load_from_disk] read 0 ALT contigs
[bsw2_aln] read 5 sequences/pairs (2691 bp) ...
[main] Version: 0.7.13-r1126
[main] CMD: bwa bwasw GRCH38chr1L3556522.fna SRR003161h20.fastq
[main] Real time: 0.306 sec; CPU: 0.251 sec
hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub$ bwa bwasw GRCH38chr1L3556522.fna SRR003161h20.fastq >SRR003161h20BWASWchr1.sam
[M::bwa_idx_load_from_disk] read 0 ALT contigs
[bsw2_aln] read 5 sequences/pairs (2691 bp) ...
[main] Version: 0.7.13-r1126
[main] CMD: bwa bwasw GRCH38chr1L3556522.fna SRR003161h20.fastq
[main] Real time: 0.312 sec; CPU: 0.253 sec
hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub$ bwa bwasw GRCH38chr2L3459909.fna SRR003161h20.fastq >SRR003161h20BWASWchr2.sam
[M::bwa_idx_load_from_disk] read 0 ALT contigs
[bsw2_aln] read 5 sequences/pairs (2691 bp) ...
[main] Version: 0.7.13-r1126
[main] CMD: bwa bwasw GRCH38chr2L3459909.fna SRR003161h20.fastq
[main] Real time: 0.293 sec; CPU: 0.239 sec
hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub$ bwa bwasw GRCH38chr2L3459909.fna SRR003161h20.fastq >SRR003161h20BWASWchr2.sam
[M::bwa_idx_load_from_disk] read 0 ALT contigs
[bsw2_aln] read 5 sequences/pairs (2691 bp) ...
[main] Version: 0.7.13-r1126
[main] CMD: bwa bwasw GRCH38chr2L3459909.fna SRR003161h20.fastq
[main] Real time: 0.287 sec; CPU: 0.238 sec
hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub$ bwa bwasw GRCH38chr2L3459909.fna SRR003161h20.fastq >SRR003161h20BWASWchr2.sam
[M::bwa_idx_load_from_disk] read 0 ALT contigs
[bsw2_aln] read 5 sequences/pairs (2691 bp) ...
[main] Version: 0.7.13-r1126
[main] CMD: bwa bwasw GRCH38chr2L3459909.fna SRR003161h20.fastq
[main] Real time: 0.293 sec; CPU: 0.240 sec
hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub$ bwa bwasw GRCH38chr1chr2L7016431.fna SRR003161h20.fastq >SRR003161h20BWASWchr1chr2.sam
[M::bwa_idx_load_from_disk] read 0 ALT contigs
[bsw2_aln] read 5 sequences/pairs (2691 bp) ...
[main] Version: 0.7.13-r1126
[main] CMD: bwa bwasw GRCH38chr1chr2L7016431.fna SRR003161h20.fastq
[main] Real time: 0.693 sec; CPU: 0.598 sec
hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub$ bwa bwasw GRCH38chr1chr2L7016431.fna SRR003161h20.fastq >SRR003161h20BWASWchr1chr2.sam
[M::bwa_idx_load_from_disk] read 0 ALT contigs
[bsw2_aln] read 5 sequences/pairs (2691 bp) ...
[main] Version: 0.7.13-r1126
[main] CMD: bwa bwasw GRCH38chr1chr2L7016431.fna SRR003161h20.fastq
[main] Real time: 0.646 sec; CPU: 0.588 sec
hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub$ bwa bwasw GRCH38chr1chr2L7016431.fna SRR003161h20.fastq >SRR003161h20BWASWchr1chr2.sam
[M::bwa_idx_load_from_disk] read 0 ALT contigs
[bsw2_aln] read 5 sequences/pairs (2691 bp) ...
[main] Version: 0.7.13-r1126
[main] CMD: bwa bwasw GRCH38chr1chr2L7016431.fna SRR003161h20.fastq
[main] Real time: 0.639 sec; CPU: 0.582 sec

只对两个串比较,时间较短,差不多查找比对时间可以减半

关于作者: dawei

【声明】:乐山站长网内容转载自互联网,其相关言论仅代表作者个人观点绝非权威,不代表本站立场。如您发现内容存在版权问题,请提交相关链接至邮箱:bqsm@foxmail.com,我们将及时予以处理。

热门文章