Hello,
I am using BFAST for indexing color space data using the command:
bfast index -f folder/human_bfast.fa -m '11111' -w 5 -A 1 -n 8
If I try to index the fasta file with a few chromosomes, say till 17, the tool successfully creates .bif file. But if I use a bigger fasta file (ie. from chr 1 till 22) for indexing then the tool throws error:
I want to create solid color space index for human fasta for all chromosomes from 1 through 22, X,Y
Please note that the step just prior to indexing (ie. creating reference in color space using the same fasta) was always successful. It creates .fa.cs.brg file successfully.
Can anybody please help me with the problem ?
Thanks in advance,
Aswathy
I am using BFAST for indexing color space data using the command:
bfast index -f folder/human_bfast.fa -m '11111' -w 5 -A 1 -n 8
If I try to index the fasta file with a few chromosomes, say till 17, the tool successfully creates .bif file. But if I use a bigger fasta file (ie. from chr 1 till 22) for indexing then the tool throws error:
Currently on [contig,pos]:
[------22,------662000]
Sorting by thread...
100.000 percent complete************************************************************
In function "RGIndexMergeHelperFromDiskContig_8": Fatal Error[ReadFileError]. Message: Could not read in tmp lower.
The file stream error was:: Bad file descriptor
***** Exiting due to errors *****
************************************************************
[------22,------662000]
Sorting by thread...
100.000 percent complete************************************************************
In function "RGIndexMergeHelperFromDiskContig_8": Fatal Error[ReadFileError]. Message: Could not read in tmp lower.
The file stream error was:: Bad file descriptor
***** Exiting due to errors *****
************************************************************
Please note that the step just prior to indexing (ie. creating reference in color space using the same fasta) was always successful. It creates .fa.cs.brg file successfully.
Can anybody please help me with the problem ?
Thanks in advance,
Aswathy
Comment