简介

Bowtie 2是一种用于长参考序列对齐测序读取的工具。它尤其擅长将大约50-100多个字符的比对到相对较长的(例如哺乳动物)基因组。Bowtie 2支持gapped, local,and paired-end alignment模式。可以同时使用多个处理器以实现更高的对齐速度。

使用方法

bowtie2[options]* -x <bt2-idx> {-1 <m1> -2 <m2> | -U <r> |--interleaved <i> | --sra-acc <acc> | b <bam>} -S[<sam>]

-x <bt2-idx> 参考基因组索引的名称。注意,这个名称不包括尾端的的.1.bt2 / .rev.1.bt2 / 等。bowtie2 首先在当前目录中查找指定的索引,如果没有找到的话,再在BOWTIE2_INDEXES 环境变量中查找。

-1 <m1> Comma-separated list of files containing mate 1s,例如 -1 flyA_1.fq,flyB_1.fq。 使用此选项指定的序列必须与 <m2> 中指定的文件对文件和读取对读取相对应。

-2 <m2> Comma-separated list of files containing mate 2s,例如 -2 flyA_2.fq,flyB_2.fq。 使用此选项指定的序列必须与 <m1> 中指定的文件对文件和读取对读取相对应。

-U <r> Comma-separated list of files containingunpaired reads to be aligned,例如lane1.fq,lane2.fq,lane3.fq,lane4.fq。

--interleaved读取interleaved的FASTQ 文件,其中前两条记录(8 行)代表配对。

输出:

--un <path> 读出没有的匹配上的 reads,并把它写入path中的文件

--al <path> 读出至少匹配一次的 reads,并把它写入path中的文件

举例

如果我有一些序列,想要去除掉这些序列中的污染

我把我的污染序列放入euk.fa这个文件

首先通过bowtie2-build 构建基因组索引

bowtie2-buildeuk.fa euk

之后 使用bowtie2去除污染

bowtie2-p 40 -N 1 -x ./ euk -f <我的需要去除污染的序列> -S output.sam --un ./没有的匹配上的序列.fasta --al ./匹配上的序列.fasta

-p表示使用的线程数目

-N Setsthe number of mismatches to allowed in a seed alignment during multiseedalignment。 可以设置为 0 或 1。将此设置得更高会使对齐变慢(通常慢得多)但会增加灵敏度。 默认值:0。

跑的时候会出现如下提示(这里意味着我的序列里边没有这个污染):

参考链接:https://bowtie-bio.sourceforge.net/bowtie2/manual.shtml#output-options

Bowtie2去除污染的使用方法相关推荐

  1. 《炬丰科技-半导体工艺》化学清洗过程中重金属污染的监测方法

    书籍:<炬丰科技-半导体工艺> 文章:化学清洗过程中重金属污染的监测方法 编号:JFKJ-21-809 作者:炬丰科技 摘要 本文通过绘制少数载流子扩散长度.体中铁浓度和表面污染(表面电荷 ...

  2. 新房装修甲醛超标的危害 专业去除室内甲醛的方法

    新房装修甲醛超标的危害 专业去除室内甲醛的方法 新装修房屋是发生甲醛中毒事故的重灾区,在入住新房前必须谨慎治理好装修甲醛超标污染.空置通风可以促进房间里的甲醛排放,从而降低新房里的甲醛浓度,通风时间越 ...

  3. java 两个数组去重复数据_Java实现数组去除重复数据的方法详解

    本文实例讲述了Java实现数组去除重复数据的方法.分享给大家供大家参考,具体如下: 前一段时间被面试问到:如果一个数组中有重复元素,用什么方法可以去重?一时间会想到用一种方法,但是后来查阅资料后发现, ...

  4. win7系统去除图标箭头的方法

    应该有不少win7电脑用户有遇到桌面图标都多出来一个小箭头的情况吧,这会让我们的桌面看起来很丑,很多用户都想要删除这个小箭头却不知道怎么下手,下面分享win7系统去除图标箭头的方法. win7系统去除 ...

  5. DeOccNet:国防科大提出阵列相机去除前景遮挡成像新方法

    点击我爱计算机视觉标星,更快获取CVML新技术 阵列相机可以从不同的视角记录当前场景,并对场景的结构进行解析,因而在战场侦察.公安监视等领域具有巨大的应用潜力. 近日,国防科技大学研究团队提出阵列相机 ...

  6. 去黑边_Vegas pro 17导入的素材有黑边咋进行去除的3种方法

    在vegas中进行视频剪辑,发现导入的素材文件不能适应创建的项目窗口大小,有难看的黑边现象,特别影响美观.视频若有黑边不去除,剪辑完素材后,直接进行视频渲染导出,最终导出的视频还是会有挥之不去的难看黑 ...

  7. 自己动手去除暴风影音2012广告方法,这样才安全

    自己动手去除暴风影音2012广告方法,这样才安全.不要再用去广告补丁了,大多数去广告补丁都加了小***,破解版也如是. 1.打开暴风影音的安装文件夹,一般是X:\ Program Files \ St ...

  8. php 去除变态空格字符方法,空格trim不掉问题解决思路

    php 去除变态空格字符方法,空格trim不掉问题解决思路 参考文章: (1)php 去除变态空格字符方法,空格trim不掉问题解决思路 (2)https://www.cnblogs.com/-mrl ...

  9. vue 列表 萌层 鼠标移入移出_vue鼠标移入添加class样式,鼠标移出去除样式(active)实现方法...

    鼠标移入添加class样式 HTML HTML绑定事件,加入或者移出class为active 流量套餐 JS 这里除了active这个class需要动态添加或者减去,其他的皆是移入移出都需要的clas ...

最新文章

  1. YARN集群维护部分问题汇总
  2. idea 中使用 jetty 插件
  3. Hadoop学习系列之PageRank
  4. mysql 导入dbm文件_DBM数据导入到mysql数据库方法
  5. Docker学习总结(9)——Docker常用命令
  6. CCNA 学习笔记(七)--交换知识(VLAN VTP TRUNK)
  7. proteus数码管不亮是什么原因_数码管灌封胶有何产品性能特点?如何正确使用?...
  8. ECSHOP商品描述和文章里不加水印,只在商品图片和商品相册加水印
  9. 用Hough投票做物体检测(续)
  10. vvv在线文档导出工具_墙裂推荐 | 在线文档编辑工具
  11. SPSS 25软件安装教程
  12. python+neo4j ,Dota2自走棋,小型知识图谱构建
  13. 迅雷欲缔造互联网“视频梦工厂”
  14. DS18B20 数字温度传感器的使用和基于RT-Thread操作系统的实现
  15. [译]const T vs. T const ——Dan Saks 【翻译】
  16. GIS空间分析(五)—— 位置与几何关系分析
  17. Nginx 反向代理的知识再温习一下
  18. 13-roman-to-integer
  19. 【VOA英语学习--0207】关于AI
  20. G1.sniper Z97主板不能进入BIOS

热门文章

  1. 使用obsidian数据同步至Git,从而解决obsidian多端同步收费的问题-windows
  2. rs结果集可以传递给另一个jsp吗_从未如此高调 宝骏RS-5独家设计解读
  3. 新版FPC摄像头测评 OV7725 OV7670 OV9650 OV9655 OV5640 OV5642 OV2640 OV3640 MT9D112
  4. 房地产新项目经济测算解决方案
  5. CentOS7.0基于hadoop集群安装配置Hive
  6. Epson机械手4轴6轴示意图
  7. java在进行流操作的close()方法
  8. Dropbox VS 华为网盘
  9. 倍福---MC_ReadParameter读取参数
  10. CUIT CTF WriteUp-Code100