This is a test version of Biostars. For the public version, visit https://www.biostars.org.
easy way to fetch out the surrounding CDS of a pair of given coordinates

Hi all,

I have a bacteria genome and the corresponding GFF annotation. Now I have a pair of coordinates and am wondering whether there is a quick way to fish out the surrounding CDS from the GFF file. Thanks in advance.

gff
• 1,310 views • link
written by zhangdengwei zhangdengwei

I finally succumb to coding by myself. The biopython module and pandas package can be helpful when doing such a thing.

• 0 views • link
written by zhangdengwei zhangdengwei

1 answer

agat_sp_filter_record_by_coordinates.pl from AGAT for example.
Then you catch CDS with an awk command: awk '{if($3=="CDS"){print $0}}' my.gff

• 0 views • link
written by Juke34 Juke34

Log in to answer this question.

More posts like this