从NCBI下载牛基因组,下载后gunzip解压

wget ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/002/263/795/GCF_002263795.1_ARS-UCD1.2/GCF_002263795.1_ARS-UCD1.2_genomic.fna.gz

安装bowtie2

sudo apt-get install bowtie2

构建索引数据库

bowtie2-build GCF_002263795.1_ARS-UCD1.2_genomic.fna cattlegenome #cattlegenome是输出,fna文件是fa格式的输入

比对。-p线程数,-x参考基因组的基本名(除了后缀名),-1和-2配对双端序列,-S输出sam文件,--un-conc 输出未注释成功的结果,双端结果会以.1和.2区分。--un-conc-gz输出gz压缩序列

bowtie2 -p 15 -x /mnt/10t/database/cattle/cattlegenome -1 1.clip.1.fq.gz -2 1.clip.2.fq.gz -S sample.sam --un-conc 1.fq 

新建sh文件夹和dehost文件夹(略)

生成所有需要的文件


for i in 123 98 126 4 128 32 134 8 130 29 131 42 119 19 139 24 125 37 144 14 117 5 142 13 127 22 129 48 138 47 140 20 132 18 149 73 120 105 148 44 122 30 124 26 143 34 133 12 145 21 152 10 118 102 147 80 121 1 151 35 135 38 136 2 150 17 146 55 141 11 137 97; do echo "bowtie2 -p 15 -x /mnt/10t/database/cattle/cattlegenome -1 /mnt/10t/mzy/rawdata/mazhiyuan/$i.clip.1.fq.gz -2 /mnt/10t/mzy/rawdata/mazhiyuan/$i.clip.2.fq.gz --un-conc-gz /mnt/10t/mzy/rawdata/mazhiyuan/dehost/$i.fq.gz" >/mnt/10t/mzy/rawdata/mazhiyuan/sh/$i.bowtie.sh; done

所有命令一起运行

for i in *.sh; do nohup bash $i & done

利用bowtie2去除宿主基因相关推荐

  1. 利用多组学整合鉴定人类疾病共享的和疾病特异性的宿主基因-微生物组关联

    利用多组学整合鉴定人类疾病共享的和疾病特异性的宿主基因-微生物组关联 2022-05-23 本文系谷歌翻译 原文:Identification of shared and disease-specif ...

  2. Cell Reports:去除宿主和胞外DNA以提高微生物基因组得率(痰液样本)

                    简介                  标题:Cell Reports:去除宿主和胞外DNA以提高微生物基因组得率(痰液样本) 杂志:Cell Reports 影响因子 ...

  3. 广东生态所孙蔚旻团队EST发表利用稳定同位素示踪-宏基因组分箱联用技术揭示砷污染土壤中的厌氧砷氧化微生物及其代谢途径...

    广东省生态环境技术研究所孙蔚旻团队ES&T发表:利用稳定同位素示踪-宏基因组分箱联用技术揭示砷污染土壤中的厌氧砷氧化微生物及其代谢途径 第一作者:张苗苗 通讯作者:孙蔚旻 通讯单位:广东省生态 ...

  4. 无法在只读编辑器中编辑_Mol Plant中国农科院作物科学研究所夏兰琴课题组成功利用 优化的引导基因编辑器在水稻中实现高效精准基因编辑...

    近日,中国农业科学院作物科学研究所夏兰琴实验室利用水稻密码子优化的nCas9 (H840A)蛋白和逆转录酶突变体(M-MLV-RT)的融合基因,开发了一种高效的植物引导编辑(Prime Editing ...

  5. 去掉java文件中的注释_利用JavaParser去除java文件中的注释

    利用JavaParser去除java文件中的注释 个人博客:记录一下在项目实施过程中的一些点 情景回顾 之前项目有个需求,就是去掉.java文件中的所有注释,常用的方法是用正则匹配.然而在网络上查找到 ...

  6. 计算机描绘的基因结构图,利用IBS软件画基因mRNA的结构图

    原标题:利用IBS软件画基因mRNA的结构图 IBS(Illustrator for Biological Sequences)是一个简单又强大的绘图软件,它可以几分钟内完成图形的绘制.同时,还拥有多 ...

  7. html特殊字符p如何屏蔽 asp,利用asp去除html标记的四种方法

    利用asp去除html标记的四种方法 发布日期:2015-10-8 16:10:45 利用asp去除html标记的四种方法 方法一 : 使用"" 如果您想要知道如何从文本中删除ht ...

  8. 如何利用ps去除图片水印和绘制图案

    欢迎点击「算法与编程之美」↑关注我们! 本文首发于微信公众号:"算法与编程之美",欢迎关注,及时了解更多此系列文章. 1 问题描述 在where2go团队中,18级的任务主要是画原 ...

  9. es like and or_广东生态所孙蔚旻团队ESamp;T发表利用稳定同位素示踪宏基因组分箱联用技术揭示砷污染土壤中的厌氧砷氧化微生物及其代谢途径...

    广东省生态环境技术研究所孙蔚旻团队ES&T发表:利用稳定同位素示踪-宏基因组分箱联用技术揭示砷污染土壤中的厌氧砷氧化微生物及其代谢途径 第一作者:张苗苗 通讯作者:孙蔚旻 通讯单位:广东省生态 ...

最新文章

  1. html 物理引擎,在物理引擎中画圆弧
  2. 重置MYSQL的root 密码
  3. android内存优化方法,Android开发内存优化注意事项和方法
  4. 当前深度神经网络模型压缩和加速都有哪些方法?
  5. Maven-学习笔记02【基础-Maven的安装和仓库种类】
  6. Android学习笔记:对Android应用进行单元测试
  7. WebMagic学习总结
  8. 诺基亚困境,力掐苹果?
  9. 【python】一道LeetCode搞懂递归算法!#131分割回文串 #以及刷LeetCode的一点点小心得 [数据结构与算法基础]
  10. 某公司电子商务网站策划方案
  11. Word中使用表格排版公式时,表格内序号纵向居中的问题。
  12. 主板24pin接口详图_老电源也兼容 简单DIY 24Pin转接线
  13. 中望cad自定义快捷键命令_1分钟成为CAD设计高手:中望CAD命令快捷键设置详解-快捷键设置...
  14. mozhe靶场——网络安全——投票系统程序设计缺陷分析
  15. json for modern c++,这是我用过的最好用的json库了
  16. css表格随内容变化,css 使表格随着内容自动适应宽度
  17. 计算机初中几年级考,2021年小升初考试时间(2021小升初考试是几月几号 )
  18. @EventListener注解使用及源码解析
  19. 软件测是中W模型和V模型 H模型X模型的优缺点
  20. 此错误(HTTP 500 内部服务器错误)意味着您正在访问的网站出现了服务器问题,此问题阻止了该网页的显示

热门文章

  1. Android系统为何不流畅?QNX道破玄机
  2. c语言实现两个大数相乘
  3. php 中使用cURL发送get/post请求,上传图片,批处理
  4. 动环监控系统生产厂家,动环监控有哪些厂家
  5. 浅谈bluebird包的使用理解
  6. 决策树中结点的特征选择方法
  7. 充电该知道的部分知识
  8. 跳出forEach循环
  9. python255多少钱_CS255留学生作业代做、代写Python语言作业、Timetabling作业代做、Python编程设计作业调试...
  10. 2023年湖北工程师职称证书如何查询真伪?启程别告诉你