Package net.sf.picard.reference

Examples of net.sf.picard.reference.ReferenceSequence


    ReferenceSequenceFile ref = ReferenceSequenceFileFactory.getReferenceSequenceFile(new File(args[1]));
    SAMRecordIterator iterator = reader.iterator();
    String refName = iterator.next().getReferenceName();
    iterator.close();
    ReferenceSequence sequence = ref.getSequence(refName);

    iterator = reader.iterator();
    Sam2CramRecordFactory factory = new Sam2CramRecordFactory(sequence.getBases());
    factory.preserveReadNames = false;
    factory.captureAllTags = false;
    factory.captureUnmappedBases = true;
    factory.captureUnmappedScores = true;
    factory.preserveReadNames = true;
View Full Code Here


    if (container.sequenceId == SAMRecord.NO_ALIGNMENT_REFERENCE_INDEX) {
      refs = new byte[] {};
    } else if (prevSeqId < 0 || prevSeqId != container.sequenceId) {
      SAMSequenceRecord sequence = cramHeader.samFileHeader
          .getSequence(container.sequenceId);
      ReferenceSequence referenceSequence = Utils
          .trySequenceNameVariants(referenceSequenceFile,
              sequence.getSequenceName());
      refs = referenceSequence.getBases();
      prevSeqId = container.sequenceId;
    }

    long time1 = System.nanoTime();
View Full Code Here

TOP

Related Classes of net.sf.picard.reference.ReferenceSequence

Copyright © 2018 www.massapicom. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.