You should align only 2 contigs that are most similar to each other. If you look closely at the PDF you posted you can probably make out which contigs are most similar to each other as pairs.
Again this is not going to help you a lot since you have 1.4M contigs.
If you don't know perl/python find a friend who can help parse the blat result file.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Originally posted by GenoMax View PostFor those larger contigs that are noted be similar by blat try using Mauve . You can get additional information from Mauve alignments: http://darlinglab.org/mauve/user-guide/files.html
I did use Mauve. Attached is a pdf with the Mauve results of largest 10 contigs for both assemblies (gerbera_alignment.pdf). However, Mauve does not work when I align the largest 10 contigs from assembly 1 to the whole assembly 2. I guess it is mostly designed for bacterial genome.Attached Files
Leave a comment:
-
For those larger contigs that are noted be similar by blat try using Mauve . You can get additional information from Mauve alignments: http://darlinglab.org/mauve/user-guide/files.html
Leave a comment:
-
Originally posted by GenoMax View Post@Milo: With 1.4M+ input sequences those viewer programs are not going to be useful (and your results are not in blast format either).
You can get an idea of coverage estimate by using Brian's suggestion in this thread: http://seqanswers.com/forums/showthread.php?t=44035 You will have to use the raw data for this. This suggestion is not directly related to question you asked, but may be worth while to do, since you are working with a unknown genome.
I took the longest 10 contigs from each assembly (>8500 bp with a max length of 120 K for assembly 1 and 104 K for assembly 2) and BLAT aligned them together. They seem to be pretty similar but I would like more in depth information about their differences. I am going to give it a try at what you suggested.
Thank you for your help and please let me know if you have any other ideas.
Thank you,
-Emilio
Leave a comment:
-
@Milo: With 1.4M+ input sequences those viewer programs are not going to be useful (and your results are not in blast format either).
You can get an idea of coverage estimate by using Brian's suggestion in this thread: http://seqanswers.com/forums/showthread.php?t=44035 You will have to use the raw data for this. This suggestion is not directly related to question you asked, but may be worth while to do, since you are working with a unknown genome.Last edited by GenoMax; 04-16-2015, 05:19 PM.
Leave a comment:
-
Perhaps you should concentrate on largest (what is the size range of the largest contigs? Are there any that are 10kb and above?) ones first and see if they are related. That would make your search space smaller.
Leave a comment:
-
Originally posted by GenoMax View PostLooks like you are a ways away from having a real assembly. You likely will need some custom scripting to get a meaningful answer since I assume your blat result file (even in PSL format) is probably pretty large.
If the blat results were in blast format and you just wanted to visualize them then http://bioinformatics.oxfordjournals...31/8/1305.full or http://www.biomedcentral.com/1471-2105/15/128 would have been useful.
Do you have a related reference genome available? What is the expected genome size for your samples?
I do no have a related reference genome. Therefore, I aligned both assemblies against each other. My specie is a diploid plant with a 2C DNA value estimated at 5.1 pg (about 5.0 Gb).
I am going to try the ones that you suggested. Is there a better way of comparing both assemblies? How can I get an alignment percentage coverage? Please let me know.
Thank you,
-Milo
Leave a comment:
-
Looks like you are a ways away from having a real assembly. You likely will need some custom scripting to get a meaningful answer since I assume your blat result file (even in PSL format) is probably pretty large.
If the blat results were in blast format and you just wanted to visualize them then http://bioinformatics.oxfordjournals...31/8/1305.full or http://www.biomedcentral.com/1471-2105/15/128 would have been useful.
Do you have a related reference genome available? What is the expected genome size for your samples?
Leave a comment:
-
Originally posted by GenoMax View PostHow big are these genomes? How many contigs are there in each?
Thank you for replying.
Assembly a is 444.5MB and it contains 1.4M contigs
Assembly b is 526.1MB and it contains 1.6M contigs
Please advice.
Thank you,
-Milo
Leave a comment:
-
Comparative Genomics - BLAT
Hello,
I have assembled two related species using Abyss. However, I would like to know if there is a significant difference between the two genomes. A way of doing this is by aligning both assemblies to one another. I used BLAT for the alignment but now I need help visualizing and interpreting the output.psl file. How can I get the total percentage coverage of the alignment? Lets say, assembly A covers x% of assembly B? Is there a better way of doing this? Please let me know. I would really appreciate your help.
Thank you,
-Milo
Latest Articles
Collapse
-
by seqadmin
The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...-
Channel: Articles
05-06-2024, 07:48 AM -
-
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 -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 06:35 AM
|
0 responses
15 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:35 AM
|
||
Started by seqadmin, 05-09-2024, 02:46 PM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
05-09-2024, 02:46 PM
|
||
Started by seqadmin, 05-07-2024, 06:57 AM
|
0 responses
18 views
0 likes
|
Last Post
by seqadmin
05-07-2024, 06:57 AM
|
||
Started by seqadmin, 05-06-2024, 07:17 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
05-06-2024, 07:17 AM
|
Leave a comment: