View Single Post
Old 07-12-2011, 07:10 AM   #2
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 838
Default

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
gringer is offline   Reply With Quote