org.broadinstitute.gatk.utils.sam.GATKSAMRecord.getFirstOfPairFlag()方法的使用及代码示例

x33g5p2x  于2022-01-20 转载在 其他  
字(2.1k)|赞(0)|评价(0)|浏览(63)

本文整理了Java中org.broadinstitute.gatk.utils.sam.GATKSAMRecord.getFirstOfPairFlag()方法的一些代码示例,展示了GATKSAMRecord.getFirstOfPairFlag()的具体用法。这些代码示例主要来源于Github/Stackoverflow/Maven等平台,是从一些精选项目中提取出来的代码,具有较强的参考意义,能在一定程度帮忙到你。GATKSAMRecord.getFirstOfPairFlag()方法的具体详情如下:
包路径:org.broadinstitute.gatk.utils.sam.GATKSAMRecord
类名称:GATKSAMRecord
方法名:getFirstOfPairFlag

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;
}

相关文章

微信公众号

最新文章

更多

GATKSAMRecord类方法