Tweet
Extract scaffold info from Mira PE 454 assembly?
Does Mira create scaffolds during paired-end read sequence assemblies? gsAssembler does, so I just assumed Mira would. But maybe I need to run downstream software (Bambus?) to extract this information?
--
Phillip
--
Phillip
-
-
No, MIRA doesn't create scaffolds; you should use an external tool like BAMBUS.
cheers,
Sven -
You can use the MIRA output for scaffolding with Bambus (see tutorial)
You can also use the MIRA output with SSPACE with only minimal manipulation:
JComment
-
Bambus Obsolete?
I am getting the mysterious "script failed" error mentioned in this thread:
Therein it is posited that there is some maximum number of contigs that can be processed. Anyone know what that maximum is or have other insights?
Also, given that no modifications to Bambus have been checked in to its SourceForge site since 2005 -- isn't this really a pre-next gen program?
--
PhillipComment
-
Bambus works for gsAssembler .ace file
On the other hand, for smaller 454 jobs it does work for Mira assemblies. And with a few modifications of the normal protocol, works directly on a gsAssembler .ace file.
Not sure if gsAssembler creates a .ace file by default. If not, use the appropriate option to have one created. Then create the .contig file from the .ace file as normal using the amostools utility:
ace2contig -i 454Contigs.ace 454Contigs.contig
The ace file created by gsAssembler uses a different paradigm to denote forward and reverse reads. So I used the following 454Contigs.mates file:
As always for bambus, the field delimiter (shown above as whitespace) for a .mates must be tabs. So a copy/paste from this post will not work without modification of the resultant file.Code:library all 1000 6000 (.......).* pair (.*)_left (.*)_right
should do the trick. Note that is two spaces prior to the "*", not one. If you only put one, that would allow perl to place a tab in between every character.Code:perl -i -pe 's/ */\t/g' 454Contigs.mates
One other catch for me. I had been skipping the .conf file, allowing goBambus to substitute its own default one. But I got the following this time:
So I created a .conf file, 454Contigs_bambus.conf with the following in it:Step 400: Running scaffolder
Error: Priority not specified: at least one library must be assigned a priority in module grommit
Grommit(/usr/local/bambus-2.33/bin/grommit -i 454Contigs_bambus.inp -o 454Contigs_bambus.out.xml -C 454Contigs_bambus.grommit.conf --append --logfile goBambus.log --debug 1) script failed
Code:# Priorities priority ALL 1 # Redundancies redundancy 2 # min group size mingroupsize 0
Then the following command:
appears to work.Code:goBambus -c 454Contigs.contig -m 454Contigs.mates -o 454Contigs_bambus -C 454Contigs_bambus.conf
Of course, gsAssembler should create scaffolds itself. But, for this project it did a very poor job, creating a few tiny scaffolds. So it is nice to have Bambus as a backup.
--
PhillipComment
-
The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...11-06-2024, 07:24 PM -
Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...10-18-2024, 07:11 AM
Comment