I test run your sample LongReads.fa
➜ HECIL git:(master) ✗ python HECIL.py -l LongRead.fa -s ShortRead.fq -len 202 -o Out
Started running alignment
[samopen] SAM header is present: 55 sequences.
-e Finished running alignment
Started running correction
-e Finished running correction
-e Generating output file with all (corrected and uncorrected) reads
-e Finished running HECIL
I notice the outfile "Corrected_LongRead.fa" is not formatted correctly, specially first sequence is missing the fasta header (attached). Did I missed something?
I test run your sample LongReads.fa ➜ HECIL git:(master) ✗ python HECIL.py -l LongRead.fa -s ShortRead.fq -len 202 -o Out Started running alignment [samopen] SAM header is present: 55 sequences. -e Finished running alignment
Started running correction -e Finished running correction
-e Generating output file with all (corrected and uncorrected) reads
-e Finished running HECIL
I notice the outfile "Corrected_LongRead.fa" is not formatted correctly, specially first sequence is missing the fasta header (attached). Did I missed something?
Corrected_LongRead.fa.txt Out.txt