生信人迷惑的一天 bam转fq

169 阅读 0 评论 112 点赞

我是靠谱客的博主动人服饰，这篇文章主要介绍生信人迷惑的一天 bam转fq，现在分享给大家，希望可以做个参考。

忙活了很久，组装基因组我发现还是有点不对劲我决定从头开始

我的HIFI数据是bam格式，今天遇到一个行家，大佬问我的bam是subreadsbam 还是ccsbam。这个问题难到我了，从老板那里拿到数据的时候我也没看名字，文件名里面也没有写。我草草给大佬截完图之后大佬说这是CCSbam，直接转fq就行，用bamTofq。我测试了一下！

bamToFastq [OPTIONS] -i <BAM> -fq <FASTQ>

结果是

	bamtofastq v1.4.1
Unknown flag: '-i'
Usage:
bamtofastq [options] <bam> <output-path>
bamtofastq (-h | --help)

就在我十分无助的时候我找到了samtools自带的bam2fq，于是我有了如下命令。

 samtools bam2fq new.bam > name.fastq && gzip -c name.fastq >name.fastq.gz

在执行这一步之前，我依旧是对我的bam文件进行了排序操作

samtools sort -n old.bam -o new.bam

这个地方，导致我的new.bam文件比old.bam文件小了0.3G，猛然间想起大佬说CCSbam不需要排序，直接转就行…？

![加群]charles_kiko@163.com 备注基因组加群。

以上就是动人服饰最近收集整理的关于生信人迷惑的一天 bam转fq的全部内容，更多相关生信人迷惑的一天内容请搜索靠谱客的其他文章。

本图文内容来源于网友提供，作为学习参考使用，或来自网络收集整理，版权属于原作者所有。