SAW 空间转录组参考基因组制作

SAW 空间转录组参考基因组制作


ensembl下载基因组序列:

wget -c https://ftp.ensembl.org/pub/release-114/fasta/mus_musculus/dna/Mus_musculus.GRCm39.dna.toplevel.fa.gz
wget -c https://ftp.ensembl.org/pub/release-114/gtf/mus_musculus/Mus_musculus.GRCm39.114.chr.gtf.gz
gunzip Mus_musculus.GRCm39.114.chr.gtf.gz
#保留染色体部分序列
seqkit grep -r -p "^[0-9]+" -p "X" -p "Y" -p "MT" Mus_musculus.GRCm39.dna.toplevel.fa.gz >Mus_musculus.GRCm39.dna.toplevel.fa

saw makeRef \
    --mode=STAR \
    --fasta=Mus_musculus.GRCm39.dna.toplevel.fa \
    --gtf=Mus_musculus.GRCm39.114.chr.gtf \
    --genome=MM_GRCm39
  • 发表于 2天前
  • 阅读 ( 22 )
  • 分类:转录组

0 条评论

请先 登录 后评论
omicsgene
omicsgene

生物信息

728 篇文章

作家榜 »

  1. omicsgene 728 文章
  2. 安生水 360 文章
  3. Daitoue 167 文章
  4. 生物女学霸 120 文章
  5. xun 89 文章
  6. rzx 85 文章
  7. 红橙子 81 文章
  8. CORNERSTONE 72 文章