召唤大神!有偿求助!

楼主

xinchao [离线]

2★★☆☆☆

发帖数:752 积分2295
1楼

生物信息相关的,CDS fasta格式下,标题序列太长,想把它改短,locos到+-号那一节不要,
求帮忙处理!会的可加我QQ 546169782 或微信 lei961740594 联系我,有偿求助!帮我修改文件内容........
>Cotton_D_gene_10030635 locus=scaffold2:766445:770682:-
ATGGAACCGAGAGAAGATGAACGGTCACTGATCAGTTTCAGTTCACCTTT
AGAAGATAGCAAAGAGAAAGATCTAGGAAAGGAAGAGATCATCGAAGAAG
TGAAAAAACAGCTATGGCTGGCAGCTCCATTAATATGTGAAAGCTTGTTG
CTGTATTGCTTGCAATTGATCTCTGTTATGTTTGTGGGTCATTTGGGTGC




有个两个perl脚本参考
#!/usr/bin/perl


$filename = "F:\基因库";


open IN, "$filename";


while(<IN>){
if(~/^>/){
s/locus=.*$//;
}
print;
}

http://muchong.com/html/201507/9166800.html



2016/2/4 18:37:24

星座神 [离线]

2★★☆☆☆

发帖数:759 积分2348
2楼

很急!很关键!


2016/2/4 20:17:04
返回本版
1

请您先 登录 再进行发帖

快速回复楼主