/*
* Copyright 2015 MiLaboratory.com
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
package com.milaboratory.core.io.sequence;
import com.milaboratory.core.io.sequence.fastq.QualityFormat;
import com.milaboratory.core.sequence.NSequenceWithQuality;
import com.milaboratory.core.sequence.UnsafeFactory;
import com.milaboratory.util.SingleIterator;
import java.nio.charset.Charset;
import java.util.Iterator;
/**
* @author Dmitry Bolotin
* @author Stanislav Poslavsky
*/
public abstract class SingleReadLazy implements SingleRead {
final long id;
final byte[] buffer;
final int descriptionFrom;
final short sequenceOffset, qualityOffset, dataLength, descriptionLength;
final boolean replaceWildcards;
NSequenceWithQuality sequenceWithQuality;
String description;
private SingleReadLazy(long id,
byte[] buffer,
int descriptionFrom,
short sequenceOffset,
short qualityOffset,
short dataLength,
short descriptionLength,
boolean replaceWildcards) {
this.id = id;
this.buffer = buffer;
this.descriptionFrom = descriptionFrom;
this.sequenceOffset = sequenceOffset;
this.qualityOffset = qualityOffset;
this.dataLength = dataLength;
this.descriptionLength = descriptionLength;
this.replaceWildcards = replaceWildcards;
}
abstract byte getQualityOffset();
@Override
public String getDescription() {
if (description == null)
description = createDescription();
return description;
}
@Override
public NSequenceWithQuality getData() {
if (sequenceWithQuality == null)
sequenceWithQuality = createNSequenceWithQuality();
return sequenceWithQuality;
}
private String createDescription() {
return new String(buffer, descriptionFrom, descriptionLength, Charset.defaultCharset());
}
private NSequenceWithQuality createNSequenceWithQuality() {
try {
return UnsafeFactory.fastqParse(buffer, descriptionFrom + sequenceOffset,
descriptionFrom + qualityOffset, dataLength, getQualityOffset(), id, replaceWildcards);
} catch (Exception e) {
throw new RuntimeException("Error while parsing read:\n" +
new String(buffer, descriptionFrom - 1, qualityOffset + dataLength), e);
}
}
@Override
public int numberOfReads() {
return 1;
}
@Override
public SingleRead getRead(int i) {
if (i != 0)
throw new IndexOutOfBoundsException();
return this;
}
@Override
public long getId() {
return id;
}
@Override
public Iterator<SingleRead> iterator() {
return new SingleIterator<>((SingleRead) this);
}
public static SingleReadLazy create(final QualityFormat format,
long id,
byte[] buffer,
int descriptionFrom,
short dataOffset,
short qualityOffset,
short dataLength,
short descriptionLength,
boolean replaceWildcards) {
if (format == QualityFormat.Phred33)
return new SingleReadLazy(id, buffer, descriptionFrom, dataOffset, qualityOffset, dataLength,
descriptionLength, replaceWildcards) {
@Override
byte getQualityOffset() {
return 33;
}
};
else if (format == QualityFormat.Phred64)
return new SingleReadLazy(id, buffer, descriptionFrom, dataOffset, qualityOffset, dataLength,
descriptionLength, replaceWildcards) {
@Override
byte getQualityOffset() {
return 64;
}
};
throw new IllegalArgumentException("Unknown quality format.");
}
}