生信人迷惑的一天 bam转fq

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

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

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

结果是

 bamtofastq v1.4.1Unknown 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相关推荐

生信人的自我修养：Linux 命令速查手册
标题:生信人的自我修养:Linux 命令速查手册目标:致力于为生信人打造一个完整的 Linux 命令速查手册作者:简佐义(jianzuoyi@qq.com) 版本:1.0 日期:2020-11-2 ...
2019linux考试题库,生信人的linux考试 2019-05-03
生信人的linux考试(20题) 1.在任意文件夹下面创建形如 1/2/3/4/5/6/7/8/9 格式的文件夹系列 mkdir 生信人的linux考试 ls cd 生信人的linux考试/ mkdi ...
生信人写程序2. Editplus添加Perl, Shell, R模板和语法高亮
https://www.editplus.com/ 前言 "工欲善其事必先利其器",生信工程师每天写代码.搭流程,而且要使用至少三门编程语言,没有个好集成开发环境(IDE,Inte ...
desktop docker 无法卸载_docker，生信人的福音！
docker的原理是什么?我也不知道呀!!! 我只是知道这个东西比虚拟机,conda什么的都要好用. 为什么? 因为它解决了一个对我们这些生信小白和编程小白来说最致命的问题,那就是软件安装的问题(软件 ...
引用另一模板的宏_生信人值得拥有的编程模板Shell
前言 "工欲善其事必先利其器",生信工程师每天写代码.搭流程,而且要使用至少三门编程语言,没有个好集成开发环境(IDE,Integrated Development Environ ...
生存分析系列教程（一）使用生信人工具盒进行生存分析
生信人工具盒是生信人团队的开发的一款软件,非常方便.下面我将演示一下如何通过这款软件进行生存分析.为了方便大家理解,形式依然是数据结构-操作-结果解读. 1. 表达矩阵与生存信息矩阵表达矩阵依然 ...
生信c语言,生信人的R使用
接下来介绍R语言: [生信技能树]生信人应该这样学R语言 R语言在你开始R之旅前,建议你看看下面这两个 1. 介绍R语言及Rstudio 了解R,Rstudio及R包;安装的包在packages中检 ...
生信人的20个R语言习题的答案
这是生信技能树关于生信人的20个R语言习题的答案: 1 安装R包数据包: ALL, CLL, pasilla, airway 软件包:limma,DESeq2,clusterProfiler 工具包 ...
python perl 比较生信_科学网—生信人写程序1. Perl语言模板及配置 - 刘永鑫的博文...
科学网对Markdown排版支持较差,对格式不满意的用户请跳转至 CSDN 或微信阅读: 如果感觉文章对您有帮助,想继续阅读同类文章,请扫描下方二维码关注"生信宝典"公众号,每天接 ...

生信人迷惑的一天 bam转fq

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

生信人迷惑的一天 bam转fq相关推荐

最新文章

热门文章