忘记密码   免费注册 查看新帖 | 论坛精华区

ChinaUnix.net

  平台 论坛 博客 认证专区 大话IT 视频 徽章 文库 沙龙 自测 下载 频道自动化运维 虚拟化 储存备份 C/C++ PHP MySQL 嵌入式 Linux系统
12
最近访问板块 发新帖
楼主: kevin_furant

求助小脚本编写 [复制链接]

论坛徽章:
0
发表于 2017-02-26 19:04 |显示全部楼层
回复 9# kevin_furant

楼上大神们的建议很不错,双端拆分的话应该是miseq / nextseq / miniseq平台的下机数据咯~

论坛徽章:
2
综合交流区版块每日发帖之星
日期:2016-07-06 06:20:00综合交流区版块每日发帖之星
日期:2016-08-16 06:20:00
发表于 2017-02-26 20:03 |显示全部楼层
提示: 作者被禁止或删除 内容自动屏蔽

论坛徽章:
0
发表于 2017-02-26 20:04 |显示全部楼层
中规中矩的写一个:
  1. #!/usr/bin/perl

  2. use warnings;
  3. use strict;

  4. open FILE_1, '<', 'lst1.txt' or die "Error 1 occured: $!";
  5. open FILE_2, '<', 'lst2.txt' or die "Error 2 occured: $!";

  6. my (@sample, @index, $sample, $index);
  7. while (<FILE_1>) {
  8.         chomp;
  9.         if (/^(\w+),(\w+)/m) {
  10.                 $sample = $1.",".$2;
  11.                 push @sample, $sample;
  12.         }       
  13. }

  14. while (<FILE_2>) {
  15.         chomp;
  16.         my @med = split;
  17.         $index = join ",", @med;
  18.         push @index, $index;
  19. }

  20. while (@sample) {
  21.         my $sig_smp = shift @sample;
  22.         my $sig_index = shift @index;
  23.         my $combine = $sig_smp.",".$sig_index;
  24.         print $combine, "\n";
  25. }
复制代码

论坛徽章:
1
15-16赛季CBA联赛之天津
日期:2017-03-21 18:26:14
发表于 2017-02-27 13:50 |显示全部楼层
谢谢,不过我自己写的也搞定了,
#! usr/bin/perl
@ARGV ||die "can not open the file!\n";
open (my $fh_1, "< $ARGV[0]");
open (my $fh_2, "< $ARGV[1]");
my %hash;
my @arr;
while (<$fh_1>){
         chomp;
        s/[ATGC]{6,8},,//;
        @arr = split /,/;       
chomp(my $lst2 = <$fh_2>);
        my ($idx) = grep s/\s+/,/, ($lst2);
push @arr, $idx;
$hash{$arr[0]} = [@arr];
print join(",",@{$hash{$arr[0]}}),"\n";
}
close $fh_1;
close $fh_2;
写得烂但是运行没发现问题。

论坛徽章:
1
15-16赛季CBA联赛之天津
日期:2017-03-21 18:26:14
发表于 2017-02-27 14:07 |显示全部楼层
回复 10# 华小飞_Perl
是illunina hiseq4000的数据

论坛徽章:
1
15-16赛季CBA联赛之天津
日期:2017-03-21 18:26:14
发表于 2017-02-27 14:09 |显示全部楼层
回复 7# moperyblue
是的,我第一次就是这么干的,因为在学perl,所以想找个东西练练手。

论坛徽章:
0
发表于 2017-02-27 20:29 |显示全部楼层
kevin_furant 发表于 2017-02-27 14:07
回复 10# 华小飞_Perl
是illunina hiseq4000的数据

原来这样 数据量有点大 要考虑脚本性能
来自安卓客户端来自客户端
您需要登录后才可以回帖 登录 | 注册

本版积分规则

  

北京皓辰网域网络信息技术有限公司. 版权所有 京ICP证:060528号 北京市公安局海淀分局网监中心备案编号:1101082001
广播电视节目制作经营许可证(京) 字第1234号 中国互联网协会会员  联系我们:
感谢所有关心和支持过ChinaUnix的朋友们 转载本站内容请注明原作者名及出处

清除 Cookies - ChinaUnix - Archiver - WAP - TOP