Question: De novo assembly of a large insertion with a mapped-read seed.
1
gravatar for Pierre Lindenbaum
5.7 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum129k wrote:

short story:

we are trying to identify the content of a large insertion localized at a known place in the human genome. Do you know an existing software that is able to start and extend a de-novo assembly from a read mapped to the REF sequence  ?

 

 

long story:

* my collaborators have sequenced the region of interest using IOn-torrent + mate-pair. As far as I can see, there are a lot of mismatches in the aligned reads.

* my BWA alignement  returns a poor depth in the region of interest (~ 2 )

* CGH shows that the insertion is a junction of a head-to-head duplication

 

 

ADD COMMENTlink modified 5.7 years ago • written 5.7 years ago by Pierre Lindenbaum129k
1
gravatar for Pierre Lindenbaum
5.7 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum129k wrote:

posting ‏@BaCh_mira 's answer on twitter here :

ADD COMMENTlink written 5.7 years ago by Pierre Lindenbaum129k
0
gravatar for Devon Ryan
5.7 years ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

Have you tried making a faux set of reads from the neighboring reference sequence and then seeing what trinity outputs when you try to assemble those along with the other reads? You could just blast the neighboring sequence against the resulting contigs to quickly narrow down which ones are candidates (assuming the results are worthwhile). Something like that might be worth a try at least.

ADD COMMENTlink written 5.7 years ago by Devon Ryan96k

moving to an answer for the discussion. I tried to run vevlet with the reads mapped in the ROI and the unmapped reads: no contig of interest was created (= large number of short contigs)

ADD REPLYlink written 5.7 years ago by Pierre Lindenbaum129k
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: 1430 users visited in the last hour