怎样去掉不满足条件的块

小猪Gina

实际运行中，我贴出的那个程序有一个致命的bug，就是质量值中也会出现“@”，所以，得改。。。

mcshell

回复 9# 小猪Gina

1. my ($data,$x,$m)=do{local $/;<DATA>};
   
2. map{/(^\@\w+\s+([ATGC]+)\s+\+\s+(?{$m=length($2);$x = qr{\w{$m}}})(??{$x})$)(?{print "$1\n"})/smg}$data;   
   
3. 
   
4. __DATA__
   
5. @A
   
6. ATCG
   
7. +
   
8. aaaa
   
9. @B
   
10. ATCG
    
11. +
    
12. aaaa
    
13. ATCG
    
14. +
    
15. aaaa
    
16. @c
    
17. ATCG
    
18. +
    
19. aaaab

复制代码

最后一个@c 我加了一个b  不会匹配到。长度不够。

jingwu110

回复 3# Perlvim
DNA数据结构？在具体点

   

小猪Gina

回复 12# mcshell


    嗯嗯~谢谢~

小猪Gina

回复 13# jingwu110


    DNA的数据结构？我不知道您想知道什么。它有好几种不同的文件类型。比如.fq文件，比如.gff3文件。不同的文件类型它的结构也不一样。。嗯，我通常是碰到一个自己不认识的文件格式就会上网去搜一下。也木有做过总结。不过您要是拿具体的文件格式来问我，也许我可以帮您回答。。。