Package fork.lib.bio.seq

Examples of fork.lib.bio.seq.FastaSequenceExtractor


        File f= fs[i];
        if(FileName.getExt(f).equals("bed")){
            GenomicRegionsBuilder gb= new BedReader(f).getGenomicRegionsBuilder();
            GenomicRegionsModifier gm= new GenomicRegionsModifier(gb);
            gb= gm.getThreePrimeEnds(150, 0); /////////// modify this
            FastaSequenceExtractor fe= new FastaSequenceExtractor(fastaf, gb, null);

            File of= new File(od+"/"+tag+"_"+FileName.getBaseName(f)+".fasta");
            fe.writeToFastaFile(of);
        }
    }
   
   
   
View Full Code Here


public void writeToFile(File out, String tit) throws Exception{
    out.getParentFile().mkdirs();
    BufferedWriter bw= new BufferedWriter(new FileWriter(out));
    BedGraphExporter.writeTitle(bw, tit);
   
    FastaSequenceExtractor se= new FastaSequenceExtractor(gen, gbout);
    SequenceExtractorEntry en;
    while( (en= se.nextEntry())!=null ){
        GenomicRegion gr= en.getGenomicRegion();
        NucleotideSequence seq= NucleotideSequenceParser.parseSequence(en.getSequence());
        double at= seq.atContent();
        bw.write( gr.toBedFormat()+"\t"+ ((double)Math.round(at*100)/100)+"\n" );
    }
View Full Code Here

    }
   
   
protected void init() throws Exception{
    GenomicRegionsBuilder gb= new BedReader(f).getGenomicRegionsBuilder();
    FastaSequenceExtractor se= new FastaSequenceExtractor(sf, gb);
    SequenceExtractorEntry en;
    while((en=se.nextEntry())!=null){
        GenomicRegion reg= en.getGenomicRegion();
        String seq= en.getSequence();
        String id= reg.getID();
        idseq.put(id , seq);
        tab.appendEmptyRow(id);
View Full Code Here

TOP

Related Classes of fork.lib.bio.seq.FastaSequenceExtractor

Copyright © 2018 www.massapicom. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.