码迷,mamicode.com
首页 > 其他好文 > 详细

【BioCode】Elm格式中提取位点信息

时间:2017-07-04 18:31:51      阅读:175      评论:0      收藏:0      [点我收藏+]

标签:iter   exp   sap   close   main   写入文件   osi   .com   block   

说明:

①Elm格式:

PLMD ID    Uniprot Accession    Position        Type        Sequence         Species          PMIDs
PlMD编号 Uniprot数据库编号         位点    翻译后修饰类型    序列信息        物种          PMID

PLMD-1      O00115        52        Ubiquitination    MIPLLLAALLCVPAGALTC Homo sapiens    21963094;23266961

技术分享


②代码说明:从上述格式中提取Position信息,当 Uniprot Accession相同时,Position位于同一行。不相同时,回车换行。每个位点之间用空格隔开。

代码:

 

package single;

import java.io.BufferedReader;
import java.io.BufferedWriter;
import java.io.FileReader;
import java.io.FileWriter;

public class Elm_site {
    // 从Elm中获得位点 每一行是一个蛋白质的的位点
    public void get_site(String path) {
        try {
            FileReader reader = new FileReader(path);
            BufferedReader br = new BufferedReader(reader);
            String str=null;
            String now="O00115";
            FileWriter fileWritter = new FileWriter("E:\\experiment--help\\linglingbao\\site.txt");
            BufferedWriter bufferWritter = new BufferedWriter(fileWritter);

            while((str=br.readLine())!=null){
                String[] temp =str.split("    ");//每一列之间使用TAB隔开的,需要用TAB来截取
                String uni_id=temp[1];
                if(uni_id.equals(now)){
                    System.out.print(temp[2]+" ");//输出
                    bufferWritter.write(temp[2]+" ");//写入文件
                    bufferWritter.flush();
                }else{
                    now=uni_id;
                    System.out.print("\n"+temp[2]+" ");
                    bufferWritter.write("\n"+temp[2]+" ");
                    bufferWritter.flush();
                }
            }
            br.close();
            reader.close();
        } catch (Exception e) {
            // TODO Auto-generated catch block
            e.printStackTrace();
        }
    }
    public static void main(String[] args) {
        Elm_site site = new Elm_site();
        String path = "E:\\experiment--help\\linglingbao\\ubiquitination.txt";
        site.get_site(path);
    }
}

 

 输出结果格式:

技术分享

 

【BioCode】Elm格式中提取位点信息

标签:iter   exp   sap   close   main   写入文件   osi   .com   block   

原文地址:http://www.cnblogs.com/yumiaomiao/p/7117350.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!