SEQanswers (
-   Bioinformatics (
-   -   Assembly Comparison (

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,


gringer 07-12-2011 07:10 AM

you can use minimus2 (or possibly minimus2BLAT) to compare / combine two different assemblies:

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:

listReadPlacedStatus [options] -b <bank path>

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

-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

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.