(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)
| 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, |
no test coverage detected