美文网首页
小麦Hisat2建库及遇到的问题解决

小麦Hisat2建库及遇到的问题解决

作者: zjl_8c97 | 来源:发表于2020-01-07 19:22 被阅读0次

一、建库

extract_splice_sites.py Triticum_aestivum.IWGSC.45.gtf > Triticum_aestivum.IWGSC.45.ss
extract_exons.py Triticum_aestivum.IWGSC.45.gtf  > Triticum_aestivum.IWGSC.45.exon
hisat2-build --ss Triticum_aestivum.IWGSC.45.ss --exon Triticum_aestivum.IWGSC.45.exon genome.fa genome_tran
##利用hisat-build构建小麦转录组比对的索引文件,--ss 可变剪切文件,--exon 外显子文件,后接参考基因组,index文件前缀为genome_tran

hisat官网:http://ccb.jhu.edu/software/hisat2/index.shtml
具体建库方法可参考:

(1)https://www.jianshu.com/p/071c1757ded1
(2)https://www.jianshu.com/p/426b5fefe1a5

二、遇到的问题

Segmentation fault (core dumped) /mnt/sda/software/hisat2-2.0.5/hisat2-build --ss Triticum_aestivum.IWGSC.45.ss --exon Triticum_aestivum.IWGSC.45.exon ../Triticum_aestivum.IWGSC.dna.toplevel.fa genome_tran

解决
通过dmesg命令可以查看发生段错误的程序名称、引起段错误发生的内存地址、指令指针地址、堆栈指针地址、错误代码、错误原因等。
我这里发现是由于该盘下内存不够导致,后换内存大的目录即运行成功。
更多段错误原因可参考:https://blog.csdn.net/YSBJ123/article/details/50035169

相关文章

网友评论

      本文标题:小麦Hisat2建库及遇到的问题解决

      本文链接:https://www.haomeiwen.com/subject/yzddactx.html