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的文件名称. 搞定!
最新文章
- ccna第四学期final一二三套题
- FreeSWITCH IVR中lua调用并执行nodejs代码
- 3399 cpu绑核
- 159. Leetcode 122. 买卖股票的最佳时机 II (贪心算法-股票题目)
- 科普文丨Java 单元测试相关概念理清
- word中图片超出页边距_如何在Word中更改页边距
- 每日英语:The Invasion of the Online Tutors
- 数据库基础知识——参考数据库基本概念6版
- 深度学习-数学-第一篇-标量,向量,矩阵,张量
- RTP包里面得到H.264数据和AAC数据的方法
- 基于FPGA VHDL 的 ASK调制与解调设计(附源码)
- 最短路 POJ2387
- HTML5设置直线长度,cad绘制直线怎样设置长度
- 什么是精细化营销和精准化营销 - whale 帷幄
- CSDN博客字体样式、字体大小、字体颜色、首行缩进调节
- u盘能不能给联想服务器做系统盘,u盘能当系统盘吗?怎么把U盘做成系统盘
- 【SNS】什么是精英?什么是草根?为什么要SNS?
- java通信:远程画板
- Java集合中的fail-fast快速失败机制
- 我在上海乐字节学习java的第二十五天(持续更新中)
热门文章
- 【修炼五】用户需求系统需求
- Tomcat修改默认主页
- 公司要抽奖活动?50行Python代码制作了一个转盘抽奖小程序
- opencv3中的图像混合操作
- 郁金香java_郁金香搜索引擎的方案
- cobalt strick 4.0 系列教程(4)---监听器和基础设施管理
- android dialog遮挡键盘,彻底解决软键盘遮挡DialogFragment
- Execution failed for task ':app:processDebugResources'. No slave process to process jobs, aborting
- python面试题-如“上海 深圳 深圳 上海“,要求输入一个匹配模式,比如: aabb,判断是否符合
- 快手发布二季度及中期业绩:Q2收入191亿元 同比劲增48.8%