Printing a sequence from a fasta file

后端 未结 4 1467
太阳男子
太阳男子 2021-01-13 20:45

I often need to find a particular sequence in a fasta file and print it. For those who don\'t know, fasta is a text file format for biological sequences (DNA, proteins, etc.

4条回答
  •  深忆病人
    2021-01-13 21:38

    Using the > as the record separator:

    awk -v seq="sequence2" -v RS='>' '$1 == seq {print RS $0}' file
    
    >sequence2
    ACTGACTGACTGACTG
    ACTGACTGACTGACTG
    

提交回复
热议问题