use of au.edu.wehi.idsv.AssemblyAttributes in project gridss by PapenfussLab.
the class NonReferenceContigAssemblerTest method should_set_assembly_attributes.
@Test
public void should_set_assembly_attributes() {
ProcessingContext pc = getContext();
pc.getAssemblyParameters().k = 4;
SoftClipEvidence sce = SCE(FWD, withSequence("ACGTGGTCGACC", Read(0, 5, "6M6S")));
List<SAMRecord> output = go(pc, true, sce);
AssemblyAttributes attr = new AssemblyAttributes(output.get(0));
assertNotNull(attr);
assertEquals(1, attr.getAssemblySupportCount());
}
use of au.edu.wehi.idsv.AssemblyAttributes in project gridss by PapenfussLab.
the class PositionalAssemblerTest method should_set_assembly_direction.
@Test
public void should_set_assembly_direction() {
ProcessingContext pc = getContext();
pc.getAssemblyParameters().k = 4;
pc.getAssemblyParameters().maxExpectedBreakendLengthMultiple = 10;
List<DirectedEvidence> e = new ArrayList<>();
// ACGTACGTACGTACGTACGT
// 12345678901234567890
// MMMM----MMMM
// MMMM>>>>
e.add(SCE(FWD, withSequence("ACGTCCGGACGG", Read(1, 1, "4M8S"))[0]));
e.add(SCE(FWD, withSequence("ACGGTTGT", Read(1, 9, "4M4S"))[0]));
e.sort(DirectedEvidenceOrder.ByStartEnd);
ArrayList<SAMRecord> output = Lists.newArrayList(new PositionalAssembler(pc, AES(pc), new SequentialIdGenerator("asm"), e.iterator(), BreakendDirection.Forward));
for (SAMRecord r : output) {
assertEquals(BreakendDirection.Forward, new AssemblyAttributes(r).getAssemblyDirection());
}
}
Aggregations