比对软件之STAR的使用方法
1 | source activate /home/lizhixin/softwares/anaconda3/envs/splicing |
建索引
1 2 | mkdir GFP_index STAR --runThreadN 6 --runMode genomeGenerate --genomeDir GFP_index --genomeFastaFiles GFP.fasta |
普通比对
1 2 3 4 5 6 7 8 | STAR --genomeDir /home/lizhixin/project/scRNA-seq/GFP_YFP_tracing/reference/all .GFPs.ref \ --runThreadN 10 \ --readFilesIn 7Ala-D60-BO-2-1_S25_L003_R1_001.fastq.gz,7Ala-D60-BO-2-1_S25_L003_R2_001.fastq.gz UE-D60-BO-2-1_S21_L003_R1_001.fastq.gz,UE-D60-BO-2-1_S21_L003_R2_001.fastq.gz \ --outFileNamePrefix HBO \ --outSAMtype BAM SortedByCoordinate \ --outSAMunmapped Within \ --outSAMattributes Standard \ --readFilesCommand gunzip -c |
主要问题:
- 默认不支持gz格式的fastq,加一行命令即可,记得带-c,否则原文件就被解压了。
- 支持多对fastq文件,格式就是fq1_1,fq1_2 fq2_1,fq2_2
二次比对
用于cufflinks和stringtie的比对
待续~
标签:
工具
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· .NET Core GC压缩(compact_phase)底层原理浅谈
· 现代计算机视觉入门之:什么是图片特征编码
· .NET 9 new features-C#13新的锁类型和语义
· Linux系统下SQL Server数据库镜像配置全流程详解
· 现代计算机视觉入门之:什么是视频
· 【译】我们最喜欢的2024年的 Visual Studio 新功能
· 个人数据保全计划:从印象笔记迁移到joplin
· Vue3.5常用特性整理
· 重拾 SSH:从基础到安全加固
· 为什么UNIX使用init进程启动其他进程?
2017-03-27 Pymonitor使用方法和常见问题