筛选染色体连续区域的策略

要求：筛选染色体上符合要求的连续的纯合突变区域。perl (/(ㄒoㄒ)/~~现在只会用perl)

笨办法，正则匹配：

举例子1：

my $info=(HOMA,871334,HET,874950,HOMA,876499,HOMA,887643);

while($info =~ /((HOMA,(\d+),){2,})/g){
print "$1\n"; ##$1是输出大括号内的内容，所有符合要求的具有连续连续至少两次纯合突变的区域被筛选出来啦
}

举例子2：

因为是按照染色体的位置顺序依次判断，预先设置一个变量来判断“区域中断点”，先设置为空，将读进去的区域第一个符合要求的变量赋给他，不同即为断点，判断之前存储的区域是否符合要求即可，然后从断点后继续筛选，将符合要求的纯合突变赋值给“断点”。。。

my @getpos=();##存储连续的homa位置信息
my $lasttype="";

foreach my $pos(@sortpos){##排好顺序的位置
my $geno=$$hash{$chr}{$pos}{'geno'};
if($geno ne $lasttype){
if($lasttype eq ""){
if($geno eq 'HOMA'){
push @getpos,$pos;
$lasttype='HOMA';
next;
}
}else{
check();##子函数检查模块
$lasttype="";
@getpos=();
}
}else{
push @getpos,$pos;
}
}

bingo!!!

转载于:https://www.cnblogs.com/gui-/p/6873906.html

筛选染色体连续区域的策略相关推荐

清除图片下默认的小间隙_PowerMILL软件应用策略(一）：模型区域清除策略
⦁:开粗参数 (主页图片就是一粗) 策略:模型区域清除策略参数: 0.1-毛胚:(有毛胚有刀路,限制刀路加工范围) 电极---以模型放大2MM 钢料---通常放大5MM 0.2-刀具: 电极:通常用 ...
SQL实现筛选出连续3天登录用户与窗口函数
阅读目录还原试题 SQL窗口函数一.窗口函数有什么用二.什么是窗口函数三.如何使用 1.专用窗口函数rank 2.其他专业窗口函数 3.聚合函数作为窗口函数 4.注意事项四.总结解题思路 ...
Excel·VBA自定义函数筛选单元格区域重复值
贴吧提问<哪位大神知道要怎么实现?>,Excel内置函数使用比较麻烦,VBA字典实现比较直观自定义函数UNIQUE_IF筛选单元格区域中的值,可以选择返回其中的唯一值或重复值,并用分隔符 ...
华为（huawei）USG6000的CLI命令行综合配置之安全区域/nat策略/安全策略配置
文章目录 2.USG6000在园区网中的相关cli配置(安全区域/nat策略/安全策略) 三.总结接上篇USG综合配置的博文. 本篇主要讲USG6000的安全区域/nat策略/安全策略的配置拓扑图 ...
pandas 更改单元格的值_懂Excel轻松入门Python数据分析包pandas(二十四)：连续区域...
此系列文章收录在公众号中:数据大宇宙 > 数据处理 >E-pd 转发本文并私信我"python",即可获得Python资料以及更多系列文章(持续更新的) 经常听别人说 ...
二值图片连续区域检测
https://www.cnblogs.com/ronny/p/img_aly_01.html 检测2行是否有重叠区域,即保证下一行的开头小于上一行的结尾同时下一行的结尾大于上一行的开头,2条件同时满 ...
Smart3D,CCC软件制作三维模型快速筛选出重点区域或大面积项目分割建模方法
现在是实景三维模型建模的时代,相信大家会经常用到Samrt3D进行建模,也会经常遇到把测区飞完了,想先挑选一部分的重点区域检查一下航片的质量,或者大面积项目建模需要.下面我就介绍一下在Smart3D中 ...
python 筛选提取连续多行_没有比这更简单的Python入门案例，用python打印你的宠物小精灵...
现在学习Python的人越来越多,但你会以怎样的方式开始编写你的第一个Python代码呢? 下面我教你一个简单的入门案例,很适合最开始的Python学习哦! 即将用到的Python基础语法是:prin ...
【操作系统】连续内存分配策略
最为简单的内存分配方法之一就是将内存分为多个固定大小的分区.每个分区只能容纳一个进程. 为此,操作系统有一个表用于记录哪些内存可用和哪些内存已用. 一开始,所有内存都可用于用户进程,因此能够作为一个大 ...

筛选染色体连续区域的策略

筛选染色体连续区域的策略相关推荐

最新文章

热门文章