shell

biobaby

大家好，我有這樣一個文件file.文件信息如下：
chr1    23    25     +      ACT
chr2    11    15     +      ACTGA
chr5    4     10    -       GCTAGCT
現(xiàn)在想要將第5列分開，每一行一個信息：如下：

chr1    3     +      A
chr1    4     +      C
chr1    5     +      T
chr2    11    +     A
chr2    12    +     C
chr2    13    +     T
chr2    14    +     G
chr2    15    +     A
chr5     4      -     T
chr5     5      -     C
chr5     6      -     G
chr5     7      -     A
chr5     8      -     T
chr5     9      -     C
chr5     10      -     G
在輸出的過程中，如果第四列是+，則順序輸出，如果第四列是-，則第五列反向輸出
謝謝大家了

jason680

回復(fù) 1# biobaby

$ awk '{OFS="\t";t=split($5,a,"");for(n=1;n<=t;++n)print $1,$2+n-1,$4,$4=="+"?a[n]:a[t-n+1]}' file
chr1        23        +        A
chr1        24        +        C
chr1        25        +        T
chr2        11        +        A
chr2        12        +        C
chr2        13        +        T
chr2        14        +        G
chr2        15        +        A
chr5        4        -        T
chr5        5        -        C
chr5        6        -        G
chr5        7        -        A
chr5        8        -        T
chr5        9        -        C
chr5        10        -        G

   

LikeLx

1. awk '{for(i=0;i++<split($NF,b,"");)print $1,$2++,$4,b[i]}'

復(fù)制代碼

偷個懶！

我是一隻羊

寫完了跟二樓一樣...不發(fā)了..換個思路再想...