召唤大神!有偿求助!
浏览量:772 回帖数:1
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
生物信息相关的,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