我有以下包含基因名称和外显子位置的字符串。
gene1
250405..250551
251490..251884
195620..195641
154254..155469
156319..156548
gene2
171224..171481
403914..403970
147436..147627
149077..149157
11635..12021
gene3
107657..107844
105642..106214
103531..103691
407044..407435
405691..405843
我想编写一个脚本,该脚本将列出所有基因名称以及输出应显示的位置
gene1_exon1 250405..250551
gene1_exon2 251490..251884
geneN_exonN X..Y
有什么想法可以做到吗?
我尝试了以下脚本
grep [a-z] temp3 > names
cat names |
while read line; do
cat temp3 |
while read position; do
if [ "$line" != "$position" ] ; then
echo "$line" ${position} >> names_and_positions
else
break
fi
done
done