Question: VarDict Amplicon-Aware Calling with barcoded seq data
gravatar for jsneaththompson
11 weeks 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(
        at com.astrazeneca.vardict.ReferenceResource.retriveSubSeq(
        at com.astrazeneca.vardict.ReferenceResource.getREF(
        at com.astrazeneca.vardict.VarDict.ampVardictNotParallel(
        at com.astrazeneca.vardict.VarDict.start(
        at com.astrazeneca.vardict.Main.main(

Does anyone have any ideas on how to approach this?

ADD COMMENTlink modified 11 weeks ago by Pierre Lindenbaum104k • written 11 weeks ago by jsneaththompson40
gravatar for Pierre Lindenbaum
11 weeks ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum104k 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 11 weeks ago by Pierre Lindenbaum104k

There's no output when running that command

ADD REPLYlink written 11 weeks ago by jsneaththompson40
Please log in to add an answer.


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