比对软件之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的比对

 

待续~

 

参考:比对软件STAR的简单使用

 

posted @   Life·Intelligence  阅读(3061)  评论(0编辑  收藏  举报
(评论功能已被禁用)
编辑推荐:
· .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使用方法和常见问题
TOP
点击右上角即可分享
微信分享提示