Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Bowtie2 with tight matching constraints on repeat sequences theoriste Bioinformatics 0 03-21-2013 11:10 AM
generating heatmaps of repeat density in genome sequences bobvanburen Bioinformatics 1 05-10-2012 08:15 PM
de novo assembly of repeat elements HESmith Bioinformatics 12 11-09-2011 04:11 PM
Determining Repeat Sequences clostridium40 Bioinformatics 0 09-02-2011 02:23 PM
repeat sequences/large files in galaxy Giles Bioinformatics 2 06-27-2011 12:08 PM

Thread Tools
Old 09-28-2014, 11:24 PM   #1
Location: US

Join Date: May 2012
Posts: 80
Default Repeat sequences-Genome assembly

I have WGS sequence of illumina paired-end reads (2x101 bp) of the plant crop. I have de novo assembled the reads using velvet.

I would like to do two things

1. To do de novo genome assembly and identify repetitive sequences. I know i can do with repeat masker, is there any software other than this?
2.Remove prior repetitive sequence from illumina reads and do de novo assembly.

Please let me know software to find repetitive sequence before assembly and after assembly.
bioman1 is offline   Reply With Quote
Old 09-29-2014, 12:15 AM   #2
Brian Bushnell
Super Moderator
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707

As I mentioned in this thread, you can use a tool such as BBNorm to get a kmer frequency spectrum of your reads, and look at the peak locations to determine the coverage range for 1-copy, 2-copy, etc peaks. Then you can bin them by kmer frequency. Illumina coverage is highly uneven so this does not work perfectly. And of course it's harder with ploidy over 1, but still possible... to some extent. Then you can just assemble the reads in the unique bin.

The repeat stuff may not assemble very well with Velvet. But another way to do it is by mapping the input reads to your assembly, and annotate/sort the contigs by coverage; contigs with 2x median coverage may be 2-copy repeats, and so forth. Velvet already annotates contig names with their estimated coverage, so you could use that instead. Either way, once you isolate the repeat contigs and unique contigs, separate them into two fasta files and map the reads to them with BBSplit to get three outputs - a pile of single-copy reads, a pile of repeat reads, and a pile of unmapped reads.
Brian Bushnell is offline   Reply With Quote
Old 09-29-2014, 10:05 AM   #3
Location: NYC

Join Date: Aug 2010
Posts: 48

the introduction of this paper contains nice brief descriptions of many of the softwares used to find and analyse genomic repeats.
ssully is offline   Reply With Quote

genome assembly, repeat sequence

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 08:44 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO