ReSeq
ReSeq copied to clipboard
Reference fasta file extension
Hi,
First, thank you for ReSeq, which is very helpful.
While trying to simulate some reads with the following command:
reseq illuminaPE --maxFragLen 500 -b reference.bam --adapterFile adapters.fa -R referencesimu.fas -r reference.fasta -1 reference_R1.fastq -2 reference_R2.fastq --numReads 1240000
I just had the following error that I found difficult to understand:
!!! bool reseq::Reference::ReadFasta(const char*) [/usr/local/ReSeq/reseq/Reference.cpp:765]: error: Could not open referencesimu.fas for reading.
It seems that the version of the library "seqan" you use does not support the ".fas" extension (it now does).
It is easy to get around the issue, but may be the error message could be more informative.
Thanks,
Frederic