进行数据分析的前提就是要获取数据,我们可以从相应的网站上下载自己需要的数据。例如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文件相关推荐

  1. 生物信息学:bedGraph文件、Bed文件转、BAM文件转化

    目录 引言 bam文件 bed文件 bedGraph文件 R语言和linux实现从bedgraph到bam的转化 从bedgraph到bed genome.size.txt bed转bam 引言 Ba ...

  2. TensorFlow:将ckpt文件固化成pb文件

    TensorFlow:将ckpt文件固化成pb文件 本文是将yolo3目标检测框架训练出来的ckpt文件固化成pb文件,主要利用了GitHub上的该项目. 为什么要最终生成pb文件呢?简单来说就是直接 ...

  3. 利用Python批量将csv文件转化成xml文件

    文章目录 一.前言 二.Python代码实现 一.前言 将 csv 格式转换成xml格式有许多方法,可以用数据库的方式,也有许多软件可以将 csv 转换成xml.但是比较麻烦,本文利用 Python ...

  4. Pcap文件转化成Pcd文件

    通过RSview将点云文件保存成了Pcap格式,但这种格式不能很好的支持PCL点云库,故不能很好的实现点云获取.滤波.分割.配准.检索.特征提取.识别.追踪.曲面重建.可视化等,因此我们需要将Pcap ...

  5. QT Designer 生成的ui文件转化成py文件以及简单使用

    QT Designer 生成的ui文件转化成py文件以及简单使用 设计 转换 使用 方法一 方法二 设计 首先用QTdesigner 设计你的界面 然后保存成.ui文件 这一步大家应该都可以做到,就不 ...

  6. sarscape 将dem文件转化成stl_STL文件,一种前处理网格划分技术??

    源:吴冠中作品 点击关注CAE仿真空间, 点亮"在看",优质内容不错过对于从事专业仿真技术的工程师而言,我们已经习惯而且能够熟练的利用诸如ANSA.Hypermesh等网格划分前处 ...

  7. 利用pandas实现json文件转化成csv文件

    补充上篇博客提到的json文件数据转换成csv文件.作为数据分析最常用文件类型json与csv,出于一定情况下,我们需要将json中字典类型的数据,转换为csv存储,这又用到python强大的pand ...

  8. Vaa3D_批量将tiff文件转化成v3draw文件

    通过这个代码可以实现输入tif.tiff文件所在的文件夹,在对应的位置生成一个文件夹然后在该文件夹中生成v3draw图像.如果不需要可以注释这部分代码. QString folder=QFileDia ...

  9. 一份basic文件转化成c文件(自己收藏)

    basic的原文件,先通过BCX转化器,转化一下,注意转化过来的东东,错误是N的多哦. 得动手自己修改: 10   REM  ************************************* ...

  10. ipynb文件转化成py文件

    在当前文件夹运行cmd, 输入 jupyter nbconvert --to script xxx.ipynb [注]xxx.ipynb是需要转换成py的文件名称. 搞定!

最新文章

  1. ccna第四学期final一二三套题
  2. FreeSWITCH IVR中lua调用并执行nodejs代码
  3. 3399 cpu绑核
  4. 159. Leetcode 122. 买卖股票的最佳时机 II (贪心算法-股票题目)
  5. 科普文丨Java 单元测试相关概念理清
  6. word中图片超出页边距_如何在Word中更改页边距
  7. 每日英语:The Invasion of the Online Tutors
  8. 数据库基础知识——参考数据库基本概念6版
  9. 深度学习-数学-第一篇-标量,向量,矩阵,张量
  10. RTP包里面得到H.264数据和AAC数据的方法
  11. 基于FPGA VHDL 的 ASK调制与解调设计(附源码)
  12. 最短路 POJ2387
  13. HTML5设置直线长度,cad绘制直线怎样设置长度
  14. 什么是精细化营销和精准化营销 - whale 帷幄
  15. CSDN博客字体样式、字体大小、字体颜色、首行缩进调节
  16. u盘能不能给联想服务器做系统盘,u盘能当系统盘吗?怎么把U盘做成系统盘
  17. 【SNS】什么是精英?什么是草根?为什么要SNS?
  18. java通信:远程画板
  19. Java集合中的fail-fast快速失败机制
  20. 我在上海乐字节学习java的第二十五天(持续更新中)

热门文章

  1. 【修炼五】用户需求系统需求
  2. Tomcat修改默认主页
  3. 公司要抽奖活动?50行Python代码制作了一个转盘抽奖小程序
  4. opencv3中的图像混合操作
  5. 郁金香java_郁金香搜索引擎的方案
  6. cobalt strick 4.0 系列教程(4)---监听器和基础设施管理
  7. android dialog遮挡键盘,彻底解决软键盘遮挡DialogFragment
  8. Execution failed for task ':app:processDebugResources'. No slave process to process jobs, aborting
  9. python面试题-如“上海 深圳 深圳 上海“,要求输入一个匹配模式,比如: aabb,判断是否符合
  10. 快手发布二季度及中期业绩:Q2收入191亿元 同比劲增48.8%