search

wwood / dirseq

Work out whether RNAseq reads in general agree with the direction of the gene predicted
MIT License
3 stars 1 forks source link

Associate sequence ID with output #3

Closed TYMichaelsen closed 6 years ago

TYMichaelsen commented 6 years ago

hi Ben,

I am using your dirseq script, which really works nicely. But, I would like to be able to merge the output with other files I've generated using a identity that links directly to the ORF/gene sequence that is mapped. ATM I only have indirect ways of doing so, like the contig, position and strand.

To be specific, I am working with output from Prokka, meaning the 9th column in the .gff file has a ID and/or locus_tag with a unique ORF identity. If this could be included in the output as a column somehow I would be very happy!

Hope my question is clear.

Best, Thomas

wwood commented 6 years ago

Hi Thomas,

I've defaulted it to "ID", but you can modify it yourself from there. Implemented in 41cda02. If you like it I'll bump the release - let me know. ta

TYMichaelsen commented 6 years ago

Hi Ben, Works like a charm, I like that feature!

wwood commented 6 years ago

Cheers.