bwa mem 报错处理：[mem_sam_pe] paired reads have different names

最新推荐文章于 2024-07-26 17:36:26 发布
最新推荐文章于 2024-07-26 17:36:26 发布
阅读量2k 收藏
点赞数
文章标签： python
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/sinat_40038704/article/details/120973808
版权

背景：

从samtools sort 默认排序后的bam文件中提取fastq序列并对其格式化，对格式化后的fastq文件重新比对到参考基因组，报错如下：”[mem_sam_pe] paired reads have different names: “A00575:297:HWHKYDMXX:1:1331:22372:31814:”, “A00575:297:HWHKYDMXX:2:2167:15004:31297:”
在这里插入图片描述

原因分析：

最开始在网上也查找了很多的提问，有说使用BBtools工具处理，最终看到有一个解答更为简单解决该问题（如何从BAM文件中提取fastq by徐洲更），比较原始fastq的两个文件，R1和R2reads name排序是统一的，最终确认为是由于reads name不统一导致该问题。

解决问题

重新使用samtools sort -n 对已排序的bam按照reads name 重新排序，重新进行fastq提取
在这里插入图片描述

UP小吴同学
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
BWA mem序列比对时出现：paired reads have different names 问题解决
weixin_47621318的博客
10-10 2194
BWA mem序列比对时出现：paired reads have different names 问题解决
bwa-mem2:下一个版本的bwa-mem
05-03
输出应与原始的bwa-mem版本0.7.17相匹配。 从提交e0ac59e开始，我们有一个git子模块safestringlib。 要获得它，请在克隆时使用--recursive或在已克隆的存储库中使用“ git submodule init”和“ git submodule upd
mapping生成sam文件时出现[mem_sam_pe] paired reads have different names错误
weixin_33744854的博客
06-20 1194
用以下命令修复： bbrename.sh in1=read1.fq in2=read2.fq out1=renamed1.fq out2=renamed2.fq 　　bbrename.sh 下载地址网上自行搜索 对于多个fq文件，可以用以下命令： while read line do nohup bbrename.sh in=${line}_R1.fq i...
BWA-MEM使用中遇到错误fail to open “**.fastq“
m0_63679102的博客
09-30 1548
bwa mem 出现 fail to open file `***.fastq'.
基因数据处理100之bwamem算法处理100万条paired-reads数据GRCH38chr1L3556522N1000000L100paired12
Keep Learning
01-11 1604
运行记录：hadoop@Master:~/xubo/ref/GRCH38L1Index/pe$ bwa mem ../GRCH38chr1L3556522.fasta GRCH38chr1L3556522N1000000L100paired1.fastq GRCH38chr1L3556522N1000000L100paired2.fastq >GRCH38chr1L3556522N1000000L1
基因数据处理6之BWA_MEM无法分配内存
Keep Learning
03-12 3891
基因数据处理之BWA_MEM无法分配内存：      建立BWA索引的时候内存不足，现在用BWA-MEM又内存不足，真耗内存 hadoop@Mcnode1:~/cloud/adam/xubo/data/data_HDFS/GRCH38/GCA_000001405.15_GRCh38/test20160310$ bwa mem GCA_000001405.15_GRCh38/GCA_
修改 bam 文件中染色体名
liangbilin的博客
10-09 5096
问题描述 在分析的过程中，有些数据的染色体命名为“chr1、chr2、…、chrX、chrY”，而有些数据的染色体命名则为“1、2、…、X、Y” （也就是不包含 chr 字符）。这里，通过代码对 bam 文件作为修改，实现染色体名的统一。 代码实现 假设我们有一个名为 test.bam 的文件，其中染色体名不包含chr字符，需要在染色体名前加上chr字符。 通过 samtools 和 shell 实现 （注：samtools reheader 需要给一个- 的参数，不给会报错）： samtools vi
测序数据量那些事
S_AGZX的博客
03-26 8636
reads、fastq、fastq.gz大小关系
SAM文件统计python最终版1.zip_DNA_dna sequencing_python SAM评价_sam
09-23
标题中的“SAM文件统计python最终版1.zip”表明这是一个使用Python编程语言处理SAM（Sequence Alignment/Map）格式文件的项目，可能是一个统计分析工具或脚本的最终版本。SAM是生物信息学领域广泛使用的文件格式，...
docker_bwa_aligner:展示概念验证 bwa 对齐工作流的 Docker 容器
07-01
docker_bwa_aligner 这是一个基于 docker 容器的系统的概念证明，用于运行简单的序列对齐工作流程。 该系统的基本部分是： 将目录 /home/bwa_user/data 导出到其他容器的数据卷。 将目录 /home/bwa_user/bwa_...
bwa-0.7.17_fightingm8b_ubuntu_
10-03
BWA提供了多种比对模式，如"BWA-MEM"（适用于长读段）和"BWA-PE"（适用于双端测序数据）。"bwamem.c"和"bwape.c"是这些模式的具体实现。此外，"ksw.c"实现了KSW2算法，用于计算比对的Smith-Waterman动态规划得分，...
samtools_primer:SAM工具
07-06
SAMtools则是一个强大的开源软件包，由Samtools项目开发，它提供了处理、分析和理解这些SAM格式文件的能力。这个工具集包括了对SAM/BAM（Compressed SAM）文件进行查看、排序、索引、统计和筛选等多种功能，对于基因...
实验记录 | scATAC-seq数据的比对（一）
weixin_40640700的博客
08-20 2892
ATAC-seq数据，用于测序的reads如图所示。 按照它的参考文件中的描述，一共有4条测序的reads。 （1）双端测序的Read 1N和Read 2N,这部分序列就是插入的序列。==>分别对应的是R1和R3的信息 （2）index序列：8bp ==>对应的是I1的信息 （3）10X barcode & spacer：16bp + 8bp =24bp ==>对应的是R2的信息 我们从序列的角度看看是怎样的情况： (base) [xxzhang@mu02 V1C]$ head..
2022.11.21【bug笔记】|bam文件报错：Cannot add sequence that already exists in SAMSequenceDictionary
11-21 636
sam文件是通过hisat2，bowtie2或者bwa将rawdata进行比对后得到的包含比对信息的数据格式。经过samtools处理后得到的bam文件经常用于后续分析，比如RNA-seq分析时，可以统计序列的插入片段也可以做后续定量，WGS流程里比对后生产的bam文件也可以去冗余获取snp位点。
BWA用法笔记
小静的博客
08-31 3303
bwa（全称Burrows-Wheeler Aligner），是用来对二代测序数据进行比对的常用工具
如何从BAM文件中提取fastq
xuzhougeng blog
02-12 2408
虽然高通量测序分析最常用的操作是将fastq比对到参考基因组得到BAM文件，但偶尔我们也需要提取BAM文件中特定区域中fastq。最开始我认为这是一个非常简单的操作，因为samtools其实已经提供了相应的工具samtools fastq. 以biostar handbook的Ebola病毒数据为例，首先获取比对得到的BAM文件。 # 建立文件夹 mkdir -p refs # 根...
全网最详细Gradio教程系列5——Gradio Client: python
最新发布
shao918516的博客
07-26 774
程序部署完成后，如何将Gradio App作为API访问使用呢，这就用到Gradio Client。本章讲解Gradio Client的三种使用方式：python、javascript和curl，受字数限制，所以分三篇博客发布。 使用Gradio Python Client非常易于将Gradio应用程序作为API使用，本节讲述gradio_client安装、如何连接Gradio应用程序、查看可用API及其使用方式、job及session等用法。 通过Gradio Python Cli
$BWA mem -t $THREADS -R "@RG\tID:re\tSM:re\tLB:lib1" ref-100-2.fa ref_1.fq > ref.sam
03-31
This command uses the BWA mem algorithm to align the reads in "ref_1.fq" to the reference genome "ref-100-2.fa". The "-t $THREADS" option specifies the number of threads to use for parallel processing. The "-R" option adds a read group header line to the output SAM file, which includes information about the sequencing library and sample. In this case, the read group ID is "re", the sample name is "re", and the library name is "lib1". The output is directed to the file "ref.sam".
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值