This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to extract multiple bits of information that appear on different lines within the same text file

I am trying to extract the sequence ID and cluster number that occur on different lines within the same text file.

The input looks like

    >Cluster 72
    0   319aa, >O311_01007... *
    >Cluster 73
    0   318aa, >1494_00753... *
    1   318aa, >1621_00002... at 99.69%
    2   318aa, >1622_00575... at 99.37%
    3   318aa, >1633_00422... at 99.37%
    4   318aa, >O136_00307... at 99.69%
    >Cluster 74
    0   318aa, >O139_01028... *
    1   318aa, >O142_00961... at 99.69%
    >Cluster 75
    0   318aa, >O300_00856... *

The desired output is the sequence ID in one column and the corresponding cluster number in the second.

>O311_01007  72
>1494_00753  73
>1621_00002  73
>1622_00575  73
>1633_00422  73
>O136_00307  73
>O139_01028  74
>O142_00961  74
>O300_00856  75

Can anyone help with this?

command line coding text processing
• 1,850 views • link
written by timothy.ghaly • timothy.ghaly

Kloetzl, you're an absolute legend. That worked perfectly.

Thank you very much.

• 1 views • link
written by timothy.ghaly • timothy.ghaly

1 answer 

cat foo.txt | awk '/^>/{n=$2}!/^>/{printf "%.11s %d\n", $3, n}'
• 1 views • link
written by kloetzl kloetzl

Log in to answer this question.

More posts like this