首页
学习
活动
专区
工具
TVP
发布

图形化开放式生信分析系统开发

专栏成员
29
文章
31470
阅读量
22
订阅数
满分室间质评之GATK Somatic SNV+Indel+CNV+SV(下)性能优化
#此处是原先Manta分析SV的步骤一,生成runWorkflow.py,因为这一不步速度很快,所以串行执行 rm -f ${result}/${sn}/runWorkflow.py python ${tools.manta} \ --normalBam ${result}/${sn}NC_marked.bam \ --tumorBam ${result}/${sn}_marked.bam \ --referenceFasta ${refs.hum} \ --exome \ --callRegions /opt/ref/projects/Illumina_pt2.bed.zip \ --runDir ${result}/${sn} # 对bam文件碱基质量校正的第二步,Normal & Tumor并行处理 ${tools.gatk} ApplyBQSR \ --bqsr-recal-file ${result}/${sn}_recal.table \ -L ${refs.interval} \ -R ${refs.hum} \ -I ${result}/${sn}_marked.bam \ -O ${result}/${sn}_bqsr.bam & ​ ​ ${tools.gatk} ApplyBQSR \ --bqsr-recal-file ${result}/${sn}NC_recal.table \ -L ${refs.interval} \ -R ${refs.hum} \ -I ${result}/${sn}NC_marked.bam \ -O ${result}/${sn}NC_bqsr.bam & ​ #原先QC步骤,获取insert size,Normal & Tumor并行 ${tools.gatk} CollectInsertSizeMetrics \ -I ${result}/${sn}_marked.bam \ -O ${result}/${sn}_insertsize_metrics.txt \ -H ${result}/${sn}_insertsize_histogram.pdf & ​ ​ ${tools.gatk} CollectInsertSizeMetrics \ -I ${result}/${sn}NC_marked.bam \ -O ${result}/${sn}NC_insertsize_metrics.txt \ -H ${result}/${sn}NC_insertsize_histogram.pdf & ​ # 运行manta SV分析 python ${result}/${sn}/runWorkflow.py -m local -j ${envis.threads} & ​ # 运行cnvkit CNV分析 ${tools.cnvkit} batch \ ${result}/${sn}_marked.bam \ --normal ${result}/${sn}NC_marked.bam \ --method hybrid \ --targets ${refs.bed} \ --annotate /opt/ref/refFlat.txt \ --output-reference ${result}/${sn}_reference.cnn \ --output-dir ${result}/ \ --diagram \ -p 0 & ​ #samtools统计测序深度 ${tools.samtools} depth -b ${refs.bed} ${result}/${sn}_marked.bam > ${result}/${sn}_marked.depth & ${tools.samtools} depth -b ${refs.bed} ${result}/${sn}NC_marked.bam > ${result}/${sn}NC_marked.depth & #samtools统计比对信息 ${tools.samtools} flagstat --threads ${envis.threads} ${result}/${sn}_marked.bam > ${result}/$
SliverWorkspace
2020-08-04
1.9K0
图形化开放式生信分析系统开发 - 7 分析报告的模板定制与自动生成
在前文中描述了生信分析流程的整个设计、图形化、服务器端运行、批量运行与过程控制(手动停止,恢复运行、覆盖运行),本文尝试解决分析报告的模板定制与自动生成。
SliverWorkspace
2020-01-17
7290
图形化开放式生信分析系统开发 - 5 生信分析流程服务器端运行
在上文图形化开放式生信分析系统开发 - 4 生信分析流程的图形化设计 讨论了生信分析pipeline的图形化,如何用图形的方式显示生信pipeline,但是pipeline脚本按照变量的形式保存之后,如何运行,在什么环境下运行?是本文要解决的问题。
SliverWorkspace
2020-01-17
7940
图形化开放式生信分析系统开发 - 4 生信分析流程图形化
在上文图形化开放式生信分析系统开发 - 3 生信分析流程的进化 讨论了生信分析pipeline的进化,从手动到自动,但仍然停留在终端命令行阶段,为了让更多非生信专业的人能够使用,就要想办法实现生信分析pipeline的图形化:提供能够快速上手的UI,简单点击鼠标就能够运行的图形化系统。
SliverWorkspace
2020-01-17
8670
没有更多了
社区活动
【纪录片】中国数据库前世今生
穿越半个世纪,探寻中国数据库50年的发展历程
Python精品学习库
代码在线跑,知识轻松学
博客搬家 | 分享价值百万资源包
自行/邀约他人一键搬运博客,速成社区影响力并领取好礼
技术创作特训营·精选知识专栏
往期视频·千货材料·成员作品 最新动态
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档