fastq文件转化成bam文件
进行数据分析的前提就是要获取数据,我们可以从相应的网站上下载自己需要的数据。例如https://www.ebi.ac.uk/ena/browser/view/PRJNA189204?show=reads。
选择自己想要下载的文件,将其下载到本机。下载的文件都是经过压缩的,就需要先解压文件。
在这里是使用unzip命令解压文件,在Ubuntu中,可以使用sudo apt install unzip命令,安装完成之后就可以使用unzip命令了。解压完成后会得到两个压缩的fastq文件,需要用gunzip命令解压。
利用bowtie2将fastq文件转化成bam文件。Ubuntu中可以用sudo apt install bowtie2安装bowtie2.使用以下命令就可以完成转换。
bowtie2 -x index -1 SRR893061_1.fastq -2 SRR893061_2.fastq -S bowtie.sam
samtools view -bS bowtie.sam -o bowtie.bam
samtools sort bowtie.bam -o bowtie.sorted.bam
samtools index bowtie.sorted.bam
fastq文件转化成bam文件相关推荐
- 生物信息学:bedGraph文件、Bed文件转、BAM文件转化
目录 引言 bam文件 bed文件 bedGraph文件 R语言和linux实现从bedgraph到bam的转化 从bedgraph到bed genome.size.txt bed转bam 引言 Ba ...
- TensorFlow:将ckpt文件固化成pb文件
TensorFlow:将ckpt文件固化成pb文件 本文是将yolo3目标检测框架训练出来的ckpt文件固化成pb文件,主要利用了GitHub上的该项目. 为什么要最终生成pb文件呢?简单来说就是直接 ...
- 利用Python批量将csv文件转化成xml文件
文章目录 一.前言 二.Python代码实现 一.前言 将 csv 格式转换成xml格式有许多方法,可以用数据库的方式,也有许多软件可以将 csv 转换成xml.但是比较麻烦,本文利用 Python ...
- Pcap文件转化成Pcd文件
通过RSview将点云文件保存成了Pcap格式,但这种格式不能很好的支持PCL点云库,故不能很好的实现点云获取.滤波.分割.配准.检索.特征提取.识别.追踪.曲面重建.可视化等,因此我们需要将Pcap ...
- QT Designer 生成的ui文件转化成py文件以及简单使用
QT Designer 生成的ui文件转化成py文件以及简单使用 设计 转换 使用 方法一 方法二 设计 首先用QTdesigner 设计你的界面 然后保存成.ui文件 这一步大家应该都可以做到,就不 ...
- sarscape 将dem文件转化成stl_STL文件,一种前处理网格划分技术??
源:吴冠中作品 点击关注CAE仿真空间, 点亮"在看",优质内容不错过对于从事专业仿真技术的工程师而言,我们已经习惯而且能够熟练的利用诸如ANSA.Hypermesh等网格划分前处 ...
- 利用pandas实现json文件转化成csv文件
补充上篇博客提到的json文件数据转换成csv文件.作为数据分析最常用文件类型json与csv,出于一定情况下,我们需要将json中字典类型的数据,转换为csv存储,这又用到python强大的pand ...
- Vaa3D_批量将tiff文件转化成v3draw文件
通过这个代码可以实现输入tif.tiff文件所在的文件夹,在对应的位置生成一个文件夹然后在该文件夹中生成v3draw图像.如果不需要可以注释这部分代码. QString folder=QFileDia ...
- 一份basic文件转化成c文件(自己收藏)
basic的原文件,先通过BCX转化器,转化一下,注意转化过来的东东,错误是N的多哦. 得动手自己修改: 10 REM ************************************* ...
- ipynb文件转化成py文件
在当前文件夹运行cmd, 输入 jupyter nbconvert --to script xxx.ipynb [注]xxx.ipynb是需要转换成py的文件名称. 搞定!
最新文章
- 用 Pyecharts 制作炫酷的可视化大屏
- 【学习笔记】高斯整数、高斯素数、费马平方和(全部相关概念及例题详解)《初等数论及其应用》
- linux 轻量化图形界面,用这13个方法,帮你做出真正轻量化的移动 App 设计
- 打包caddy为docker镜像
- C++用string 定义字符串数组
- Scala数组:使用()代替java的[]
- 客座编辑:杜小勇(1963-),男,博士,中国人民大学信息学院教授、博士生导师。...
- 单词查找树(信息学奥赛一本通-T1337)
- Elementui 自定义loading
- 使用动态SQL获取具有非空字段的自定义表结构
- Ubuntu 16.04安装Guake Terminal终端(使用一键唤醒功能)
- BZOJ4435 : [Cerc2015]Juice Junctions
- spring 的web.xml的加载顺序
- Java实现文件及文件夹的删除
- 基于PHP+MySQL游戏商城销售网站的设计与实现
- linux下 VSCode快捷键
- 8uftp,怎么使用8uftp
- Android 开启热点 适配8.0以上
- 人工智能-语音交互-NLP自然语言(三) 语法-敲黑板
- 2022年中国智能家居产业链图谱 | 产业链全景图