建立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文件有問題怎么解決這一問題有了更深刻的體會，具體使用情況還需要大家實踐驗證。這里是億速云，小編將為大家推送更多相關知識點的文章，歡迎關注！