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

Join Date: May 2011
Posts: 838

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