本文整理了Java中org.broadinstitute.gatk.utils.sam.GATKSAMRecord.getFirstOfPairFlag()
方法的一些代码示例,展示了GATKSAMRecord.getFirstOfPairFlag()
的具体用法。这些代码示例主要来源于Github
/Stackoverflow
/Maven
等平台,是从一些精选项目中提取出来的代码,具有较强的参考意义,能在一定程度帮忙到你。GATKSAMRecord.getFirstOfPairFlag()
方法的具体详情如下:
包路径:org.broadinstitute.gatk.utils.sam.GATKSAMRecord
类名称:GATKSAMRecord
方法名:getFirstOfPairFlag
暂无
代码示例来源:origin: broadgsa/gatk-protected
@Override
public int compare(final ReadCost o1, final ReadCost o2) {
final String s1 = o1.read.getReadName() + (o1.read.getReadPairedFlag() ? (o1.read.getFirstOfPairFlag() ? "/1" : "/2") : "");
final String s2 = o2.read.getReadName() + (o2.read.getReadPairedFlag() ? (o2.read.getFirstOfPairFlag() ? "/1" : "/2") : "");
return s1.compareTo(s2);
}
};
代码示例来源:origin: broadgsa/gatk-protected
continue; // read is non-informative or MQ0
if (a.getAlleleIfInformative().equals(refAllele, true) && el.getKey().getReadPairedFlag()) {
if (el.getKey().getReadNegativeStrandFlag() == el.getKey().getFirstOfPairFlag())
REF_F2R1++;
else
if (el.getKey().getReadNegativeStrandFlag() == el.getKey().getFirstOfPairFlag())
ALT_F2R1++;
else
代码示例来源:origin: broadgsa/gatk
public Integer map(RefMetaDataTracker tracker, ReferenceContext ref, AlignmentContext context) {
for ( final PileupElement p : context.getBasePileup() ) {
final GATKSAMRecord read = p.getRead();
final int offset = p.getOffset();
final boolean firstOfPair = ! read.getReadPairedFlag() || read.getFirstOfPairFlag();
if ( firstOfPair && read.getMappingQuality() >= MIN_MAPPING_QUAL && p.getQual() >= MIN_BASE_QUAL ) {
final byte readBase = p.getBase();
final byte refBase = ref.getBase();
final int cycle = offset;
if ( BaseUtils.isRegularBase(readBase) && BaseUtils.isRegularBase(refBase) ) {
final TableKey key = new TableKey(read.getReadGroup().getReadGroupId(), cycle);
if ( ! table.containsRowID(key) ) {
table.set(key, "cycle", cycle);
table.set(key, "readgroup", read.getReadGroup().getReadGroupId());
table.set(key, "counts", 0);
table.set(key, "mismatches", 0);
}
table.increment(key, "counts");
if (readBase != refBase)
table.increment(key, "mismatches");
}
}
}
return null;
}
内容来源于网络,如有侵权,请联系作者删除!