- 论坛徽章:
- 5
|
本帖最后由 klainogn 于 2013-11-15 15:58 编辑
- awk '{match($6,/[0-9]+S/);len=substr($6,RSTART,RLENGTH)+0;str=$1"_"gsub("-","",$9)+1}NR==FNR{gsub(/[0-9]*D|[0-9]*S86M|[0-9]*S/,"",$6);split($6,a,/[A-Z]/);for(i in a) start+=a[i];print str"\n"substr($10,start,len)"\n+\n"substr($11,start,len)"\n";start=0;next}{str2=$6;gsub(/[0-9]*I|[0-9]*S86M|[0-9]*S/,"",$6);split($6,a,/[A-Z]/);for(i in a) start+=a[i];print str,str2,$4+start;start=0}' file file
- FCD2FEHACXX:8:2301:1875:46038#79_118_1
- GT
- +
- B?
- FCD2FEHACXX:8:1104:6915:94857#79_118_2
- AG
- +
- ??
- FCD2FEHACXX:8:2307:11230:38678#79_118_2
- CCCTC
- +
- A@@>A
- FCD2FEHACXX:8:2301:1875:46038#79_118_1 16M1D70M2S 112452
- FCD2FEHACXX:8:1104:6915:94857#79_118_2 2S86M 73823
- FCD2FEHACXX:8:2307:11230:38678#79_118_2 44M4I13M4D8M1I13M5S 37911841
复制代码 |
|