原始 BAM 文件和 sort 之后 BAM 文件的行数，是一样的。
SEQanswers：BAM is compressed. Sorting helps to give a better compression ratio because similar sequences are grouped together.

bam转回fq时报错: github查找同问题结果1 2 3

 *****WARNING: Query 17 is marked as paired, but its mate does not occur next to it in your BAM file.  Skipping.
*****WARNING: Query 13 is marked as paired, but its mate does not occur next to it in your BAM file.  Skipping.
*****WARNING: Query 223 is marked as paired, but its mate does not occur next to it in your BAM file.  Skipping.

sort -n 后 warning行数变少1,333,109,095变为11,212,985
？？？：

nohup samtools sort -n S1_T_SRR1273943.bam -o ./S1_T_SRR1273943.sortedByName.bam >log.S1.bam.sortbyname 2>&1 &nohup bedtools bamtofastq  -i S1_T_SRR1273943.sortedByName.bam -fq S1_T_1.fq -fq2 S1_T_2.fq > log_S1_sortedbyname 2>&1 &

用sortbam得到的fq后续分析继续报错

reads: 0 |ERROR: The mate1 read name did not match the mate2 read name. Resynchr onization support needs to be implemented.

？暂时解决方法，提取序列之后可以按照read name排序，然后提取。为何对于排序后的用bedtools bamtofastq得到的结果会后续报错，而samtools fastq暂时没有

samtools sort -n  bam -o sorted.bam| samtools fastq -1 read_1.fq -2 read_2.fq -s singleton.fq -

samtools fastq

一般而言BAM文件都是按照位置信息排序，想要找到配对的reads，要么是根据read的编号进行排序（这个方法要求额外的内存和存储空间），或者就是在提取的时候记录当前的read的ID，再找到另一端ID后释放内存空间。

BAM中reads名称和fq中reads名称差异，mate1 read name----mate2 read name？

bam文件处理转fq相关推荐

bam文件转fq.gz文件
bam转fq过程 1.fq文件格式 2.bam文件格式 3.转换思路 3.1 软件bedtools自带功能 3.2 自己写代码 3.3 代码示例 4.参考资料 1.fq文件格式 fastq格式是一 ...
如何高效地从BAM文件中提取fastq
在一年前,我写过一篇文章,叫做如何从BAM文件中提取fastq,之前也发现了从BAM里面提取Fastq是有些麻烦,只不过最后通过samtools的子命令实现了数据提取,实现功能之后也没有再去思考如何提 ...
bam获取序列_如何从BAM文件中提取fastq
虽然高通量测序分析最常用的操作是将fastq比对到参考基因组得到BAM文件,但偶尔我们也需要提取BAM文件中特定区域中fastq.最开始我认为这是一个非常简单的操作,因为samtools其实已经提供了 ...
bam获取序列_如何高效地从BAM文件中提取fastq
在一年前,我写过一篇文章,叫做如何从BAM文件中提取fastq,之前也发现了从BAM里面提取Fastq是有些麻烦,只不过最后通过samtools的子命令实现了数据提取,实现功能之后也没有再去思考如何提 ...
linux怎么查看一个bam文件,生信分析过程中这些常见文件的格式以及查看方式你都知道吗？...
原标题:生信分析过程中这些常见文件的格式以及查看方式你都知道吗? 生信分析过程中,会与很多不同格式的文件打交道,除了原始测序数据 fastq 之外,还需要准备基因组文件 fasta 格式和基因注释文件 ...
Samtools应用指南-处理Sam与Bam文件
安装去官网下载想要的版本 tar jxvf samtools-1.9.tar.bz2 cd samtools-1.9 ./configure --prefix=全路径/samtools-1.9 ma ...
samtools 检测bam文件的完整度
检测bam文件的完整度 samtools view T_recal.bam|head samtools view T_recal.bam|tail for i in *.bam ;do (samtoo ...
linux bam文件格式介绍,Sam和bam文件说明
SAM文件 SAM(Sequence Alignment/Map)格式是一种通用的比对格式,用来存储reads到参考序列的比对信息. SAM是一种序列比对格式标准,由sanger制定,是以TAB为分割 ...
bam文件读取_bam格式文件处理大全(一）
sam文件是短序列比对生成的文件,是二代测序中最核心的文件.在RNAseq,变异检测等分析中,都需要首先生成sam文件格式.bam文件是sam格式的二进制格式,转换为二进制之后,可以减小文件的存储.掌 ...

bam文件处理转fq

BAM中reads名称和fq中reads名称差异，mate1 read name----mate2 read name？

bam文件处理转fq相关推荐

最新文章

热门文章

bam文件处理 转fq

BAM中reads名称和fq中reads名称差异，mate1 read name----mate2 read name？

bam文件处理 转fq相关推荐

最新文章

热门文章

bam文件处理转fq

bam文件处理转fq相关推荐