gffread icon indicating copy to clipboard operation
gffread copied to clipboard

extract cds sequence +- given bases

Open vivekbhr opened this issue 5 years ago • 1 comments

Hi @gpertea

I want to use gffread to get CDS sequences +- N bases. Something like --w-add <N> but for CDS --x-add <N>. Would it be easy to implement this?

Best, Vivek

vivekbhr avatar Mar 19 '20 11:03 vivekbhr

@vivekbhr waiting this feature request to be addressed you can use agat_sp_extract_sequences.pl from AGAT running:

agat_sp_extract_sequences.pl --gff inputgff --fasta input.fa --up <N> --down <N>

P.S: AGAT and gffread are complementary

Juke34 avatar Apr 09 '21 16:04 Juke34