This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Formatting problem gene table conversion

Hey all, a question about a problem in formatting. I have this table, i would like to filter it based on the genes that have matches (if they have, the matches are present in the lines below them until the next gene).

>gene1
match1
match2
match3
>gene2
>gene3
match1
match2
match3
>gene4
>gene5
>gene6

Desired output:

>gene1 match1
>gene1 match2
>gene1 match3
>gene3 match1
>gene3 match2
>gene3 match3
awk bash python
• 1,569 views • link
updated by ReWeeda ReWeeda • written by lessismore lessismore

lessismore : Please accept (green check mark) answers for this and your past questions (you can select multiple answers) to validate them.

• 0 views • link
written by GenoMax GenoMax

2 answers 

awk '/^>/{G=$0;next;}{printf("%s\t%s\n",G,$0);}'  input.txt
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Alternatively, if you are not familiar with awk, you can use a simple script like the following one written in python:

with open('file_name_or_path_here','r') as file_:
for line in file_:
     line=line.rstrip()
     if line[0]=='>':
          gene = line
     else:
          print ('\t'.join([gene,line])+'\n')
• 0 views • link
written by ReWeeda ReWeeda

Log in to answer this question.

More posts like this