private void checkFile(String method) { if (!hasFileId()) { setFileId(""); // throw new IllegalArgumentException("Can not " + method + " without file."); } }
private VariantOverlappingStatus processVariantFile(Variant variant, Set<Integer> missingSamples, Put put, Integer fileId, VariantOverlappingStatus overlappingStatus, String gt) { LinkedHashMap<String, Integer> samplePosition = getSamplePosition(fileId); List<List<String>> samplesData = new ArrayList<>(samplePosition.size()); for (int i = 0; i < samplePosition.size(); i++) { samplesData.add(Collections.singletonList(gt)); } VariantBuilder builder = Variant.newBuilder( variant.getChromosome(), variant.getStart(), variant.getEnd(), variant.getReference(), variant.getAlternate()) .setStudyId(String.valueOf(studyConfiguration.getStudyId())) .setFormat("GT") .setFileId(fileId.toString()) .setSamplesPosition(samplePosition) .setSamplesData(samplesData); studyConverter.convert(builder.build(), put, missingSamples, overlappingStatus); return overlappingStatus; }
.addAlternate("<*>") .setStudyId(String.valueOf(studyConfiguration.getStudyId())) .setFileId(fileId.toString())
.setEnd(variantContext.getEnd()) .setStudyId(studyId) .setFileId(fileId);
.setAlternates(alternates) .setStudyId(studyId.toString()) .setFileId(fileId.toString()) .setFilter("PASS_" + fileId) .setQuality(fileId * 100.0)
public static Variant newVariant(String chromosome, int start, int end, String reference, String alternate, Integer fileId, Integer studyId) { int pad = (fileId - 1) * NUM_SAMPLES; return Variant.newBuilder() .setChromosome(chromosome) .setStart(start) .setEnd(end) .setReference(reference) .setAlternate(alternate) .setStudyId(studyId.toString()) .setFileId(fileId.toString()) .setFilter("PASS_" + fileId) .setQuality(fileId * 100.0) .addAttribute("AD", "1,2") .setFormat("GT", "DP", "GQX", "AD") .addSample("S" + (1 + pad), "./.", String.valueOf(11 + pad), "0.7", "1,2") .addSample("S" + (2 + pad), "1/1", String.valueOf(12 + pad), "0.7", "1,2") .addSample("S" + (3 + pad), "0/0", String.valueOf(13 + pad), "0.7", "1,2") .addSample("S" + (4 + pad), "1/0", String.valueOf(14 + pad), "0.7", "1,2").build(); }
@SuppressWarnings("unchecked") public static List<Variant> createFile2Variants(String chromosome, String fileId, String studyId) { List<Variant> variants = new LinkedList<>(); variants.add(Variant.newBuilder(chromosome, 1000, 1000, "A", "C") .setStudyId(studyId) .setFileId(fileId) .setFormat("GT", "DP", "GQX") .addSample("NA19600", "./.", "1", "0.7") .addSample("NA19660", "1/1", "2", "0.7") .addSample("NA19661", "0/0", "3", "0.7") .addSample("NA19685", "1/0", "4", "0.7") .build()); variants.add(Variant.newBuilder(chromosome, 1004, 1004, "A", "C") .setStudyId(studyId) .setFileId(fileId) .setFormat("GT", "DP", "GQX") .addSample("NA19600", "0/1", "1", "0.7") .addSample("NA19660", "0/0", "2", ".") .addSample("NA19661", "1/0", "3", "0.7") .addSample("NA19685", "0/0", "4", "..") .build()); return variants; }
variants.add(Variant.newBuilder(chromosome, 999, 999, "A", "C") .setStudyId(studyId) .setFileId(fileId) .setFormat("GT", "DP", "GQX") .addSample("NA19600", "./.", "11", "0.7") .setFileId(fileId) .setFormat("GT", "DP", "GQX") .addSample("NA19600", "./.", "11", "0.7") .setFileId(fileId) .setFormat("GT", "DP", "GQX") .addSample("NA19600", "0/1", "11", "0.7")