use of ubic.gemma.core.analysis.sequence.RepeatScan in project Gemma by PavlidisLab.
the class ArrayDesignRepeatScanCli method processArrayDesign.
private void processArrayDesign(ArrayDesign design) {
ArrayDesign thawed = this.thaw(design);
// no taxon is passed to this method so all sequences will be retrieved even for multi taxon arrays
Collection<BioSequence> sequences = ArrayDesignSequenceAlignmentServiceImpl.getSequences(thawed);
RepeatScan scanner = new RepeatScan();
Collection<BioSequence> altered;
if (this.inputFileName != null) {
altered = scanner.processRepeatMaskerOutput(sequences, inputFileName);
} else {
altered = scanner.repeatScan(sequences);
}
AbstractCLI.log.info("Saving...");
bsService.update(altered);
if (this.inputFileName != null) {
this.audit(thawed, "Repeat scan data from file: " + inputFileName + ", updated " + altered.size() + " sequences.");
} else {
this.audit(thawed, "Repeat scan done, updated " + altered.size() + " sequences.");
}
AbstractCLI.log.info("Done with " + thawed);
}
use of ubic.gemma.core.analysis.sequence.RepeatScan in project Gemma by PavlidisLab.
the class ArrayDesignRepeatScanTaskImpl method execute.
@Override
public TaskResult execute() {
ArrayDesign ad = taskCommand.getArrayDesign();
ad = arrayDesignService.thaw(ad);
Collection<BioSequence> sequences = ArrayDesignSequenceAlignmentServiceImpl.getSequences(ad);
RepeatScan scanner = new RepeatScan();
scanner.repeatScan(sequences);
return new TaskResult(taskCommand, new ModelAndView(new RedirectView("/", true)));
}
Aggregations