Bioinformatics Stack Exchange is a question and answer site for researchers, developers, students, teachers, and end users interested in bioinformatics. It only takes a minute to sign up.
Sign up to join this community
SAMtools returned an error message for my gzipped genome FASTA:
Indexed block-compressed FASTA file cannot be handled
The source code for the message is here.
// check if the it is a valid block-compressed file
if (IOUtil.isBlockCompressed(path, true)) {
throw new SAMException("Indexed block-compressed FASTA file cannot be handled: " + path);
}
What's a "block-compressed" file?
A block compression usually refers to compressing your file into a series of small blocks (with a tool like bgzip). This allows indexing in that the index can record which record lives in which block so that the whole file does not need to be decompressed and therefore can be accessed randomly. This isn't true of files compressed as a single block using gzip.
As far as I was aware samtools only works with uncompressed FASTA files.
