Inkit wiki

用户工具

站点工具

您在这里: Wiki » 统计 » Bioinformatic » majiq

侧边栏

导航

常用社区



统计

Notebook 环境搭建

手册

统计学问题
Bioinformatic
RNA-Seq
SPSS
统计词典
写作相关
Stock


Ports
navbar

stat:bioinform:majiq

这是本文档旧的修订版！

/mnt/smb/notebooktian/RNASeq/260429AGING

### References * MAJIQ [RNA splicing analysis using heterogeneous and large RNA-seq datasets - nature communications](https://www.nature.com/articles/s41467-023-36585-y) * MAJIQ → cryptic 3'splice site

* A3SS 图B， 有拖尾感的sashimi [Cancer-Associated SF3B1 Hotspot Mutations Induce Cryptic 3′ Splice Site Selection through Use of a Different Branch Point](https://www.cell.com/cell-reports/fulltext/S2211-1247(15)01078-5)

### MAJIQ

https://bitbucket.org/biociphers/majiq_academic/src/main/

https://biociphers.bitbucket.io/majiq-docs/

alias majiq='docker run –rm -e MAJIQ_LICENSE_FILE=/data/majiq_license_academic_official.lic -v “$(pwd)/datamj:/data” -w /data majiq_image majiq' alias voila='docker run –rm -e MAJIQ_LICENSE_FILE=/data/majiq_license_academic_official.lic -v “$(pwd)/datamj:/data” -p 5000:5000 -w /data majiq_image voila'

alias majiq3='docker run –rm -e MAJIQ_LICENSE_FILE=/data/majiq_license_academic_official.lic -v “$(pwd)/datamj:/data” -w /data majiq_image majiq-v3'

* 放入 ~/majiq_license_academic_official.lic ``` gAAAAABlW6-jYooOgeYvzqoEZ9Z6jMLcFAfydEzRxJWwJ7bl9z71xmQCS343USkuLpu88lpTwTOGiGqpJ0uMZV36nto1SCRmdP_LWSugZyS7-E1wDVOzyV8kQcUHzNtPQVMVaLUQXmpYGfXGHEVo-XUPVG3v3waurw== ```

```bash apt install build-essential apt install libhts-dev

```

majiq build 这是最耗时的一步。它会扫描你的 BAM 文件并识别所有的剪接事件（LSVs）。 majiq build /data/mm39_115.gff3 -c ./majiq.config -o ./build_output

-v “$(pwd):/data”，所以你在外面（宿主机）看到的文件如果是 ~/project/data/sample1.bam，在配置文件里应该直接写 sample1（如果你就在 data 目录下运行的话），或者相对于 /data 的路径。

### Steps v3 ```bash alias voila='docker run –rm -e MAJIQ_LICENSE_FILE=/data/majiq_license_academic_official.lic -v “$(pwd)/datamj:/data” -p 5000:5000 -w /data majiq_image voila' alias majiq3='docker run –rm -e MAJIQ_LICENSE_FILE=/data/majiq_license_academic_official.lic -v “$(pwd)/datamj:/data” -w /data majiq_image majiq-v3'

PATH_GFF3=“/data/mm39_115.gff3” PATH_BAMS=“/data/bams” PATH_ANNO_SG=“/data/mj/anno_sg.zarr” PATH_SJ=“/data/mj/sjs” PATH_SG_ZAAR=“/data/mj/sg.zarr” PATH_PSICOV=“/data/mj/psi” PATH_SGC=“/data/mj/sgc” PATH_VOILA=“/data/mj/voila”

majiq3 gff3 $PATH_GFF3 $PATH_ANNO_SG majiq3 sj $PATH_BAMS/OLDSED_1.bam $PATH_ANNO_SG $PATH_SJ/OLDSED_1.sj majiq3 sj $PATH_BAMS/OLDSED_2.bam $PATH_ANNO_SG $PATH_SJ/OLDSED_2.sj majiq3 sj $PATH_BAMS/OLDSED_3.bam $PATH_ANNO_SG $PATH_SJ/OLDSED_3.sj majiq3 sj $PATH_BAMS/YNGSED_1.bam $PATH_ANNO_SG $PATH_SJ/YNGSED_1.sj majiq3 sj $PATH_BAMS/YNGSED_2.bam $PATH_ANNO_SG $PATH_SJ/YNGSED_2.sj majiq3 sj $PATH_BAMS/YNGSED_3.bam $PATH_ANNO_SG $PATH_SJ/YNGSED_3.sj

majiq3 build $PATH_ANNO_SG $PATH_SG_ZAAR –sjs $PATH_SJ/OLDSED_1.sj $PATH_SJ/OLDSED_2.sj $PATH_SJ/OLDSED_3.sj –sjs $PATH_SJ/YNGSED_1.sj $PATH_SJ/YNGSED_2.sj $PATH_SJ/YNGSED_3.sj

majiq3 psi-coverage $PATH_SG_ZAAR $PATH_PSICOV/OLDSED.psicov $PATH_SJ/OLDSED_1.sj $PATH_SJ/OLDSED_2.sj $PATH_SJ/OLDSED_3.sj majiq3 psi-coverage $PATH_SG_ZAAR $PATH_PSICOV/YNGSED.psicov $PATH_SJ/YNGSED_1.sj $PATH_SJ/YNGSED_2.sj $PATH_SJ/YNGSED_3.sj

majiq3 sg-coverage $PATH_SG_ZAAR $PATH_SGC/OLDSED.sgc $PATH_SJ/OLDSED_1.sj $PATH_SJ/OLDSED_2.sj $PATH_SJ/OLDSED_3.sj majiq3 sg-coverage $PATH_SG_ZAAR $PATH_SGC/YNGSED.sgc $PATH_SJ/YNGSED_1.sj $PATH_SJ/YNGSED_2.sj $PATH_SJ/YNGSED_3.sj

majiq3 deltapsi –splicegraph $PATH_SG_ZAAR –output-voila $PATH_PSICOV/OLDSED_v_YNGSED.psicov –output-tsv $PATH_PSICOV/OLDSED_v_YNGSED.tsv -psi1 $PATH_PSICOV/OLDSED.psicov -psi2 $PATH_PSICOV/YNGSED.psicov

majiq3 quantify $PATH_PSICOV/OLDSED.psicov –min-experiments 0.01 –splicegraph $PATH_SG_ZAAR –output-tsv $PATH_PSICOV/OLDSED.tsv –overwrite majiq3 quantify $PATH_PSICOV/YNGSED.psicov –min-experiments 0.01 –splicegraph $PATH_SG_ZAAR –output-tsv $PATH_PSICOV/YNGSED.tsv –overwrite

majiq3 heterogen –stats infoscore mannwhitneyu ttest tnom –splicegraph $PATH_SG_ZAAR –output-voila $PATH_VOILA/OLDSED_v_YNGSED.voila –output-tsv $PATH_PSICOV/OLDSED_v_YNGSED.tsv -psi1 $PATH_PSICOV/OLDSED.psicov -psi2 $PATH_PSICOV/YNGSED.psicov

voila view –host 0.0.0.0 –port 5000 –enable-type-indexing –psicov-grouping-file /data/majiq.tsv $PATH_SG_ZAAR $PATH_VOILA/OLDSED_v_YNGSED.voila $PATH_SGC/

```

### Steps majiq build /data/mm39_115.gff3 -c /data/majiq.config -o /data/build_output majiq psi-coverage \ 

  /data/build_output/splicegraph.zarr \
  /data/psi_coverage_output \
  /data/build_output/YNGSED_1.sj /data/build_output/YNGSED_2.sj /data/build_output/YNGSED_3.sj \
  /data/build_output/OLDSED_1.sj /data/build_output/OLDSED_2.sj /data/build_output/OLDSED_3.sj \
  --nthreads 8

majiq deltapsi \

  1. grp1 /data/psi_coverage_output \
  2. grp2 /data/psi_coverage_output \
  3. -select-grp1-prefixes YNGSED_1 YNGSED_2 YNGSED_3 \
  4. -select-grp2-prefixes OLDSED_1 OLDSED_2 OLDSED_3 \
  5. -names YNGSED OLDSED \
  6. -splicegraph /data/build_output/splicegraph.zarr \
  7. -output-voila /data/deltapsi_output/YNGSED_vs_OLDSED.voila \
  8. -output-tsv /data/deltapsi_output/YNGSED_vs_OLDSED.tsv \
  9. -nthreads 8

majiq sg-coverage /data/build_output/splicegraph.zarr /sgc/YNGSED.sgc /data/build_output/YNGSED_1.sj /data/build_output/YNGSED_2.sj /data/build_output/YNGSED_3.sj

voila view /data/build_output/ –host 0.0.0.0 –port 5000 –enable-type-indexing –psicov-grouping-file /data/majiq.tsv

stat/bioinform/majiq.1777556392.txt.gz · 最后更改: 由 inkit

页面工具

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki