建立refFlat.txt文件有问题怎么解决

这篇文章主要讲解了“建立refFlat.txt文件有问题怎么解决”，文中的讲解内容简单清晰，易于学习与理解，下面请大家跟着小编的思路慢慢深入，一起来研究和学习“建立refFlat.txt文件有问题怎么解决”吧！

转录组索引

统计基因区域覆盖比例命令：

/share/work/biosoft/java/jre1.8.0_73/bin/java -jar /share/work/biosoft/picard/picard-tools-2.1.0/picard.jar CollectRnaSeqMetrics REF_FLAT=/share/work/database/ref/Saccharomyces_cerevisiae/Yeastgenome/refFlat.txt INPUT=/share/nas1/renzx/project/zx-20211019-394-1_ref_15jiaomu/ano/2.hisat2_stringtie/1.Mapping/3ts-30min/3ts-30min.bam OUTPUT=/share/nas1/renzx/project/zx-20211019-394-1_ref_15jiaomu/ano/2.hisat2_stringtie/2.Map_stat/3ts-30min.map_base.stat CHART_OUTPUT=/share/nas1/renzx/project/zx-20211019-394-1_ref_15jiaomu/ano/2.hisat2_stringtie/2.Map_stat/3ts-30min.Normalized_Distance_Along_Transcript.pdf STRAND_SPECIFICITY=NONE VALIDATION_STRINGENCY=LENIENT

refFlat.txt 文件行数过少导致出错：

查看原因：

生成该文件报错，导致生成的基因区域过少

gtfToGenePred -genePredExt saccharomyces_cerevisiae.gtf refFlat.tmp.txt
no exons defined for group HRA1, feature noncoding_exon (perhaps try -ignoreGroupsWithoutExons)

改进：

gtfToGenePred -genePredExt -ignoreGroupsWithoutExons saccharomyces_cerevisiae.gtf refFlat.tmp.txt

感谢各位的阅读，以上就是“建立refFlat.txt文件有问题怎么解决”的内容了，经过本文的学习后，相信大家对建立refFlat.txt文件有问题怎么解决这一问题有了更深刻的体会，具体使用情况还需要大家实践验证。这里是亿速云，小编将为大家推送更多相关知识点的文章，欢迎关注！