非模式物种10X空间转录组单细胞参考基因组构建

非模式物种10X空间转录组单细胞参考基因组构建


以下是示例代码:

#download from website: https://orchidstra2.abrc.sinica.edu.tw/orchidstra2/pagenome/padownload.php
wget -c https://orchidstra2.abrc.sinica.edu.tw/orchidstra2/pagenome/padownload/P_aphrodite_genomic_scaffold_v1.0.fa.gz
wget -c  https://orchidstra2.abrc.sinica.edu.tw/orchidstra2/pagenome/padownload/P_aphrodite_genomic_scaffold_v1.0_gene.gff.gz
#解压
gunzip *gz
#如果没有GTF文件可以转换一下
gffread  P_aphrodite_genomic_scaffold_v1.0_gene.gff  -T -o P_aphrodite_genomic_scaffold_v1.0_gene.gtf

空间转录组

#可选操作,选择蛋白编码的基因:
#spaceranger mkgtf  P_aphrodite_genomic_scaffold_v1.0_gene.gtf  P_aphrodite_genomic_scaffold_v1.0_gene.filtered.gtf  --attribute=gene_biotype:protein_coding
#10X官方参考基因组索引建立参考:https://www.10xgenomics.com/cn/support/software/space-ranger/latest/advanced/custom-references
spaceranger mkref  --nthreads 20 --genome=P_aphrodite --fasta=P_aphrodite_genomic_scaffold_v1.0.fa --genes=P_aphrodite_genomic_scaffold_v1.0_gene.gtf


单细胞:

#可选操作,选择蛋白编码的基因:
#cellranger mkgtf  P_aphrodite_genomic_scaffold_v1.0_gene.gtf  P_aphrodite_genomic_scaffold_v1.0_gene.filtered.gtf  --attribute=gene_biotype:protein_coding
#10X官方参考基因组索引建立参考:https://www.10xgenomics.com/cn/support/software/space-ranger/latest/advanced/custom-references
cellranger mkref  --nthreads 20 --genome=P_aphrodite --fasta=P_aphrodite_genomic_scaffold_v1.0.fa --genes=P_aphrodite_genomic_scaffold_v1.0_gene.gtf



  • 发表于 3天前
  • 阅读 ( 41 )
  • 分类:转录组

相关问题

0 条评论

请先 登录 后评论
omicsgene
omicsgene

生物信息

707 篇文章

作家榜 »

  1. omicsgene 707 文章
  2. 安生水 352 文章
  3. Daitoue 167 文章
  4. 生物女学霸 120 文章
  5. xun 82 文章
  6. rzx 79 文章
  7. 红橙子 78 文章
  8. CORNERSTONE 72 文章