jbrowse基因组浏览器部署

一、jbrowse部署

1.1 jbrowse部署

复制代码
sudo mkdir /var/www/jbrowse;
sudo chown `whoami` /var/www/jbrowse;
# cd into it
cd /var/www/jbrowse;

1.2 nginx配置

复制代码
  server {
        listen       80 default_server;
        listen       [::]:80 default_server;
        server_name  _;
        #root         /var/www/html/jbrowse;

        # Load configuration files for the default server block.
        include /etc/nginx/default.d/*.conf;

        location / {
        root  /var/www/html/jbrowse;
        autoindex on;
        index  index.html index.htm  index.jsp;
        }

        error_page 404 /404.html;
            location = /40x.html {
        }

fetch a JBrowse release zip file

复制代码
curl -O http://jbrowse.org/releases/JBrowse-1.8.0.zip

# unzip it and cd into it
unzip JBrowse-1.8.0.zip
cd JBrowse-1.8.0

二 数据导入

1、 ca文件处理

samtools工具安装

复制代码
 wget https://github.com/samtools/samtools/releases/download/1.9/samtools-1.9.tar.bz2
  tar -jxvf samtools-1.9.tar.bz2 
  cd samtools-1.9/
  ./configure --prefix=/usr/local/samtools #编译
  make 
  make install

samtools faidx zdd0006.genome.v1.fa #生成ca文件的索引
执行命令后会自动生成文件zdd0006.genome.v1.fa.fai

2、vcf格式文件处理

确保你的 VCF 文件是经过 BGZF 压缩的(以 .vcf.gz 结尾),并且已经创建了相应的索引文件(.tbi 文件)。如果没有,你可以使用 bgzip 和 tabix 工具来进行压缩和索引。

tabix命令安装

复制代码
git clone https://github.com/samtools/htslib.git
 git submodule update --init --recursive

然后进入项目make 和make install

数据处理

复制代码
 bgzip -c 680_sample_filter_ID.vcf > 680_sample_filter_ID.vcf.gz
 tabix -p vcf 680_sample_filter_ID.vcf.gz

track.conf文件配置

复制代码
[root@VM-4-17-opencloudos data]# cat tracks.conf 
[GENERAL]
refSeqs=Gmax_275_v2.0.fa.fai
[tracks.Gmax_275_v2.0]
urlTemplate=Gmax_275_v2.0.fa
storeClass=JBrowse/Store/SeqFeature/IndexedFasta
type=Sequence

refSeqs=SN14.genome.fa.fai
[tracks.SN14_genome]
urlTemplate=SN14.genome.fa
storeClass=JBrowse/Store/SeqFeature/IndexedFasta
type=Sequence

refSeqs=zdd0006.genome.v1.fa.fai
[tracks.zdd0006_genome]
urlTemplate=zdd0006.genome.v1.fa
storeClass=JBrowse/Store/SeqFeature/IndexedFasta
type=Sequence


[tracks.680_sample_filter]
# settings for what data is shown in the track
storeClass     = JBrowse/Store/SeqFeature/VCFTabix
urlTemplate    = 680_sample_filter_ID.vcf.gz

 
# settings for how the track looks
category = VCF
type = JBrowse/View/Track/CanvasVariants
key  = 680_sample
#GFF格式
[tracks.Gmax_275_Wm82]
urlTemplate = Gmax_275_Wm82.a2.v1.gene.gff3
storeClass = JBrowse/Store/SeqFeature/GFF3
type = CanvasFeatures
key = Gmax_275_Wm82
相关推荐
pftn1 小时前
从0到1搭建具身智能数据飞轮:Ray + Kafka + LanceDB 全链路实战
大数据
Web3_Daisy2 小时前
什么是 Robinhood Chain?
大数据·人工智能·web3·区块链
无忧智库3 小时前
智慧银行反欺诈大数据管控平台建设方案:从规则拦截到图谱智能,重构银行实时风控中枢(PPT)
大数据
Elasticsearch3 小时前
Elastic Observability 中的 Prometheus 指标:你的 PromQL 可以保持不变地运行
elasticsearch
Elasticsearch3 小时前
为什么 Elasticsearch 平台是你的 AI 技术栈中缺失的一块拼图
elasticsearch
湘美书院--湘美谈教育4 小时前
湘美谈教育湘美书院大湘西文学系列:AI时代的武侠小说怎么写
大数据·人工智能·深度学习·机器学习·生活
VortMall4 小时前
『平台去经营化』平台治理能力全新重构|VortMall微服务商城系统v1.3.10
java·大数据·微服务·商城系统·开源商城·vortmall·去经营化
workflower4 小时前
模型中心主义路径
大数据·人工智能·深度学习·机器学习·自动化·制造
Elasticsearch4 小时前
数据平台赌注:为什么金融 AI 计划停滞,以及成功者如何扩展规模
elasticsearch
Urbano4 小时前
卫衣全工序自动化智造科普:替代工位、设备选型与产能升级方案
大数据·人工智能·自动化