Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
I understood that with the old Newbler repeat reads where simply not used in the assembly. While with the upgrade they are now used once. I have no idea how these end up in the final contigs. We only did such an experiment once, and purely looking at contig number and contig length, binning seemed to work very well.
-
Originally posted by Tuxido View PostThat would depend on your genome size and your average read length (i.e. it depends on your coverage).
So starting from I guess 50x coverage it might be worthwhile to give it a try. I would just try a few different bin settings and see which one gives you the best contigs. (I only did this once myself and used 18 bins on 500Mb of sequence data for a 4.5Mb genome, which improved assembly considerably).
Leave a comment:
-
Originally posted by DNAjunk View Post>flxlex:
I have compared the contigs to a reference genome. In some cases the phrap and in others the 454 contig is in agreement with the reference.
Of course, I have used quality scores also for phrap (extracted from the sff file). What information does the "flowgrams" contain?
sffinfo yourfile.sff
you will see among the ouput something like this:
Flowgram: 1.01 0.04 1.01 0.06 0.06 0.98 0.05 1.04
2.33 1.16 1.10 0.06 0.21 0.89 0.07 2.00 0.12 0.08
1.89 0.10 0.45 1.02 1.84 0.17 0.92 0.34 0.09 0.99
With flow order TACG, this means that at the first flow, T, the signal intensity was just over 1, meaning most likely 1 T. The next flow, A, gave no signal. etc
newbler uses this information for the assembly/quality scoring
Leave a comment:
-
estimate genome size
>Tuxido:
Thanks for the link and advice!
It's a de novo assemly, and I don't have any reference or knowledge of the genome size. Is there a way to estimate the genome size after a 454 run with 800k reads and an average read length of 354bp?
Leave a comment:
-
That would depend on your genome size and your average read length (i.e. it depends on your coverage).
So starting from I guess 50x coverage it might be worthwhile to give it a try. I would just try a few different bin settings and see which one gives you the best contigs. (I only did this once myself and used 18 bins on 500Mb of sequence data for a 4.5Mb genome, which improved assembly considerably).
There's also a paper about this from Bas Dutilh
http://bioinformatics.oxfordjournals...bstract/btp377 but he uses it when working with metagenomes.
Leave a comment:
-
Thank you for your kind replies!
I appreciate any advice and suggestion because I am very new to bioinformatics and just have run my first assembly.
>flxlex:
I have compared the contigs to a reference genome. In some cases the phrap and in others the 454 contig is in agreement with the reference.
Of course, I have used quality scores also for phrap (extracted from the sff file). What information does the "flowgrams" contain?
>Tuxido:
In another study, I am doing a de novo assembly. The sequencer made 800'000 reads. I have executed runAssembly on the whole set of reads. The output was contigs ranging upto 16'000bp. Would you still bin the reads, and if so, how many into each bin?
Leave a comment:
-
We have good experience with the performance of 454 Newbler assembler. If you're doing denovo assembly however it might be a good idea to use an iterative approach. Dividing your reads into bins, assembling these separately and then assembling the contigs that you got. Especially useful if you have a bit too much coverage.
The notrim option turns off the additional read trimming in the assembly (See the 454TrimStatus.txt that is generated with the assembly).
Leave a comment:
-
I would be really careful in using phrap for assembly of 454 reads. The program was designed for assembly of Sanger reads, which have very different kind of errors than 454 reads. Even though the contigs may be longer in the phrap assembly, I would do some checks to see if they are correct (do you have a reference genome to compare the contigs to? PCR?).
runAssembly (newbler) is so far (on of) the best assembly program for 454 reads (some people might disagree?). Be sure to feed it sff files, as they contain more information (flowgrams) than just fasta and quality files.
Leave a comment:
-
454 runAssembly contig lengths
Hi
I've made assemblies of several 454 reads (in sff and fasta format) by using two different programs: runAssembly (reads in sff format) and phrap (reads in fasta format).
The contig length of the runAssembly program is much smaller than the one of phrap.
Does anybody how to set the parameters/options of the runAssembly program in order to modify the contig length?
I have tried the "notrim" option but it didn't have any effect on the contig length.
Thanks for your advice!
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
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 -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-25-2024, 11:49 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
04-25-2024, 11:49 AM
|
||
Started by seqadmin, 04-24-2024, 08:47 AM
|
0 responses
17 views
0 likes
|
Last Post
by seqadmin
04-24-2024, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
62 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
Leave a comment: