我对RNA-seq领域真的很陌生,正试图将我的数据与STAR的病毒基因组进行比对。然而,我经常收到下面提到的错误。我试图更改节点数量(10(、线程数量(48(,并有大约200G的内存来运行我的作业,但这似乎没有帮助。。。有什么建议吗?
代码:12月24日19:09:40。。。。。启动STAR运行
12月24日19:09:40。。。。。加载基因组
12月24日19:09:41。。。。。已启动映射分段故障(堆芯转储(’
谢谢!
您是如何创建STAR索引的?我这么问是因为它取决于基因组的大小(我想你的基因组很小,因为它是病毒(,也许你使用的是默认值,而不是最佳值。
根据手册说明:
对于小基因组,参数
--genomeSAindexNbases
必须按比例缩小最小值(14,log2(基因组长度(/2-1(。例如,对于1兆碱基的基因组,这是相等的到9,对于100千碱基基因组,这等于7。
STAR --runThreadN 16
--runMode genomeGenerate
--genomeDir </path/to/your/starIndex/directory>
--genomeFastaFiles </your/reference/file>
--genomeSAindexNbases <your_numerical_value>