@@ -152,7 +152,11 @@ public Output performAnalysisPerSampleRemote(Readset rs, File inputBam, Referenc
152152
153153 wrapper .execute (inputBam , referenceGenome .getWorkingFastaFile (), siteTable , primerTable , genbank , metrics , blastDbs .get (referenceGenome .getGenomeId ()));
154154
155- output .addOutput (siteTable , "Tag-PCR Integration Sites" );
155+ if (siteTable .exists ())
156+ {
157+ output .addOutput (siteTable , "Tag-PCR Integration Sites" );
158+ }
159+
156160 if (designPrimers )
157161 {
158162 output .addOutput (primerTable , "Tag-PCR Primer Table" );
@@ -198,17 +202,24 @@ public Output performAnalysisPerSampleRemote(Readset rs, File inputBam, Referenc
198202 getPipelineCtx ().getLogger ().error ("Readset did not have information on input files" );
199203 }
200204
201- output .addSequenceOutput (siteTable ,
202- "Putative Integration Sites: " + rs .getName (),
203- "Tag-PCR Integration Sites" , rs .getReadsetId (),
204- null ,
205- referenceGenome .getGenomeId (),
206- "Reads: " + metricMap .get ("NumReadsSpanningJunction" ) +
207- "\n Junction hit rate (of alignments): " + pf .format (Double .parseDouble (metricMap .get ("PctReadsSpanningJunction" ))) +
208- (hitRate == null ? "" : "\n Junction hit rate (of total reads): " + pf .format (hitRate )) +
209- "\n Integration Sites: " + metricMap .get ("TotalIntegrationSitesOutput" ) +
210- "\n Alignments Matching Insert: " + metricMap .get ("FractionPrimaryAlignmentsMatchingInsert" )
211- );
205+ if (siteTable .exists ())
206+ {
207+ output .addSequenceOutput (siteTable ,
208+ "Putative Integration Sites: " + rs .getName (),
209+ "Tag-PCR Integration Sites" , rs .getReadsetId (),
210+ null ,
211+ referenceGenome .getGenomeId (),
212+ "Reads: " + metricMap .get ("NumReadsSpanningJunction" ) +
213+ "\n Junction hit rate (of alignments): " + pf .format (Double .parseDouble (metricMap .get ("PctReadsSpanningJunction" ))) +
214+ (hitRate == null ? "" : "\n Junction hit rate (of total reads): " + pf .format (hitRate )) +
215+ "\n Integration Sites: " + metricMap .get ("TotalIntegrationSitesOutput" ) +
216+ "\n Alignments Matching Insert: " + metricMap .get ("FractionPrimaryAlignmentsMatchingInsert" )
217+ );
218+ }
219+ else
220+ {
221+ getPipelineCtx ().getLogger ().info ("Site output not found" );
222+ }
212223
213224 return output ;
214225 }
0 commit comments