Thanks for the clarification, that sounds good! So, considering I'm trying to merge two sets of contigs returned by Newbler and AbySS, this approach can make feasible to compare directly the coverage of the datasets, as Newbler comes out with a bunch of files (very useful about coverage and so on), while AbySS output needs more processing.
Does the "less stringent" assembly phase take into account the repeat prediction and try to infer this information in the assembly itself?
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Hi RLB_84, that's a very good question.
In addition to the contig fasta files, Zorro takes as input the reads file (a subsample of WGS reads). The reads are used only to allow us to identify repeats in the contigs. I will explain technically
1-Zorro counts the occurences of 22-mers in the reads file supplied by the user
2-22-mer words that are unique in the genome should occur proportionally to the genome coverage
3-22-mer words that represent repeats should occur at least twice the peak coverage
4-we select the 22-mer words that occur at least twice the mode of the distribution. These 22-mer words are used to mask the contig files using bowtie.
This technique is used by many ab initio repeat detection software. We do not need to screen repeat libraries and, even if newbler (or other software) has collapsed the repeats, we coulod still detect them.
Leave a comment:
-
hi, I read about your pipeline and it seems quite interesting, as I'm a newbie user of MUMmer. So I hope I could try your assembler as soon as possible. By the way, just one question about "masking": when I consider a contig dataset returned by Newbler, I should keep in mind that many of them can be repeated and I can't figure it from the contig fasta file (I need the other Newbler outputs). So, how does your assembler deal with the probable occurrence of a contig in more than one genome locus?
Leave a comment:
-
Zorro: The Masked Assembler (first public release)
ZORRO is an hybrid sequencing technology assembler. It takes 2 sets of pre-assembled contigs and merge them into a more contiguous and consistent assembly. We have already tested Zorro with Illumina Solexa and 454 from some of organisms varying from 3Mb to 100Mb. The main caracteristic of Zorro is the treatment before and after assembly to avoid errors.
The ZORRO project is maintained by Gustavo Lacerda, Ramon Vidal and Marcelo Carazzole and were first used in this Yeast assembly: Genome structure of a Saccharomyces cerevisiae strain widely used in bioethanol production
ZORRO needs to be better documented and has not undergone enough testing. If you want to discuss the pipeline you can join the mailing list: zorro-google group
ZORRO PIPELINE
Zorro is based on the minimus2 pipeline (AMOS package) and uses MuMMer,
AMOS and bowtie in its internals. Zorro takes 2 contigs fasta files as
input (representing assembled contigs from a whole genome assembly)
and one fasta file containing some of the reads used for assembly
(only 10X coverage is enough, more will slow down the pipeline and
consume more resources).
Zorro initial phase detect inconsistencies in the assemblies and split
the contigs where they occur. Next, zorro counts k-mers (default k=22)
in the reads and use the k-mer count table to detect and mask repeats
in both assembly1 and assembly2. After repeat masking, zorro uses nucmer
to detect overlaps between assembly1 and assembly2 (no overlaps between
contigs from the same assembly are allowed). All overlaps found in this
phase are expected to be between unique regions (because repeats are
masked). The overlaps are used to layout and generate consensus for the
merged contigs, using AMOS tools. The merged contigs are built using the
unmasked contigs, so the final merged assembly should include the repeat
regions.
Another round of assembly, less stringent, tries to merge contigs that
were not included in the first Zorro phase. All the contigs are outputted
to <prefix>.ZORRO.fasta. We recommend the use of SSPACE to scaffold the
ZORRO contigs.
Zorro Website: www.lge.ibi.unicamp.br/zorro
Latest Articles
Collapse
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
-
by seqadmin
Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...-
Channel: Articles
03-22-2024, 06:39 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
30 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
32 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
28 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
53 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Leave a comment: