Question: IndelRealigner and RealignerTargetCreator for amplicon sequencing
0
gravatar for kulvait
2.4 years ago by
kulvait220
Czech Republic
kulvait220 wrote:

Hi,

I am trying to incorporate indel realignment step to my processing pipeline. I use myeloidampliconpanel from Illumina that contains ~1000 amplicons and in total ~100Kb of genome coverage. Prior to the indel realignment step I clean my bams so that no reads are present outside of my amplicons.

Unfortunately when I construct intervals file from bam file and some indel databases (whole genome) I will get the interval file that covers whole genome. I do not understand why it constructs intervals in areas where there is zero coverage?

There is no documentation of interval file. Since I have relatively small genomic area IndelRealigner should do the more work there than in whole genome project. I guess I can somehow pool all the intervals from all my bam files to create list of all possible indels (including those present in all my files) and then run RealignerTargetCreator with this file.

Do anybody of you know what is the correct format of intervals file? I mean if there might be two indels chr1:2-3 and chr1:2-4 should i have interval file with

chr1:2-3
chr1:2-4

or

 chr1:2-4

or even

chr1:2
chr1:3
chr1:4

Thanks, Vojtěch.

ngs amplicon-seq indel • 1.0k views
ADD COMMENTlink modified 2.4 years ago by WouterDeCoster37k • written 2.4 years ago by kulvait220
1
gravatar for WouterDeCoster
2.4 years ago by
Belgium
WouterDeCoster37k wrote:

The -L flag also accepts a bed file with (tab-separated) chromosome - begin - end fields. I would also suggest to use the -ip flag (interval padding) to add extra space up and downstream, to make sure the region isn't too small.

ADD COMMENTlink written 2.4 years ago by WouterDeCoster37k

Thank you, it did the trick.

ADD REPLYlink written 2.4 years ago by kulvait220
0
gravatar for Zaag
2.4 years ago by
Zaag720
Amsterdam
Zaag720 wrote:

You could restrict the construction of the intervals with the -L option to your amplicon intervals.

It should be chr1:2-4

ADD COMMENTlink written 2.4 years ago by Zaag720

Thank you, would it be then wise to just add all amplicon coordinates to interval file?

ADD REPLYlink written 2.4 years ago by kulvait220
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: 1077 users visited in the last hour