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

Method assembleFragments

src/rnabloom/RNABloom.java:5407–5439  ·  view source on GitHub ↗
(RNABloom assembler, boolean forceOverwrite,
            String outdir, String name, FastxFilePair[] fqPairs,
            String[] forwardReadPaths, String[] reverseReadPaths,
            long sbfSize, long pkbfSize, int sbfNumHash, int pkbfNumHash, int numThreads,
            int bound, int minOverlap, int sampleSize, int maxErrCorrItr, boolean extendFragments,
            int minKmerCoverage, boolean keepArtifact)

Source from the content-addressed store, hash-verified

5405 }
5406
5407 private static void assembleFragments(RNABloom assembler, boolean forceOverwrite,
5408 String outdir, String name, FastxFilePair[] fqPairs,
5409 String[] forwardReadPaths, String[] reverseReadPaths,
5410 long sbfSize, long pkbfSize, int sbfNumHash, int pkbfNumHash, int numThreads,
5411 int bound, int minOverlap, int sampleSize, int maxErrCorrItr, boolean extendFragments,
5412 int minKmerCoverage, boolean keepArtifact) throws FileFormatException, IOException, InterruptedException {
5413
5414 final File fragsDoneStamp = new File(outdir + File.separator + STAMP_FRAGMENTS_DONE);
5415
5416 if (forceOverwrite || !fragsDoneStamp.exists()) {
5417 FragmentPaths fragPaths = new FragmentPaths(outdir, name);
5418 fragPaths.deleteAll();
5419
5420 assembler.setupKmerScreeningBloomFilter(sbfSize, sbfNumHash);
5421 assembler.setupFragmentPairedKmersBloomFilter(pkbfSize, pkbfNumHash);
5422
5423 Quartiles fragStats = assembler.assembleFragmentsMultiThreaded(fqPairs,
5424 forwardReadPaths, reverseReadPaths, fragPaths,
5425 bound, minOverlap, sampleSize, numThreads,
5426 maxErrCorrItr, extendFragments, minKmerCoverage, keepArtifact);
5427
5428 String fragStatsFile = outdir + File.separator + name + ".fragstats";
5429 String graphFile = outdir + File.separator + name + ".graph";
5430
5431 assembler.updateGraphDesc(new File(graphFile));
5432 assembler.writeQuartilesToFile(fragStats, fragStatsFile);
5433
5434 touch(fragsDoneStamp);
5435 }
5436 else {
5437 System.out.println("WARNING: Fragments were already assembled for \"" + name + "!");
5438 }
5439 }
5440
5441 private static long[] splitFastaByLength(String inFasta, String outLongFasta, String outShortFasta, int lengthThreshold) throws IOException {
5442 FastaReader fin = new FastaReader(inFasta);

Callers 1

mainMethod · 0.95

Calls 7

deleteAllMethod · 0.95
updateGraphDescMethod · 0.80
writeQuartilesToFileMethod · 0.80
touchMethod · 0.80

Tested by

no test coverage detected