Hello,
I have a question regarding the RNA-star alignment program. I'm trying to map thousands of reads (~7000 sequences) against a fasta file containing a handful of sequences (~20 sequences). In order to do that in RNA-star, I had to generate a genome index based on the small fasta file that contains ~20 sequences:
$STAR --runMode genomeGenerate --genomeDir $direc --genomeFastaFiles $filename --runThreadN 4
Next I run the following command to do the mapping:
$STAR --genomeLoad NoSharedMemory --outFileNamePrefix $direc/STAR. --genomeDir $direc --readFilesIn file.fasta
After running the above I get the following segmentation fault error:
Jul 22 11:44:49 ..... Started STAR run
Jul 22 11:44:50 ..... Started mapping
Segmentation fault (core dumped)
I think this has something to do with the size of the genome index I'm generating. This is because I was able to run the above commands successfully when I had a larger number of sequences used to generate the genome index. Has anyone run into this issue before? If yes, any ideas/suggestions on how to resolve it?