MCPcopy Create free account
hub / github.com/BirolLab/RNA-Bloom / assembleTranscriptsSE

Method assembleTranscriptsSE

src/rnabloom/RNABloom.java:5546–5590  ·  view source on GitHub ↗
(RNABloom assembler, boolean forceOverwrite,
            String outdir, String name, String txptNamePrefix,
            long sbfSize, int sbfNumHash, int numThreads, 
            String[] forwardReadPaths, String[] reverseReadPaths,
            int minTranscriptLength, boolean keepArtifact, boolean keepChimera,
            float minKmerCov, boolean reduceRedundancy, boolean writeUracil,
            boolean usePacBioPreset, String minimapOptions, boolean destroyBf)

Source from the content-addressed store, hash-verified

5544 }
5545
5546 private static void assembleTranscriptsSE(RNABloom assembler, boolean forceOverwrite,
5547 String outdir, String name, String txptNamePrefix,
5548 long sbfSize, int sbfNumHash, int numThreads,
5549 String[] forwardReadPaths, String[] reverseReadPaths,
5550 int minTranscriptLength, boolean keepArtifact, boolean keepChimera,
5551 float minKmerCov, boolean reduceRedundancy, boolean writeUracil,
5552 boolean usePacBioPreset, String minimapOptions, boolean destroyBf) throws IOException, InterruptedException {
5553
5554 final File txptsDoneStamp = new File(outdir + File.separator + STAMP_TRANSCRIPTS_DONE);
5555 final String transcriptsFasta = outdir + File.separator + name + ".transcripts" + FASTA_EXT;
5556 final String shortTranscriptsFasta = outdir + File.separator + name + ".transcripts.short" + FASTA_EXT;
5557
5558 if (forceOverwrite || !txptsDoneStamp.exists()) {
5559 Timer timer = new Timer();
5560 assembler.setupKmerScreeningBloomFilter(sbfSize, sbfNumHash);
5561
5562 assembler.assembleSingleEndReads(forwardReadPaths,
5563 reverseReadPaths,
5564 transcriptsFasta,
5565 shortTranscriptsFasta,
5566 numThreads,
5567 minTranscriptLength,
5568 keepArtifact,
5569 keepChimera,
5570 txptNamePrefix,
5571 minKmerCov,
5572 writeUracil);
5573 System.out.println("Transcripts assembled in " + timer.elapsedDHMS());
5574
5575 touch(txptsDoneStamp);
5576
5577 System.out.println("Assembled transcripts at `" + transcriptsFasta + "`");
5578 }
5579 else {
5580 System.out.println("WARNING: Transcripts were already assembled for \"" + name + "\"!");
5581 }
5582
5583 if (destroyBf) {
5584 assembler.destroyAllBf();
5585 }
5586
5587 if (reduceRedundancy) {
5588 assembleTranscriptsNR(assembler, outdir, name, forceOverwrite, numThreads, keepArtifact, usePacBioPreset, minimapOptions);
5589 }
5590 }
5591
5592 private static void assembleTranscriptsPE(RNABloom assembler, boolean forceOverwrite,
5593 String outdir, String name, String txptNamePrefix,

Callers 1

mainMethod · 0.95

Calls 6

elapsedDHMSMethod · 0.95
assembleTranscriptsNRMethod · 0.95
touchMethod · 0.80
destroyAllBfMethod · 0.80

Tested by

no test coverage detected