Search Results

Perfect thank you.It works fine

---------- Post updated at 02:33 PM ---------- Previous update was at 02:30 PM ----------


strange I checked it again here is the result...

thank you for your help,but I got this out put output
Transcript Name GO
POPTR_0016s06290.1 ||98654
POPTR_2158s00200.1 ||11324
POPTR_0004s22390.1 ||12897...

I have file which contains gene lines something like this

Transcript Name GO
POPTR_0016s06290.1 98654
POPTR_2158s00200.1 11324
POPTR_0004s22390.1 12897...

perfect,Thank you it works fine

---------- Post updated at 08:56 AM ---------- Previous update was at 08:55 AM ----------

Great It works nicely,I meant "justlooks" perl code it gave me *.3 out...

hi justlooks your code output contains

>POPTR_0005s00710.3|PACid:18208427
ACAAAACAAAAACAAAAGAGAAAGAGAGAAAGAAAGCACTTTCTTTTGTGTATTGTATTGGAAGGCAGCTTGTTGTGTTA...

but your code is not working POPTR_1855s00201 this file contains set of different values.If I use it just return following

POPTR_1855s00200.1|PACid:18205976...

It works perfectly,Thank you.I have one more question.
I have fatsa file something like this

>POPTR_1446s00200.20|PACid:18205974...

Thank you it works perfectly,your code is faster.I tried with following code earlier

awk '/mRNA/{split($(NF-2),a,"=");print;next}{sub(/t=.*/,"t="a[2],$1)}1 ' FS=\; OFS=\; input > tmp1
awk...

your out put missing with .1 I need to replace POPTR_0001s00200.1,which is name= value with PAC ids.

---------- Post updated at 08:26 AM ---------- Previous update was at 08:22 AM ----------

Hi...

Thanks rdcwayx you are really talented,I just want to remove spaces between ; because your output gave following Parent=POPTR_0001s00230.1; PACid=18235601

---------- Post updated at 04:04 AM...

Here is the expected input and output,but your code not generate following output

sscaffold_1 phytozome6 gene 12632 13612 . + . ID=POPTR_0001s00200;Name=POPTR_0001s00200...

Yes It supports FS,Thank you for your help but It I want unique sub sets depend on the POPTR_ for example,those are gene names

scaffold_997 phytozome6 gene 1687 2351 . - ...

yes It is different
scaffold_1 phytozome6 gene 12632 13612 . + . ID=POPTR_0001s00200;Name=POPTR_0001s00200
scaffold_1 phytozome6 mRNA 12632 13612 ...

Thank you for your reply but I get following

scaffold_1 phytozome6 gene 2330052 2335284 . - . ID=.1;Name=PT_0001s02940
scaffold_1 phytozome6 mRNA ...

I have used following script but still I cant replace the Parent=PAC

awk '{if(substr($9,17,10)=="Name=PT") n=substr($9,22,16) gsub("Parent=PAC:"substr($9,47,9),"Parent="n".1");
...

I have following text

scaffold_1 phytozome6 gene 12632 13612 . + . ID=PT_0001s00200;Name=PT_0001s00200
scaffold_1 phytozome6 mRNA 12632 13612 ....