批量下载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. 5500以太网模块pcb布局_以太网PCB布局布线
  2. 百度献礼高校开学季:AI Studio教育版上线!
  3. CyclicBarrier 和CountDownLatch使用详解
  4. 802.11ac标准简介
  5. CS231n Convolutional Neural Networks for Visual Recognition------Python Tutorial
  6. 国科大prml12-半监督学习
  7. Oracle ITL(Interested Transaction List)理解
  8. License for package Android SDK Build-Tools 28.0.3 not accepted.
  9. 野外帐篷露营避难有感
  10. Android开发之网络
  11. 数据结构经典问题——出栈顺序 转载至:canlynet微博
  12. taro 如何使用dom_taro创建的小程序项目怎么引用ant design mobile UI库
  13. 多级评论单表结构设计
  14. 用友软件计算机时间格式,如何正确设置系统日期格式?_速达软件_用友,速达,管家婆,微软,金蝶,方正,博世通,数据博士,进销存,财务软件-飞鸿软件帮助中心 -...
  15. 电源开关电源200W 12V 24V,电源架构PFC+LLC+同步整流,高效率高功率因数
  16. 20230220学习总结02
  17. 205.自定义菜单[个人公众号未认证无法使用]
  18. 单台电脑jmeter压力测试最大值
  19. connect 连接超时
  20. CentOS7-启动网卡-查看IP地址-查看端口-管理端口

热门文章

  1. 智能风扇语音控制新方案,高品质语音识别芯片NRK3303
  2. Java高并发 - JUC
  3. CAN/CANFD数据诊断记录仪解决车辆数据采集难点
  4. 中国1,4-丁烯二醇市场调研与投资预测报告(2021版)
  5. android设置水平布局,Android开发借助LinearLayout实现垂直水平居中布局
  6. 微信小程序关闭按钮样式设计
  7. 电磁式继电器使用注意事项
  8. Codeforces 156B. Suspects
  9. GA求函数极值问题MATLAB代码学习
  10. 良心建议,生活中远离这种人