以下是示例代码:
#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
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!