Question: Align HTS on microsatellite
1
gravatar for Nicolas Rosewick
4.9 years ago by
Belgium, Brussels
Nicolas Rosewick9.3k wrote:

Hi,

I've a bunch of HTS data enriched for specific microsatellite (MS) regions. My idea was to create a bwa index based on these microsatellite sequences as :

Here for one of the MS region, produced sequences as:

seq1: [ Upstream flanking region ][ Repeat ][ Downstream flanking region ]
seq2: [ Upstream flanking region ][ Repeat ][ Repeat ][ Downstream flanking region ]
seq3: [ Upstream flanking region ][ Repeat ][ Repeat ][ Repeat ][ Downstream flanking region ]
seq4: [ Upstream flanking region ][ Repeat ][ Repeat ][ Repeat ][ Repeat ][ Downstream flanking region ]
.
seqN: [ Upstream flanking region ][ Repeat ]...[ Repeat ][ Downstream flanking region ]

And then align the reads (paired-end 2x100bp) on it, and count how many reads for each different variant of MS. Will bwa manage the alignment; and not report multiple alignment for one read ?

> edit : Repeats can be mono / di / tri nucleotides

Maybe someone has an other idea / aligner / method to compute such thing ?

Thanks

 

microsatellite dna-seq • 922 views
ADD COMMENTlink modified 4.9 years ago • written 4.9 years ago by Nicolas Rosewick9.3k
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: 2285 users visited in the last hour
_