Question: Align HTS on microsatellite
1
gravatar for Nicolas Rosewick
3.1 years ago by
Belgium, Brussels
Nicolas Rosewick7.4k 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 • 658 views
ADD COMMENTlink modified 3.1 years ago • written 3.1 years ago by Nicolas Rosewick7.4k
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: 1253 users visited in the last hour