perl程序实现对fasta序列反向互补

最新推荐文章于 2022-11-27 20:29:36 发布

小余吐泡泡

最新推荐文章于 2022-11-27 20:29:36 发布

阅读量1k

点赞数

分类专栏： Perl

本文链接：https://blog.csdn.net/ncyujijian/article/details/107620732

版权

本文介绍如何使用Perl编程语言实现对fasta格式的DNA序列进行反向互补转换。通过理解DNA双链结构，我们可以创建算法来计算序列的反向互补，这对于生物信息学分析至关重要。

摘要由CSDN通过智能技术生成

#!/usr/bin/perl
use strict;
use warnings;


open IN,"$ARGV[0]" or die "Can't open fasta file!";
while (<IN>) {
   
	chomp;
	if (/>\S+/)

最低0.47元/天解锁文章

确定要放弃本次机会？

福利倒计时

: :

立减 ¥

普通VIP年卡可用

立即使用

小余吐泡泡

关注关注

0
点赞
踩
3

收藏

觉得还不错? 一键收藏
0
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫

专栏目录

生物信息脚本练习（2）求反向互补序列

07-30

3106

这是一个简单的生信脚本练习，对一段碱基序列切片，然后求反向互补序列。注意其中join 函数的用法#求反向互补序列 lll = [] with open("Program1.txt","r") as f: seq = f.readlines() seq = seq[1] print(seq) print("\n") seq = seq[50:101] pr

perl实现序列反向互补

little^raccoon

03-23

3445

#!/usr/bin/perl -w print"Input sequence:\n"; chomp(my $seq = <STDIN>); $ seq =~ tr/atcguATCGU/tagcaTAGCA/; print “Reverse complement sequence:\n” print scalar reverse $seq; #或者使用下面这种形式输出，因为reve...

参与评论您还未登录，请先登录后发表或查看评论

反向互补序列

weixin_30315905的博客

12-10

588

#!/usr/bin/perl -w use strict; use Getopt::Long; use Bio::SeqIO; my %opts; GetOptions(\%opts, "i=s", "o=s"); if (!(defined $opts{i} && defined $opts{o})){ print "\n"; &...

fasta格式文件转换为反向互补链

weixin_30764771的博客

09-22

721

1 #!/usr/bin/env perl 2 3 use strict; 4 use warnings; 5 6 local $/ =">"; 7 while(<>){ 8 chomp; 9 my $lineSepPos = index($_," 10 "); 11 my $header = substr($_,0,$lineSep...

perl应用：ＤＮＡ互补序列的获取

GRC

10-24

5603

sub revcom { # A subroutine to compute the reverse complement of DNA sequence # 一个获取DNA互补序列的子程序 my($dna)=@_; print $dna."\n"; my ($revcom)=reverse($dna); $revcom=~tr/ACGTacgt/TGCAtgca/; return

[笔记]python对FASTA文件的处理

热门推荐

Cccrush的博客

04-27

2万+

这学期选了生信的选修课—perl/python在生物信息学中的应用把结课作业的代码整理出来主要是python对FASTA文件的读取和数据处理FASTA文件数据处理FASTA文件读取：只含一个基因序列将FASTA文件的基因序列读取到一个列表中，列表中的每个元素为每一行基因序列构成的字符串f=open('/home/miaoyr/perl_practice/test1_file/DTNBP1.fast...

linux提取基因名称和序列,一种批量提取基因组基因信息并翻译比对分析序列的方法与流程...

weixin_39631094的博客

05-14

2902

技术特征：1.一种批量提取基因组基因信息并翻译比对分析序列的方法，其特征在于，将某一物种的转录本id或者基因id，依据供试基因组cds文件、蛋白质文件、gff文件和染色体fasta文件信息，通过6个perl脚本程序，实现目标转录本或基因在基因组上的位置、长度、正反义链结构信息的提取，并在染色体fasta文件上提取该转录本或基因的cds或基因序列，在基因组蛋白文件上提取该转录本的蛋白质序列；最后对所...

Perl的一些初级练习题

qq_45091883的博客

05-15

3353

Perl形势不太好，我是必修课要学的，大家可以去学习Python，不过Perl处理文本还是很可以的。有错误欢迎私信我 1.读入一个文件，将序列变成一行输出 #!/usr/bin/perl -w use strict; use 5.026; open IN,"<$ARGV[0]",or die $!; open OUT,">out.txt",or die $!; while(<IN>){ chomp; print OUT $_; } print OUT "\n"; close

使用python读取和分析fasta文件

qq_18369669的博客

12-05

1万+

分享一些处理fasta文件的python函数

perl 统计字符串中特定字符出现的次数

hs6605015的博客

04-29

1949

perl 统计字符串中特定字符出现的次数方法一： #!/usr/bin/perl -w use strict; my $str = 'abcsfsaf#sdagasdga#sdgaghoiiopjh#dsfjkopahkl;fjh#dsjfklpafj '; my $count = $str =~ tr/#/#/; print $count; 统计了字符串$str 中出现的#的个数，在这里用到了转换操作符tr///，简单快速，当之首选。附上s/// 和tr/// 的区别： s///：替

perl碱基计数的三种不同的方法，for循环+while循环+tr///

GRC

09-26

5758

第一种方式，我们运用的for循环，然后用了if......else.......结构 for ($position=0;$position<length $DNA;++$position) { $base=substr($DNA,$position,1); if ($base eq 'A') { $count_A=$count_A+1; } elsif ($base eq 'T

perl实现计算GC含量

little^raccoon

03-22

2462

#!/usr/bin/perl -w unless(@ARGV==2){ die "Usage: perl $0 <input_fasta> <output> error:$!\n"; } my($input, $output)=@ARGV; open INPUT, $input; open OUTPUT, ">$output"; #===============...

perl中tr的用法

qq_42119386的博客

06-17

1822

perl中tr的用法 (2014-05-29 13:57:30) 转载▼ 标签：杂谈分类：Perl 找了不少关于tr的用法，发现最后面附的那个帖子算是比较详细的。个人测试的结果总结如下：（1）/c表示把匹配不上的字符进行替换. $temp="AAAABCDEF"; $count=$temp=~tr/A/H/c; print "$temp\t$count\n"; 结果：AAAAHHHHH 5 ...

[Python|生信]从Fasta文件出发获取序列的基本信息

最新发布

F_zero_T_hero的博客

11-27

2065

背景最近参加了个生信的面试，记录一下有意思的面试题。题目描述要求从提供的*.fasta文件出发：获得序列的反向互补序列，并统计信息：序列条数，碱基总数，N50，N90，GC 含量。提取每条序列上 32bp-332bp、780bp-992bp 的序列。统计单碱基重复 4 次及以上的序列在每条序列上出现的次数。如 AAAAA 或者 TTTT 等如 “AAATTTTTTTCCCCAAAAAAA”，结果如下： the 4th character T repeat 7 timesthe 11th ch

生信（二）反向互补序列

生信了

10-23

9051

关键词：reverse; complement; sequence; **如何得到一段基因序列的反向互补序列？**这是基因测序领域经常遇到的问题。其实答案很简单，许多现成的软件都有这个功能。但是作为一个生信人，当然可以自己实现一个了。首先想到的也是最基础的方法就是利用**多个if…else…**的语句进行判断选择。这种方法太笨拙，写出来的代码很不好看。今天我们分享几种好一点的方法。 Shell版...

**生信自学记录1——获取Fastq格式的反向互补序列**

weixin_43392168的博客

10-14

1537

生信自学记录1——获取Fastq格式的反向互补序列总共分为三步 1.读取基因序列的str格式，返回反向互补序列str 2.打开fastq格式的文本提取基因序列，返回互补序列list 3.读取互补序列的list替换原基因序列，返回反向互补后的fastq格式的文本具体代码如下： 1.读取基因序列的str格式，返回反向互补序列str def to_reverse(str):#获得目的序列的反向互补序...