I am trying to use soap2.21 release to build and align reads to Human genome(UCSC hg18).
I downloaded individual chromosomes FASTA files and cat them into a single Human Genome file. The final file size is around 3GB.
Soap index built is run using the command:
./2bwt-builder humanGenome.fa
The above command generated 13 files. Now, I tried aligning some reads using this index file..
./soap –D index/humanGenome.fa.index –a reads/reads1_human.fastq –o output/reads_aligned_soap
The above command resulted in the following error on the console:
Begin Program SOAPaligner/soap2
Mon Feb 13 23:09:37 2012
Reference: index/humanGenome.fa.index
Query File a: reads/reads1_human.fastq
Output File: output/reads_aligned_soap
Load Index Table ...
Segmentation fault
Could someone throw some light on what could be the possible issue here?
Thanks and Regards,
gany
I downloaded individual chromosomes FASTA files and cat them into a single Human Genome file. The final file size is around 3GB.
Soap index built is run using the command:
./2bwt-builder humanGenome.fa
The above command generated 13 files. Now, I tried aligning some reads using this index file..
./soap –D index/humanGenome.fa.index –a reads/reads1_human.fastq –o output/reads_aligned_soap
The above command resulted in the following error on the console:
Begin Program SOAPaligner/soap2
Mon Feb 13 23:09:37 2012
Reference: index/humanGenome.fa.index
Query File a: reads/reads1_human.fastq
Output File: output/reads_aligned_soap
Load Index Table ...
Segmentation fault
Could someone throw some light on what could be the possible issue here?
Thanks and Regards,
gany
Comment