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

Method assembleSingleEndReads

src/rnabloom/RNABloom.java:4844–4884  ·  view source on GitHub ↗
(String[] forwardReadPaths,
                                    String[] reverseReadPaths,
                                    String outFasta,
                                    String outFastaShort,
                                    int numThreads,
                                    int minTranscriptLength,
                                    boolean keepArtifact,
                                    boolean keepChimera,
                                    String txptNamePrefix,
                                    float minKmerCov,
                                    boolean writeUracil)

Source from the content-addressed store, hash-verified

4842 }
4843
4844 public void assembleSingleEndReads(String[] forwardReadPaths,
4845 String[] reverseReadPaths,
4846 String outFasta,
4847 String outFastaShort,
4848 int numThreads,
4849 int minTranscriptLength,
4850 boolean keepArtifact,
4851 boolean keepChimera,
4852 String txptNamePrefix,
4853 float minKmerCov,
4854 boolean writeUracil) throws IOException, InterruptedException {
4855
4856 //boolean assemblePolya = minPolyATailLengthRequired > 0;
4857 /*@TODO support prioritized assembly of polya reads */
4858
4859 boolean includeNaiveExtensions = true;
4860 boolean extendBranchFreeFragmentsOnly = false;
4861
4862 FastaWriter fout = new FastaWriter(outFasta, false, writeUracil);
4863 FastaWriter foutShort = new FastaWriter(outFastaShort, false, writeUracil);
4864 TranscriptWriter writer = new TranscriptWriter(fout, foutShort, minTranscriptLength, maxTipLength, writeUracil);
4865
4866 SingleEndReadsIterator readsItr = new SingleEndReadsIterator(forwardReadPaths, reverseReadPaths);
4867
4868 TranscriptAssemblyWorker[] workers = new TranscriptAssemblyWorker[numThreads];
4869 Thread[] threads = new Thread[numThreads];
4870 for (int i=0; i<numThreads; ++i) {
4871 workers[i] = new TranscriptAssemblyWorker(readsItr, writer, includeNaiveExtensions, extendBranchFreeFragmentsOnly, keepArtifact, keepChimera, false, minKmerCov);
4872 threads[i] = new Thread(workers[i]);
4873 threads[i].start();
4874 }
4875
4876 for (Thread t : threads) {
4877 t.join();
4878 }
4879
4880 fout.close();
4881 foutShort.close();
4882
4883 System.out.println("Parsed " + NumberFormat.getInstance().format(readsItr.numReadsParsed) + " reads.");
4884 }
4885
4886 public void assembleTranscriptsMultiThreaded(FragmentPaths fragPaths,
4887 String outFasta,

Callers 1

assembleTranscriptsSEMethod · 0.80

Calls 3

closeMethod · 0.95
joinMethod · 0.80
startMethod · 0.65

Tested by

no test coverage detected