SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Assembly Comparison (http://seqanswers.com/forums/showthread.php?t=12665)

strob 07-12-2011 05:42 AM

Assembly Comparison
 
Dear all,

Let's say I have performed several assemblies on the same dataset and I would like to compare the results. I know I can look at N50, average contig length, etc. But I was more thinking of comparing the obtained contigs: how much contigs are the same, how much contigs of assembly1 are imbedded in contigs from assembly2, ...

Thanks for the feedback,

Steven

gringer 07-12-2011 07:10 AM

you can use minimus2 (or possibly minimus2BLAT) to compare / combine two different assemblies:
https://sourceforge.net/apps/mediawi...title=Minimus2

for looking at how the assemblies are merged, you can use hawkeye or the listReadPlacedStatus command-line script (these are all part of AMOS). Here's the documentation for the listReadPlacedStatus command:

.USAGE.
listReadPlacedStatus [options] -b <bank path>

.DESCRIPTION.
Prints the status and containing contig(s) of each read as:
iid eid code numcontigs contig iid list

Codes are:
S: Singleton read
P: Placed in a single contig
D: Placed in multiple contigs

.OPTIONS.
-h Display help information
-s Disregard bank locks and write permissions (spy mode)
-v Display the compatible bank version

-S Just list singleton reads
-D Just list duplicate reads
-P Just list placed reads

-E Just list read EIDs
-I Just list read IIDs

.KEYWORDS.
amos bank, reads


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

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.