Question: STAR error (segmentation fault) in Aligning reads
0
gravatar for spriyansh29
9 months ago by
spriyansh2930
spriyansh2930 wrote:

I was running STAR on my local computer with i5 8 generation processor and 8 gigs of RAM. I was implementing a variant calling pipeline on simulated data of the ebola virus. I started with forward reads (.fq) and reverse reads (.fq) along with the reference genome of ebola(.fasta). I successfully used GATK to make reference genome index. While I was doing alignment with STAR (STAR --genomeDir ../genomedir/ --readFilesIn ../ebola_mut_reads1.fq ../ebola_mut_reads2.fq --runThreadN 4), I encountered an error.

Error is as follows

Dec 17 15:00:05 ..... started STAR run

Dec 17 15:00:05 ..... loading genome

Dec 17 15:00:06 ..... started mapping

Segmentation fault (core dumped)

star variant alignment ngs reads • 314 views
ADD COMMENTlink modified 4 months ago • written 9 months ago by spriyansh2930
1

Did you create the ebola index with the STAR aligner? You should have used something like STAR --mode genomeGenerate --genomeDir /path/to/genomeDir --genomeFastaFiles /path/to/ebola.fasta ... .

ADD REPLYlink written 9 months ago by michael.ante3.6k

Yes, I did, I used it along with GATK and samtools as follows,

samtools faidx ebola_ref.fasta

gatk CreateSequenceDictionary --REFERENCE ebola_ref.fasta --OUTPUT=ebola_ref.dict

STAR --runMode genomeGenerate --genomeDir ../genomedir/ --genomeFastaFiles --genomeSAindexNbases 6 ebola_ref.fasta --runThreadN 4

ADD REPLYlink written 9 months ago by spriyansh2930

Usually, STAR needs a lot of RAM to align and even more to build the index. Thus, having the STAR index files properly generated, you should be able to run the alignment.

Are there any clues in the log files from STAR, what went wrong? Have you tried it with R1 reads only?

ADD REPLYlink written 9 months ago by michael.ante3.6k

I have used with R1 and R2 reads together as they are forward/reverse reads. Also, where I could find the log for STAR, I have installed it using 'make STAR'

ADD REPLYlink written 9 months ago by spriyansh2930
0
gravatar for spriyansh29
4 months ago by
spriyansh2930
spriyansh2930 wrote:

It was a hardware constraint, as STAR requires a minimum of 30Gigs of RAM.

ADD COMMENTlink written 4 months ago by spriyansh2930

Certainly not for Ebola genome? Or did you also have human reads in your data?

ADD REPLYlink modified 4 months ago • written 4 months ago by genomax89k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1859 users visited in the last hour