文本格式转换

mswsg

文本格式如下

1. ID          Ta.3
   
2. TITLE       Beta-amylase
   
3. GENE        amy1
   
4. GENE_ID     543318
   
5. LOCUSLINK   543318
   
6. SCOUNT      407
   
7. SEQUENCE    ACC=JP820949.1; NID=g359554132; SEQTYPE=mRNA
   
8. SEQUENCE    ACC=JV869389.1; NID=g391051274; SEQTYPE=mRNA
   
9. //
   
10. ID          Ta.68263
    
11. TITLE       Transcribed locus, weakly similar to XP_003570811.1 PREDICTED: protein SPA1-RELATED 3-like [Brachypodium distachyon]
    
12. SCOUNT      2
    
13. SEQUENCE    ACC=AK334954.1; NID=g241987697; SEQTYPE=HTC
    
14. SEQUENCE    ACC=DR734341.1; NID=g70961145; LID=16675; SEQTYPE=EST
    
15. //
    
16. ID          Ta.68264
    
17. TITLE       TSA: Triticum aestivum Unigene64601 mRNA sequence
    
18. SCOUNT      3
    
19. SEQUENCE    ACC=AK334170.1; NID=g241986909; SEQTYPE=HTC
    
20. SEQUENCE    ACC=DR734784.1; NID=g70961588; LID=16675; SEQTYPE=EST
    
21. SEQUENCE    ACC=GAEF01020581.1; NID=g478393045; SEQTYPE=mRNA
    

复制代码

结果如下

1. Ta.3 JP820949.1
   
2. Ta.3 JV869389.1
   
3. Ta.68263 AK334954.1
   
4. Ta.68263 DR734341.1
   
5. Ta.68264 AK334170.1
   
6. Ta.68264 DR734784.1
   
7. Ta.68264 GAEF01058385.1
   

复制代码

ly5066113

回复 1# mswsg


try:

1. awk -F'[ =;]+' '/^ID/{s=$2}/^SEQUENCE/{print s,$3}' file

复制代码

mswsg

十分感谢！回复 2# ly5066113


   

elu_ligao

1. awk -F '[[:space:]=;]+' '/^ID/{id=$2}/^SEQUENCE/{print id,$3}'

复制代码

rdcwayx

完全用shell写的

1. #! /usr/bin/env bash
   
2. 
   
3. while read -r A B
   
4. do
   
5.     [[ "$A" =~ ^ID ]] && id="$B"
   
6.     [[ "$A" =~ ^SEQUENCE ]] && [[ "$B" =~ ACC=([^;]*) ]] && echo "$id" ${BASH_REMATCH[1]}
   
7. done < file
   

复制代码

1. $ bash a.sh
   
2. 
   
3. Ta.3 JP820949.1
   
4. Ta.3 JV869389.1
   
5. Ta.68263 AK334954.1
   
6. Ta.68263 DR734341.1
   
7. Ta.68264 AK334170.1
   
8. Ta.68264 DR734784.1
   
9. Ta.68264 GAEF01020581.1
   

复制代码

blackold

回复 2# ly5066113


    Tim 来抢红包。

zhaopingzi

Tim回复的太快了，1秒啊

ly5066113

回复 6# blackold


男人嘛，下手就是要狠一点。

laliheyi

回复 8# ly5066113


    请问你是如何做到秒杀的

ly5066113

回复 9# laliheyi


你没发现所有的悬赏帖都是秒杀的么？
不是我秒杀了楼主，是楼主秒杀了我。