Question: VarDict Amplicon-Aware Calling with barcoded seq data
0
gravatar for jsneaththompson
6 days ago by
jsneaththompson40 wrote:

I'm trying out different variant calling software on non-matched somatic samples, where the pcr has been performed using an Agilent kit called Haloplex HS.

At the moment I'm trying VarDict, and I want to try to run VarDict in amplicon-aware mode. The problem is, amplicon-aware mode requires a bed file with start and end positions both with primers and without primers, and Haloplex HS doesn't use standard pcr primers. As can be seen in this image, Haloplex has a molecular barcodes and indexes in addition to a generic pcr primer.

Is it possible to run amplicon-aware calling on barcoded-data? I've tried adding two columns to my amplicons bed file which are the original amplicon start/end points moved one base back/forward, but I get the following error:

Exception in thread "main" htsjdk.samtools.SAMException: Malformed query; start point 139389396 lies after end point 138394717
        at htsjdk.samtools.reference.IndexedFastaSequenceFile.getSubsequenceAt(IndexedFastaSequenceFile.java:215)
        at com.astrazeneca.vardict.ReferenceResource.retriveSubSeq(ReferenceResource.java:31)
        at com.astrazeneca.vardict.ReferenceResource.getREF(ReferenceResource.java:50)
        at com.astrazeneca.vardict.VarDict.ampVardictNotParallel(VarDict.java:5577)
        at com.astrazeneca.vardict.VarDict.start(VarDict.java:69)
        at com.astrazeneca.vardict.Main.run(Main.java:138)
        at com.astrazeneca.vardict.Main.main(Main.java:25)

Does anyone have any ideas on how to approach this?

ADD COMMENTlink modified 6 days ago by Pierre Lindenbaum102k • written 6 days ago by jsneaththompson40
0
gravatar for Pierre Lindenbaum
6 days ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum102k wrote:

what's the output of

awk -F '\t' '(int($3)<int($2))' input.bed

if there is any ouput, then your bed is malformed (chromEnd>chromStart)

ADD COMMENTlink written 6 days ago by Pierre Lindenbaum102k

There's no output when running that command

ADD REPLYlink written 5 days ago by jsneaththompson40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1431 users visited in the last hour