This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problem using SED for change IDs in gff file

Hi,

I want to use SED to change the names ID in a GFF file. I just want to change "ILINNGPN_01481" that appears only once for "tfpI" in the whole file. I am using the simplest line:

sed 's/ILINNGPN_01481/tfpI/g' bovis_PROKKA.gff > bovis_PROKKA_edited.gff

The IDs changes happened well, but suddenly the CDS coords also change for no reason.

looking for ILINNGPN_01482 in the file (Before):

NODE_16_length_55139_cov_21.7549        Prodigal:2.6    CDS     25370   26155   .       -       0       ID=ILINNGPN_01482;inference=ab initio prediction:Prodigal:2.6;locus_tag=ILINNGPN_01482;product=hypothetical protein

looking for tfpI in the edited file (After):

NODE_16_length_55139_cov_21.7549        Prodigal:2.6    CDS     24581   25048   .       +       0       ID=tfpI;inference=ab initio prediction:Prodigal:2.6,similar to AA sequence:GCF_002014975.1_ASM201497v1_genomic.gbff:WP_078275419.1;locus_tag=tfpI;product=pilin

Have a clue for what's going on here?

Thanks!

sed gff
• 1,948 views • link
updated by Ram Ram • written by luciabilbao03 • luciabilbao03

what is the output of

cat -n bovis_PROKKA.gff | sed 's/ILINNGPN_01481/tfpI/g' | grep -n tfpI
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

1 answer

Sed command has ILINNGPN_01481

You are looking at ILINNGPN_01482 in the before file

• 0 views • link
written by rbagnall rbagnall

Thanks rbagnall! rigth then I posted this I realize that! That's what was happening

• 0 views • link
written by luciabilbao03 • luciabilbao03

Log in to answer this question.

More posts like this