cgat-developers / cgat-apps

cgat-apps repository
Other
33 stars 14 forks source link

[bed2fasta]: Use of `--use-strand` and `--output-mode=segments` together gives incorrect result. #83

Open IanSudbery opened 3 years ago

IanSudbery commented 3 years ago

If one tries to use --output-mode=segments together with the --use-strand switch, to get transcript sequences from bed12 files, the output is incorrect for negative strand transcripts. The sequences of the indevidual exons are corrects reverse complemented, but then they are concatenated in the wrong order.

I will try to fix this if I get time, shouldn't be hard, I'm just busy.