批量下载SRA文件:

提交到服务器版本

vim for.shcript

#BSUB -n 16

#BSUB -q priority

#BSUB -o for.log

for i in `cat sra_acc.txt`;do

x=$(echo $i | cut -b1-6)

y=$(echo $i | cut -b1-3)

ascp -T -i ~/.aspera/connect/etc/asperaweb_id_dsa.openssh -k 1 -T

-l200m anonftp@ftp-private.ncbi.nlm.nih.gov:/sra/sra-instant/reads/ByRun/sra/$y/$x/$i/$i.sra

.

done

运行命令:bsub < for.script

Shell版本:

vim for.sh.........

for i in `cat sra_acc.txt`;do

x=$(echo $i | cut -b1-6)

y=$(echo $i | cut -b1-3)

ascp -T -i ~/.aspera/connect/etc/asperaweb_id_dsa.openssh -k 1 -T

-l200m anonftp@ftp-private.ncbi.nlm.nih.gov:/sra/sra-instant/reads/ByRun/sra/$y/$x/$i/$i.sra

.

done

运行命令:bash for.sh

其中sra_acc.txt为SRA数据的列表文件。

批量质控脚本:

提交到服务器上的脚本

vim fastq-dump.script

#BSUB -n 16

#BSUB -o fastq-dump.log

for i in *sra

do

echo $i

/gss1/biosoft/sratoolkit.2.8.2-1-centos_linux64/bin/fastq-dump

--split-3 --gzip $i

done

~

运行命令:bsub

shell版本:

vim fastq-dump.sh

for i in *sra

do

echo $i

/gss1/biosoft/sratoolkit.2.8.2-1-centos_linux64/bin/fastq-dump

--split-3 --gzip $i

done

运行命令:bash fastq-dump.sh

批量转比对脚本:

#BSUB -q priority

#BSUB -o hisat.log

for i in *.fastq.gz;do

names=${i%_*}

UB -n 16

#BSUB -q priority

##BSUB

如何利用linux分析转录组数据库,转录组分析(8)----批量处理脚本了解一下相关推荐

  1. linux 进入redis 数据库,Linux下Redis数据库的安装方法与自动启动脚本分享

    安装Redis (1) 下载Redis wget http://redis.googlecode.com/files/redis-2.2.11.tar.gz tar xzvf redis-2.2.11 ...

  2. 阿里云分析型数据库使用教程

    产品概述 随着企业IT和互联网系统的发展,越来越多的数据被产生了.数据的量的积累带来了质的飞跃,使得数据系统从业务系统的一部分演变得愈发独立,通过对数据的分析和挖掘产生自己独特的价值. 在业务系统中, ...

  3. (52)性能分析 ---CPU的性能分析

    性能分析 -CPU的性能分析 服务器中一个硬件 什么是服务器? 由 计算机硬件 + 操作系统 + 应用服务构成,提供能力输出. 服务器的性能分析,需求分析哪些? 硬件性能分析 操作系统性能分析 项目性 ...

  4. 转录组分析_高级转录组分析和R数据可视化

    封面来源:https://www.zhihu.com/question/304747766 常规转录组是我们最常接触到的一种高通量测序数据类型,其实验方法成熟,花费较低,是大部分CNS必备的技术,以后 ...

  5. 热点综述 | 单细胞+空间转录组的整合分析方法总结

    目前scRNA-seq将每个转录物与单个细胞相关联,但关于这些转录物在组织中的位置信息丢失了:相反的,空间转录组学技术知道转录物的位置,却不知道是哪个细胞产生了转录物.因此,scRNA-seq与空间转 ...

  6. 单细胞测序分析及单细胞转录组

    单细胞测序技术的应用与数据分析.单细胞转录组为主题,精心设计了具有前沿性.实用性和针对性强的理论课程和上机课程.培训邀请的主讲人均是有理论和实际研究经验的人员.学员通过与专家直接交流,能够分享到这些顶 ...

  7. 单细胞转录组数据整合分析专题研讨会(2019.11)

    2019年10月9日,单细胞转录组再等Nature.题为Decoding human fetal liver haematopoiesis的研究,对受孕后4周至17周的人胚胎肝脏.卵黄囊.肾脏和皮肤组 ...

  8. SCS【7】单细胞转录组之轨迹分析 (Monocle 3) 聚类、分类和计数细胞

    点击关注,桓峰基因 桓峰基因公众号推出单细胞系列教程,有需要生信分析的老师可以联系我们!首选看下转录分析教程整理如下: Topic 6. 克隆进化之 Canopy Topic 7. 克隆进化之 Car ...

  9. oracle服务器的操作系统,Oracle Linux 操作系统及数据库的时区机制分析

    Oracle Linux 操作系统及数据库的时区机制分析 1. /etc/localtime 这个文件记录的是系统的时区,缺省的数据库由此获得时区信息 这个文件是二进制文件,修改该文件的方法是拷贝/u ...

最新文章

  1. React Native常用组件之ListView
  2. 985硕博导师:如果你没空看书,我强烈建议你做这件事
  3. mac git命令按tab键自动补全
  4. 拓扑排序--关键路径
  5. spring23:Aspectj实现异常通知@AfterThrowing
  6. HiveSQL运行优化参数配置
  7. 使用jieba 提取 关键词
  8. C#-常用对象-思维导图
  9. 《WinForm开发系列之控件篇》Item3 BindingSource (暂无)
  10. GUI学习之二十一——QSlider、QScroll、QDial学习总结
  11. npm install -g @vue/cli时 -4048 npm ERR! Error: EPERM: operation not permitted, lstat报错的几种解决方案
  12. 小程序源码:全新实用工具证件照制作微信小程序源码下载支持多种证件生成与制作
  13. 阿里云域名解析详细步骤讲解
  14. 神经元如何将视觉世界映射到人脑?
  15. vs2017 中项目的publish,即“发布”到底是什么?
  16. HEU The message
  17. python 异步io 写excel_python异步IO编程(二)
  18. Android原生系统开发如何优雅的提供系统级的API供第三方程序调用?
  19. sprintf 函数的详解
  20. js将数字转换成万、亿、万亿

热门文章

  1. pandas excelwriter遇到的坑
  2. jsp页面计算梯形面积
  3. 在web页面上快速生成二维码的三种实用方法
  4. ISP算法学习之LSC(镜头阴影校正)
  5. 好听的计算机用户名称,好听的电脑英文用户名
  6. iOS 马甲包过包注意事项
  7. python 画矩形树图
  8. MyJupyter,一款支持Python和Java的可移动Jupyter软件包
  9. 【定量分析、量化金融与统计学】纵向数据分析(1)——截面数据、时间序列数据和纵向数据
  10. 呼叫中心的语音导航功能