SAMComparator.java

/*
 *                  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.io.comparators;

import java.io.BufferedReader;
import java.io.IOException;
import java.io.InputStream;
import java.io.InputStreamReader;
import java.util.Collection;
import java.util.HashSet;
import java.util.Set;

import com.google.common.collect.Sets;

import fr.ens.biologie.genomique.eoulsan.Globals;
import fr.ens.biologie.genomique.eoulsan.util.EnhancedBloomFilter;

/**
 * This class allow compare two SAM file with use BloomFilter.
 * @since 2.0
 * @author Sandrine Perrin
 */
public class SAMComparator extends AbstractComparatorWithBloomFilter {

  public static final String COMPARATOR_NAME = "SAMComparator";
  private static final Collection<String> EXTENSIONS = Sets.newHashSet(".sam");

  final Set<String> tagsToNotCompare;

  private int numberElementsCompared;

  @Override
  public boolean compareFiles(final EnhancedBloomFilter filter,
      final InputStream is) throws IOException {

    final BufferedReader reader =
        new BufferedReader(new InputStreamReader(is, Globals.DEFAULT_CHARSET));
    String line = null;
    this.numberElementsCompared = 0;

    while ((line = reader.readLine()) != null) {
      this.numberElementsCompared++;

      // Header
      if (line.charAt(0) == '@') {

        // Skip specified tag in header sam file
        if (!this.tagsToNotCompare.contains(getTag(line))) {

          if (!filter.mightContain(line)) {
            // Save line occurs fail comparison
            setCauseFailComparison(line);

            reader.close();
            return false;
          }
        }
      } else {
        // Line
        if (!filter.mightContain(line)) {
          // Save line occurs fail comparison
          setCauseFailComparison(line);

          reader.close();
          return false;
        }
      }
    }
    reader.close();

    // Check count element is the same between two files
    if (this.numberElementsCompared != filter.getAddedNumberOfElements()) {
      setCauseFailComparison("Different count elements "
          + this.numberElementsCompared + " was "
          + filter.getAddedNumberOfElements() + " expected.");
      return false;
    }

    return true;
  }

  //
  // Other methods
  //

  private static String getTag(final String samHeaderLine) {

    if (samHeaderLine.length() == 0) {
      return "";
    }

    final int pos = samHeaderLine.indexOf('\t');

    if (pos == -1) {
      return samHeaderLine.substring(1);
    }

    return samHeaderLine.substring(1, pos);
  }

  @Override
  public String getName() {

    return COMPARATOR_NAME;
  }

  @Override
  public Collection<String> getExtensions() {
    return EXTENSIONS;
  }

  @Override
  public int getNumberElementsCompared() {
    return this.numberElementsCompared;
  }

  //
  // Constructor
  //

  /**
   * Public constructor public FastqComparator(final boolean useSerializeFile) {
   * super(useSerializeFile); } /** Public constructor
   * @param useSerializeFile true if it needed to save BloomFilter in file with
   *          extension '.ser'
   */
  public SAMComparator(final boolean useSerializeFile) {
    super(useSerializeFile);

    this.tagsToNotCompare = new HashSet<>();
  }

  /**
   * Public constructor, specify all headers tags not used to compare.
   * @param useSerializeFile true if it needed to save BloomFilter in file with
   *          extension '.ser'
   * @param headersTags all headers tags
   */
  public SAMComparator(final boolean useSerializeFile,
      final String... headersTags) {
    super(useSerializeFile);

    if (headersTags == null) {
      throw new NullPointerException("headersTags is null");
    }

    this.tagsToNotCompare = Sets.newHashSet(headersTags);
  }

}