首页 > 其他分享 >00 去宿主DNA--bowtie2

00 去宿主DNA--bowtie2

时间:2023-07-04 21:01:21浏览次数:33  
标签:Acc 00 DNA -- SRR 索引 wz 10t bowtie2

1、构建索引

bowtie2-build maizev4.fa maizev4

使用 Bowtie 2 的 bowtie2-build 命令来构建一个maizev4基因组的索引。bowtie2-build 是 Bowtie 2 中用于创建索引的工具,而 Bowtie 2 是一个用于短序列比对的快速、内存高效的工具。下面是这条命令的分析:

  • bowtie2-build: Bowtie 2 中用于创建索引的命令。
  • maizev4.fa: 输入的 FASTA 格式的基因组文件。该文件包含了maizev4基因组的所有序列信息。
  • maizev4: 为创建的索引文件指定的基本名。构建好的索引文件将以此为前缀。

执行完这条命令后,将会生成一系列以 "maizev4" 为前缀的 Bowtie 2 索引文件。这些索引文件可用于后续的序列比对任务。

2、DNA序列比对和分析

bowtie2 -x /10t/wz/SRR_Acc/index/index -1 /10t/wz/SRR_Acc/PRJNA/SRR8185562_1.fastq -2 /10t/wz/SRR_Acc/PRJNA/SRR8185562_2.fastq --very-sensitive-local -p 35 -S /10t/wz/SRR_Acc/bowtie2/SRR8185562.sam --un-conc /10t/wz/SRR_Acc/bowtie2/SRR8185562_unmap 帮我分析一下这条语句
  • bowtie2: 这是Bowtie2程序的名称,用于DNA序列比对。
  • -x /10t/wz/SRR_Acc/index/index-x选项后面是索引文件的路径,即待比对的参考基因组的索引。
  • -1 /10t/wz/SRR_Acc/PRJNA/SRR8185562_1.fastq-1选项后面是第一条PE(paired-end)测序数据的路径。
  • -2 /10t/wz/SRR_Acc/PRJNA/SRR8185562_2.fastq-2选项后面是第二条PE测序数据的路径。
  • --very-sensitive-local: 这是Bowtie2的参数选项,使用非常敏感的本地模式进行比对。
  • -p 35-p选项后面是并行线程数,此处设置为35,表示使用35个线程进行比对。
  • -S /10t/wz/SRR_Acc/bowtie2/SRR8185562.sam-S选项后面是输出比对结果的SAM格式文件的路径。
  • --un-conc /10t/wz/SRR_Acc/bowtie2/SRR8185562_unmap--un-conc选项后面是未比对上的read对的输出文件路径,这里是将未比对上的paired reads输出到指定文件中。

通过执行这条命令,Bowtie2将使用参考基因组的索引进行敏感本地比对,并输出比对结果到SAM文件中。同时,未比对上的paired reads会被输出到指定的文件中。此外,该命令还使用35个线程来加速比对过程。

 

   

标签:Acc,00,DNA,--,SRR,索引,wz,10t,bowtie2
From: https://www.cnblogs.com/wzbzk/p/17526982.html

相关文章

  • 【docker】docker的环境变量怎么写
    背景docker启动的时候,有时候需要修改一些参数,就是那些--envkey=vvvvvvvv这种的,那这个key是怎样定的呢?比如现在有一个eureka,在springboot项目里面的bootstrap.yml文件有这么一段eureka:instance:leaseRenewalIntervalInSeconds:10##心跳检测时长pr......
  • 一文读懂:完整的支付系统整体架构!
    支付产品模块是按照支付场景来为业务方提供支付服务。这个模块一般位于支付网关之后,支付渠道之前。它根据支付能力将不同的支付渠道封装成统一的接口,通过支付网关来对外提供服务。所以,从微服务的角度,支付产品本身也是一个代理模式的微服务,它透过支付网关响应业务方请求,进行一......
  • 最全Java锁详解:独享锁/共享锁+公平锁/非公平锁+乐观锁/悲观锁
    最全Java锁详解:独享锁/共享锁+公平锁/非公平锁+乐观锁/悲观锁在Java并发场景中,会涉及到各种各样的锁如公平锁,乐观锁,悲观锁等等,这篇文章介绍各种锁的分类:公平锁/非公平锁可重入锁独享锁/共享锁乐观锁/悲观锁分段锁自旋锁 最全Java锁详解:独享锁/共享锁+公平锁/非......
  • Python基础37 基于tcp、udp套字编程、粘包现象、struct模块
    基于tcp协议的套接字编程(sochet编程)什么是socket?通常翻译为套接字,socket是在应用层和传输层之间的一个抽象层,它把tcp/ip层复杂的操作抽象为几个简单的接口供应用层调用已实现进程在网络中。套接字分类:AF_UNIX:用在局域网中AF_INET:用在互联网中客户端和服务端启动顺序......
  • JAVA生产环境性能监控与调优看这一篇就够了
     JVM的内存结构 JVM的参数类型标准参数(各版本中保持稳定)-help-server-client-version-showversion-cp-classpathX参数(非标准化参数)-Xint:解释执行-Xcomp:第一次使用就编译成本地代码-Xmixed:混合模式,JVM自己决定是否编译成本地代码示例:java-version(默认是......
  • TryHackMe | KenobiWriteup
    TryHackMe|KenobiWriteupTask1DeploythevulnerablemachineScanthemachinewithnmap,howmanyportsareopen?运行下nmap一眼就是7个Task2EnumeratingSambaforsharesSamba是适用于和Unix的标准Windows互操作性程序套件。它允许最终用户访问和使......
  • 钡铼技术多功能RTU S475多功能RTU改变养殖行业现殖效率
    在养殖行业中,对环境参数的精确监测与控制至关重要。然而,传统的监测方法往往存在诸多痛点,如数据采集不准确、传输速度慢、可视化效果差等。为了解决这些问题,钡铼技术公司推出了其旗舰产品——S475多功能RTU,该产品在养殖行业监测中展现出了显著的优势。钡铼S475多功能RTU是一款......
  • 键盘行
    给你一个字符串数组words,只返回可以使用在美式键盘同一行的字母打印出来的单词。键盘如下图所示。美式键盘中:第一行由字符"qwertyuiop"组成。第二行由字符"asdfghjkl"组成。第三行由字符"zxcvbnm"组成。示例1:输入:words=["Hello","Alaska","Dad","Peace"]输......
  • 微软确认了 Windows 11 存在 CPU 占用异常的问题
    导读微软近日确认了Windows11系统中所存在的一个错误,该问题影响了运行Windows1121H2和22H2版本的设备,换句话说就是所有发布的Windows11操作系统。该问题是在用户安装了KB5026372(2023年5月补丁)后产生的,该问题破坏了文件资源管理器的一个特定部分,导致CPU......
  • 数据库问题之“字符编码问题 Cause: java.sql.SQLException: Incorrect string value:
     1)表1和表2的产品名称[数据库字段]字符编译方式不一致①问题 org.springframework.jdbc.UncategorizedSQLException:Errorupdatingdatabase.Cause:java.sql.SQLException:Incorrectstringvalue:'\xF0\x9F\x8E\x81\xE7\x88...'forcolumn'product_name'atr......