/* * Eoulsan development code * * This code may be freely distributed and modified under the * terms of the GNU Lesser General Public License version 2.1 or * later and CeCILL-C. This should be distributed with the code. * If you do not have a copy, see: * * http://www.gnu.org/licenses/lgpl-2.1.txt * http://www.cecill.info/licences/Licence_CeCILL-C_V1-en.txt * * Copyright for this code is held jointly by the Genomic platform * of the Institut de Biologie de l'École normale supérieure and * the individual authors. These should be listed in @author doc * comments. * * For more information on the Eoulsan project and its aims, * or to join the Eoulsan Google group, visit the home page * at: * * http://outils.genomique.biologie.ens.fr/eoulsan * */ package fr.ens.biologie.genomique.eoulsan.bio.readsfilters; import static fr.ens.biologie.genomique.eoulsan.util.StatUtils.mean; import fr.ens.biologie.genomique.eoulsan.EoulsanException; import fr.ens.biologie.genomique.eoulsan.bio.ReadSequence; /** * This class define a filter based on mean quality of a read. * @since 1.0 * @author Maria Bernard * @author Laurent Jourdren */ public class QualityReadFilter extends AbstractReadFilter { public static final String FILTER_NAME = "quality"; private double qualityThreshold = -1.0; @Override public boolean accept(final ReadSequence read) { if (read == null) { return false; } return mean(read.qualityScores()) > this.qualityThreshold; } @Override public String getName() { return "quality"; } @Override public String getDescription() { return "Quality threshold ReadFilter"; } @Override public void setParameter(final String key, final String value) throws EoulsanException { if (key == null || value == null) { return; } if ("threshold".equals(key.trim())) { try { this.qualityThreshold = Double.parseDouble(value.trim()); } catch (NumberFormatException e) { return; } if (this.qualityThreshold < 0.0) { throw new EoulsanException( "Invalid qualityThreshold: " + this.qualityThreshold); } } else { throw new EoulsanException( "Unknown parameter for " + getName() + " read filter: " + key); } } @Override public void init() { if (this.qualityThreshold < 0.0) { throw new IllegalArgumentException( "Quality threshold is not set for " + getName() + " read filter."); } } @Override public String toString() { return this.getClass().getSimpleName() + "{qualityThreshold=" + this.qualityThreshold + "}"; } // // Constructor // /** * Public constructor. */ public QualityReadFilter() { } }