Bug: some extracted sequences are empty
Files: extract_subregions_from_bed_files.sh
?
descr:
$ cat t.bed
madrigal|scaffold_105785 22581 19141
$ ../extract_subregions_from_bed_files.sh -b t.bed
Extract from madrigal (scaffold_105785: 22581-19141).
WARN: seq_id 'scaffold_105785' => end '19141' is lower than start '22581'.
strand will be '-'.End of extract_subregions_from_bed_files.sh.
$ cat t.fa
>madrigal|scaffold_105785|19141-22581|-
A
The output sequence is wrong !