我按照标准流程进行转录组分析:
hisat2 -p 8 --no-mixed --dta -x $indexs -1 $sample1 -2 $sample2 -S $out
samtools sort -@ 8 $d'.bam' $d'.sorted'
samtools view -bS $d'.sam' > $d'.bam'
stringtie -p 8 -G duckbase.refseq.v4.NCBI2016.gtf -A $d'.tab' -o $d'gtf' $d'.sorted.bam'
为什么参考基因组中有20000多个基因,而进行stringtie -A输出的基因表达量tab文件中只有10000多个基因呢,是软件版本的问题还是我参数设置的不对呢?
转录组分析流程中表达基因的基因数量问题
- 写回答
- 好问题 0 提建议
- 追加酬金
- 关注问题
- 邀请回答
-
0条回答 默认 最新
悬赏问题
- ¥15 基于卷积神经网络的声纹识别
- ¥15 Python中的request,如何使用ssr节点,通过代理requests网页。本人在泰国,需要用大陆ip才能玩网页游戏,合法合规。
- ¥100 为什么这个恒流源电路不能恒流?
- ¥15 有偿求跨组件数据流路径图
- ¥15 写一个方法checkPerson,入参实体类Person,出参布尔值
- ¥15 我想咨询一下路面纹理三维点云数据处理的一些问题,上传的坐标文件里是怎么对无序点进行编号的,以及xy坐标在处理的时候是进行整体模型分片处理的吗
- ¥15 CSAPPattacklab
- ¥15 一直显示正在等待HID—ISP
- ¥15 Python turtle 画图
- ¥15 stm32开发clion时遇到的编译问题