My amos bank appears to have errors in it. I have detected these problems using both analyzeSNPs and hawkeye. My assembly is hybrid and was made by the celera assembler.
I used: toAmos -f gsb/antrc230_0.1_hybrid.frg -a gsb/gs35/9-terminator/gs35.asm -o - -utg | bank-transact -m - -b gsb-utg.bnk -c
I have included screen shots of a unitig with no reads showing in hawkeye (I realise that contigs can have missing reads because of surrogates), the basic statistics of that unitig, a unitig with consensus that does not match the reads, the corresponding section of the contig with the same reads with correct consensus and a unitig with some reads that seem misplaced.
The assembly used has not yet been published and thus I have not included it but I may be able to provide it confidentially.
Log in to post a comment.