Examples with DiscordantReadPair au.edu.wehi.idsv.DiscordantReadPair used on opensource projects

Example 6 with DiscordantReadPair

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

the class NonReferenceContigAssemblerTest method should_truncate_anchored_polyA_expansion_bwd.

@Test
public void should_truncate_anchored_polyA_expansion_bwd() {
    ProcessingContext pc = getContext();
    pc.getAssemblyParameters().k = 4;
    SoftClipEvidence sce = SCE(BWD, withSequence("AAAAGTAC", Read(0, 100, "4S4M")));
    DiscordantReadPair dp = (DiscordantReadPair) NRRP(SES(11, 100), withSequence("CAAAAG", DP(0, 110, "6M", false, 1, 1, "6M", true)));
    List<SAMRecord> output = go(pc, true, sce, dp);
    assertEquals(1, output.size());
    assertEquals("AAAAGTAC", S(output.get(0).getReadBases()));
}

Example 7 with DiscordantReadPair

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

the class NonReferenceContigAssemblerTest method should_call_overlapping_sc_rp.

@Test
public void should_call_overlapping_sc_rp() {
    ProcessingContext pc = getContext();
    pc.getAssemblyParameters().k = 4;
    SoftClipEvidence sce = SCE(FWD, withSequence("ACGTGGTCGACC", Read(0, 50, "6M6S")));
    DiscordantReadPair e = (DiscordantReadPair) NRRP(SES(10, 200), withSequence("GACCTCTACT", DP(0, 25, "10M", true, 1, 1, "10M", false)));
    List<SAMRecord> output = go(pc, true, sce, e);
    assertEquals(1, output.size());
    assertEquals("ACGTGGTCGACCTCTACT", S(output.get(0).getReadBases()));
}

Example 8 with DiscordantReadPair

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

the class NonReferenceContigAssemblerTest method should_truncate_anchored_polyA_expansion_fwd.

@Test
public void should_truncate_anchored_polyA_expansion_fwd() {
    ProcessingContext pc = getContext();
    pc.getAssemblyParameters().k = 4;
    SoftClipEvidence sce = SCE(FWD, withSequence("GTACAAAA", Read(0, 10, "4M4S")));
    DiscordantReadPair dp = (DiscordantReadPair) NRRP(SES(11, 100), withSequence("CAAAAT", DP(0, 1, "6M", true, 1, 1, "6M", false)));
    List<SAMRecord> output = go(pc, true, sce, dp);
    assertEquals(1, output.size());
    assertEquals("GTACAAAA", S(output.get(0).getReadBases()));
}

Example 9 with DiscordantReadPair

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

the class NonReferenceContigAssemblerTest method should_call_simple_bwd_RP.

@Test
public void should_call_simple_bwd_RP() {
    ProcessingContext pc = getContext();
    pc.getAssemblyParameters().k = 4;
    MockSAMEvidenceSource ses = SES(100, 200);
    DiscordantReadPair e = (DiscordantReadPair) NRRP(ses, withSequence("ACGTGGTCGACC", DP(0, 450, "12M", false, 1, 1, "12M", true)));
    List<SAMRecord> output = go(pc, true, e);
    assertEquals(1, output.size());
    SoftClipEvidence ass = SoftClipEvidence.create(AES(ses), BWD, output.get(0));
    assertEquals(e.getBreakendSummary().localBreakend(), ass.getBreakendSummary());
    assertEquals("ACGTGGTCGACC", S(ass.getBreakendSequence()));
    assertEquals("", S(ass.getAnchorSequence()));
}

Example 10 with DiscordantReadPair

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

the class NonReferenceContigAssemblerTest method should_call_multiple_nonoverlapping_contigs.

@Test
public void should_call_multiple_nonoverlapping_contigs() {
    ProcessingContext pc = getContext();
    pc.getAssemblyParameters().k = 4;
    SoftClipEvidence sce = SCE(FWD, withSequence("ACGTGGTCGACC", Read(0, 5, "6M6S")));
    DiscordantReadPair e = (DiscordantReadPair) NRRP(SES(100, 200), withSequence("GACCTCTACT", DP(0, 25, "10M", true, 1, 1, "10M", false)));
    List<SAMRecord> output = go(pc, true, sce, e);
    assertEquals(2, output.size());
    assertEquals("ACGTGGTCGACC", S(output.get(0).getReadBases()));
    assertEquals("NN" + "GACCTCTACT", S(output.get(1).getReadBases()));
}