Question: VarDict Amplicon-Aware Calling with barcoded seq data
gravatar for jsneaththompson
9 months ago by
jsneaththompson50 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 9 months ago by Pierre Lindenbaum112k • written 9 months ago by jsneaththompson50
gravatar for Pierre Lindenbaum
9 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum112k 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 9 months ago by Pierre Lindenbaum112k

There's no output when running that command

ADD REPLYlink written 9 months ago by jsneaththompson50
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: 1526 users visited in the last hour