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

最新推荐文章于 2023-07-05 11:57:19 发布
最新推荐文章于 2023-07-05 11:57:19 发布
阅读量1.8k 收藏
点赞数
分类专栏： 基因数据处理 文章标签： 基因数据处理20之BWASW算法ref分
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xubo245/article/details/51221035
版权 


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

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

KeepLearningBigData
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
sam to bam文件出错
fethin的博客
08-18 8824
开始学习生物信息分析，在ubuntu上运行了samtools view，结果显示错误，如下： # samtools view -@ 60 -bS *.sam > ra24p.bam [samopen] no @SQ lines in the header. 尝试改变参数也不行 就查看了sam文件 #head -10 ra24p.sam [M::bwa_idx_load_
BWA MEM比对人类参考基因组详解
热门推荐
睿智如水
11-10 3万+
在获得下机数据后，做的第一步是质控。质控工具有很多，这里就不做一一介绍了。这里讲如何使用BWA MEM将质控合格的数据比对到参考基因组上。         BWA是一款基于BWT的快速比对工具，其由三个算法组成。这三个算法分别是：BWA backtrack, BWA SW and BWA MEM。其中，BWA MEM是最新的，其更快更准确，更适合用于人重数据分析。对于上述三种算法，首先需要使用索
Samtools view转sam为bam报错[main_samview] truncated file.
jiangyi_pig的博客
11-20 6573
报错分析：在shell里加nohup &bwa比对出来的sam文件里面有除注释行和比对结果以外的不明信息，转bam格式时会报错 [W::sam_read1] parse error at line 1 [main_samview] truncated file. head查看前三行内容为 [M::bwa_idx_load_from_disk] read 0 ALT contigs @SQSN:NC_050096.1LN:308452471 @SQSN:NC_0...
bwa mem比对结果错误，sam文件不能被samtools识别的原因之一
fethin的博客
09-21 7454
解决[M::bwa_idx_load_from_disk] read 0 ALT contigs的问题 总结如下，bwa mem比对结果错误，sam文件不能被samtools识别的原因之一是bwa安装的问题！ 写这个初级的帖子，为后来人遇到同样问题的人，在百度搜索的时候能够找到能解决问题的帖子！
2023.07.05【bug笔记】|bwa比对报错：Segmentation fault
最新发布
07-05 368
Segmentation fault 这个报错还有一种情况会与参考基因组有联系，即参考基因组过大时，由于内存不足导致溢出。有更多问题可以加VX：bbplayer2021 （木青）进群交流，备注 申请加入生信交流群。
BWA比对及Samtools提取目标序列
weixin_56827666的博客
05-13 4620
今天想看一下自己的序列里面会不会有某细菌基因组存在，主要使用BWA和Samtools：bwa主要用于将低差异度的短序列与参考基因组进行比对。主要包含三种比对算法：backtrack、SW和MEM，第一种只支持短序列比对（<100bp），后两种支持长序列比对(70bp~1M)，并支持分割比对（split alignment）。MEM算法是最新的也是官方推荐的。BWA-MEM 是一种新的比对算法，用于将测序 reads 或者组装后 contigs 比对至大型参考基因组，例如人参考基因组。它会自动选择局部比
BWA-MEM使用中遇到错误fail to open “**.fastq“
m0_63679102的博客
09-30 1349
bwa mem 出现 fail to open file `***.fastq'.
基因数据处理21之BWASW算法ref分块建立索引然后比对（ref切分为四段，read为250条）
Keep Learning
04-22 1553
1.时间分析 对ref为单条染色体进行比对，第一次比对在3-5s不等，对chr1-4比对，在20s左右 连续比对多次后，对单染色体比对降到1s左右，chr1-4降到2s左右 不懂为什么比一次比对时间比较长，后面几次比对时间变短 运行代码： hadoop@Master:~/cloud/adam/xubo/data/GRCH38Sub$ bwa bwasw GRCH38chr1L35
基因数据处理24之BWASW算法ref分块建立索引然后比对（ref切分为四段，read为25000条，4线程）
Keep Learning
04-22 2281
比默认的要块1倍左右 1 [M::bwa_idx_load_from_disk] read 0 ALT contigs [bsw2_aln] read 25000 sequences/pairs (14470082 bp) ... [main] Version: 0.7.13-r1126 [main] CMD: bwa bwasw -t 4 GRCH38chr1L3556522.fna SRR0
BWT完整算法
06-03
BWT的完整算法，包括SA，Occ等数组的建立。用于在基因链中快速匹配基因。
基因数据处理64之bwamem处理500bp和1000bp的记录
Keep Learning
12-20 1501
xubo@xubo:~/xubo/data/alignment/cs-bwamem$ bwa mem bwaindex/GRCH38chr1L3556522.fasta g38l500N10000.fq >g38l500N10000.bwamem.sam [M::bwa_idx_load_from_disk] read 0 ALT contigs [M::process] read 10000 se
基因数据处理23之BWASW算法ref分块建立索引然后比对（ref切分为四段，read为25000条）
Keep Learning
04-22 1517
结论： read为25000条对单条染色体进行比对需要110秒左右，对chr1-4比对需要160秒左右 运行结果： 1 [M::bwa_idx_load_from_disk] read 0 ALT contigs [bsw2_aln] read 17040 sequences/pairs (10000385 bp) ... [bsw2_aln] read 7960 sequences/pai
基因数据处理14之BWA三种方式bwa、BWA_SW、BWA_MEM使用
Keep Learning
03-13 1万+
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.
基因数据处理115之重新运行SparkBWA本地模式
Keep Learning
06-03 993
更多代码请见：https://github.com/xubo245 基因数据处理系列 1.记录 hadoop@Master:~/disk2/xubo/tools/SparkBWA/build$ ./pairedGRCH38L1Local.sh [Stage 3:&gt; (...
BWA-MEM算法结构分析
轩令的技术博客
07-21 1万+
一、BWA-MEM函数框架 1 读入 bwt、options、reads； 2 利用mem_chain生成chain； 3 利用mem_chain_flt过滤掉部分chain； 4 利用mem_chain2aln生成比对结果元数据。 1.第一步：数据输入 加载已经生成的bwt表。接口的参数文件名为：xx.fasta；实际中包含具有以下几个后缀名的文件.amb,.ann,.bw...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值