Java Reference類代碼示例

import com.google.api.services.genomics.model.Reference; //導入依賴的package包/類
private Map<String, Reference> createDummyReferenceMap(String[] contig) {
    Map<String,Reference> fake = new LinkedHashMap<>();
    for (String s : contig) {
        Reference r = new Reference();
        String id = "id-"+s;
        r.setName(s);
        r.setId(id);
        r.setLength(100L);
        fake.put(s, r);
    }
    return fake;
}
 

import com.google.api.services.genomics.model.Reference; //導入依賴的package包/類
/**
 * Collect a list of references mentioned in this Readgroupset and get their meta data.
 * @throws GeneralSecurityException 
 * @throws IOException 
 */
private Map<String, Reference> getReferences(ReadGroupSet readGroupSet) 
    throws IOException, GeneralSecurityException {
  Set<String> referenceSetIds = Sets.newHashSet();
  if (!Strings.isNullOrEmpty(readGroupSet.getReferenceSetId())) {
    LOG.info("Found reference set from read group set: " + 
        readGroupSet.getReferenceSetId());
    referenceSetIds.add(readGroupSet.getReferenceSetId());
  }
  if (readGroupSet.getReadGroups() != null) {
    LOG.info("Found read groups");
    for (ReadGroup readGroup : readGroupSet.getReadGroups()) {
      if (!Strings.isNullOrEmpty(readGroup.getReferenceSetId())) {
        LOG.info("Found reference set from read group: " + 
            readGroup.getReferenceSetId());
        referenceSetIds.add(readGroup.getReferenceSetId());
      }
    }
  }
  
  Map<String, Reference> references = Maps.newHashMap();
  for (String referenceSetId : referenceSetIds) {
    LOG.info("Getting reference set " + referenceSetId);
    ReferenceSet referenceSet = getApi().referencesets().get(referenceSetId).execute();
    if (referenceSet == null || referenceSet.getReferenceIds() == null) {
      continue;
    }
    for (String referenceId : referenceSet.getReferenceIds()) {
      LOG.fine("Getting reference  " + referenceId);
      Reference reference = getApi().references().get(referenceId).execute();
      if (!Strings.isNullOrEmpty(reference.getName())) {
        references.put(reference.getName(), reference);
        LOG.fine("Adding reference  " + reference.getName());
      }
    }
  }
  return references;
}
 

import com.google.api.services.genomics.model.Reference; //導入依賴的package包/類
public static final SAMRecord makeSAMRecord(Read read,
                                            ReadGroupSet readGroupSet, List<Reference> references,
                                            boolean forceSetMatePositionToThisPosition) {
  return makeSAMRecord(read, makeSAMFileHeader(readGroupSet, references));
}
 

import com.google.api.services.genomics.model.Reference; //導入依賴的package包/類
@Override Iterable<Reference> getResponses(SearchReferencesResponse response) {
  return response.getReferences();
}
 

import com.google.api.services.genomics.model.Reference; //導入依賴的package包/類
public ReadIteratorResource(ReadGroupSet readGroupSet, List<Reference> references,
    UnmappedReads<Read> unmappedReads, 
    Iterable<Read> iterable) {
  super(readGroupSet, references, unmappedReads, iterable, 
      new GenomicsConverter());
}
 

import com.google.api.services.genomics.model.Reference; //導入依賴的package包/類
@Override
public ReadIteratorResource getReads(String readsetId, 
    String sequenceName, int sequenceStart, int sequenceEnd) 
        throws IOException, GeneralSecurityException {
  LOG.info("Getting readset " + readsetId + ", sequence " + sequenceName + 
      ", start=" + sequenceStart + ", end=" + sequenceEnd);    
  final Genomics stub = getApi();

  try {
    ReadGroupSet readGroupSet = stub.readgroupsets().get(readsetId).execute();
    String datasetId = readGroupSet.getDatasetId();
    LOG.info("Found readset " + readsetId + ", dataset " + datasetId);
    // Fix up PP in Programs
    
     if (readGroupSet.getReadGroups() != null) { 
      for (ReadGroup readGroup : readGroupSet.getReadGroups()) {
        if (readGroup.getPrograms() != null) {
          for (Program PG : readGroup.getPrograms()) {
            if (PG.getPrevProgramId() != null && PG.getPrevProgramId().length() == 0) {
              PG.setPrevProgramId(null);
            }
          }
        }
      }   
    }
    
    final Map<String, Reference> references = getReferences(readGroupSet);
    final Reference reference = references.get(sequenceName);
    if (reference != null) {
        LOG.info("Reference for sequence name " + sequenceName + " is found, length="
            + String.valueOf(reference.getLength()));
    } else {
      LOG.warning("Reference for sequence name " + sequenceName + " not found");
    }
    LOG.info("Searching for reads in sequence " + sequenceName + " " + 
        String.valueOf(sequenceStart) + "-" + String.valueOf(sequenceEnd));
    com.google.cloud.genomics.gatk.common.UnmappedReads<Read> unmappedReads = null;
    if (sequenceName.isEmpty()) {
      unmappedReads = getUnmappedMatesOfMappedReads(readsetId); 
    }
    Paginator.Reads searchReads = Paginator.Reads.create(stub, ShardBoundary.Requirement.OVERLAPS);
    SearchReadsRequest readRequest = new SearchReadsRequest()
      .setReadGroupSetIds(Arrays.asList(readsetId))
      .setReferenceName(sequenceName)
      .setPageSize(2048);
    if (sequenceStart != 0) {
      readRequest.setStart(Long.valueOf(sequenceStart));
    }
    if (sequenceEnd != 0) {
      readRequest.setEnd(Long.valueOf(sequenceEnd));
    }
    Iterable<Read> reads = searchReads.search(readRequest); 
    
    return new ReadIteratorResource(readGroupSet, 
        Lists.newArrayList(references.values()), unmappedReads, reads);
  } catch (GoogleJsonResponseException ex) {
    LOG.warning("Genomics API call failure: " + ex.getMessage());
    if (ex.getDetails() == null) {
      throw ex;
    }
    throw new IOException(ex.getDetails().getMessage());
  }
}
 

import com.google.api.services.genomics.model.Reference; //導入依賴的package包/類
/**
 * Generates a SAMFileHeader from a ReadGroupSet and Reference metadata
 */
@Override
public SAMFileHeader makeSAMFileHeader(ReadGroupSet readGroupSet,
    List<Reference> references) {
  return ReadUtils.makeSAMFileHeader(readGroupSet, references);
}
 

import com.google.api.services.genomics.model.Reference; //導入依賴的package包/類
/**
 * Gets the references for a given referenceSetId using the Genomics API.
 *
 * @param referenceSetId The id of the referenceSet to query.
 * @param auth The OfflineAuth for the API request.
 * @return The list of references in the referenceSet.
 * @throws IOException
 */
public static Iterable<Reference> getReferences(String referenceSetId, OfflineAuth auth)
    throws IOException {
  Genomics genomics = GenomicsFactory.builder().build().fromOfflineAuth(auth);
  return Paginator.References.create(
      genomics).search(new SearchReferencesRequest().setReferenceSetId(referenceSetId));
}