Search in sources :

Example 11 with DiscordantReadPair

use of au.edu.wehi.idsv.DiscordantReadPair in project gridss by PapenfussLab.

the class PathNodeIteratorTest method should_assemble_adjacent_sc_rp.

@Test
public void should_assemble_adjacent_sc_rp() {
    SoftClipEvidence sc = SCE(FWD, withSequence("GTGGTCGACC", Read(0, 50, "4M6S")));
    DiscordantReadPair rp = (DiscordantReadPair) NRRP(SES(10, 200), withSequence("GACCT", DP(0, 25, "5M", true, 1, 1, "5M", false)));
    List<KmerPathNode> result = asCheckedKPN(5, 1000, sc, rp);
    assertCompleteGraph(result, 5);
    assertEquals(2, result.size());
    assertEquals(1, result.stream().mapToInt(n -> n.next().size()).sum());
    assertEquals(1, result.stream().mapToInt(n -> n.prev().size()).sum());
}
Also used : SoftClipEvidence(au.edu.wehi.idsv.SoftClipEvidence) DiscordantReadPair(au.edu.wehi.idsv.DiscordantReadPair) Test(org.junit.Test)

Aggregations

DiscordantReadPair (au.edu.wehi.idsv.DiscordantReadPair)11 Test (org.junit.Test)11 SoftClipEvidence (au.edu.wehi.idsv.SoftClipEvidence)10 ProcessingContext (au.edu.wehi.idsv.ProcessingContext)7 SAMRecord (htsjdk.samtools.SAMRecord)7 SAMEvidenceSource (au.edu.wehi.idsv.SAMEvidenceSource)2 AssemblyEvidenceSource (au.edu.wehi.idsv.AssemblyEvidenceSource)1 DirectedEvidence (au.edu.wehi.idsv.DirectedEvidence)1 DirectedEvidenceOrder (au.edu.wehi.idsv.DirectedEvidenceOrder)1 IdsvVariantContext (au.edu.wehi.idsv.IdsvVariantContext)1 IntermediateFilesTest (au.edu.wehi.idsv.IntermediateFilesTest)1 NonReferenceReadPair (au.edu.wehi.idsv.NonReferenceReadPair)1 VariantCaller (au.edu.wehi.idsv.VariantCaller)1 VariantContextDirectedBreakpoint (au.edu.wehi.idsv.VariantContextDirectedBreakpoint)1 AutoClosingIterator (au.edu.wehi.idsv.util.AutoClosingIterator)1 ImmutableList (com.google.common.collect.ImmutableList)1 Iterables (com.google.common.collect.Iterables)1 Lists (com.google.common.collect.Lists)1 MoreExecutors (com.google.common.util.concurrent.MoreExecutors)1 SortOrder (htsjdk.samtools.SAMFileHeader.SortOrder)1