zoukankan      html  css  js  c++  java
  • RepBaseRepeatMaskerEdition下载 | RepeatMasker

    开源的生物信息世界居然有这么个需要注册才能下载的工具,开源世界不是怎么方便怎么来吗?

    这个注册真的麻烦,这里上传了一个可以使用的版本。

    RepBaseRepeatMaskerEdition-20170127.tar.gz

    想转成fasta可以用如下命令:

    awk -f awk.sh RMRBSeqs.embl > raw.fasta
    cat raw.fasta | cut -f1-9 -d' '| sed -e 's/ID   //g' | sed -e 's/repeatmasker;//g'| sed -e 's/ //g' > RepBase.fasta
    

     

    其中awk.sh文件内容为:

    /^ID/   {printf(">%s;",$0); next;}
    /^(PT|PA)/  {printf(" %s;",$0); next;}
    /^/// {printf("
    "); next;}
    /^    / {printf("
    %s",substr($0,5)); next;}
        {
        /* ignore default */
        }
    END   {
        printf("
    ");
        }
    

      

  • 相关阅读:
    9
    8
    7
    6
    5
    第四周
    作业14-数据库
    作业13-网络
    作业12-流与文件
    作业11-多线程
  • 原文地址:https://www.cnblogs.com/leezx/p/8654838.html
Copyright © 2011-2022 走看看