[maker-devel] failing contigs

Mikael Brandström Durling mikael.durling at slu.se
Fri Feb 22 15:21:16 MST 2013


You are most likely right. When you mention that I remember that I stumbled across that particular problem in another analysis pipeline based on blast I wrote long ago. Interesting that such an arbitrary limit has not been lifted. 
Thanks for the hint. I will try mangling the fasta ids to something shorter. 
Mikael

---
Sent from a crippled computer (a.k.a a phone)

22 feb 2013 kl. 22:40 skrev "Daniel Ence" <dence at genetics.utah.edu>:

> Hi Mikael, 
> 
> Sometimes really long fasta ID's (like the one in your example) cause blast to fail. I'm not super familiar with how short they need to be, but I know I've had this issue in the past. If you can try shortening the IDs, and then retry maker, I think that will solve your problem.
> 
> Thanks,
> Daniel
> 
> Daniel Ence
> Graduate Student
> Eccles Institute of Human Genetics
> University of Utah
> 15 North 2030 East, Room 2100
> Salt Lake City, UT 84112-5330
> ________________________________________
> From: maker-devel-bounces at yandell-lab.org [maker-devel-bounces at yandell-lab.org] on behalf of Mikael Brandström Durling [mikael.durling at slu.se]
> Sent: Friday, February 22, 2013 6:40 AM
> To: maker-devel at yandell-lab.org
> Subject: [maker-devel] failing contigs
> 
> Hi,
> 
> while at trying to find why things don't work for me, I have encountered another problem. Sometimes contigs fail with these lines found in the logs:
> 
> WARNING: Cannot find >0_0_CCPB1339.b1_CCOZ_CCPA_CCPB_CCPC_EXTA, trying to re-index the fasta.
> […]
> stop here:0_0_CCPB1339.b1_CCOZ_CCPA_CCPB_CCPC_EXTA
> doing blastx repeats
> ERROR: Failed while collecting blastn reports
> ERROR: Chunk failed at level:1, tier_type:3
> FAILED CONTIG:scaffold_01
> 
> polishig ESTs
> ERROR: Chunk failed at level:4, tier_type:0
> FAILED CONTIG:scaffold_01
> 
> If I then grep for that specific EST in the fasta file, the sequence is there. Sometimes the issue is resolved by deleting the mpi_blastdb and restarting MAKER, sometimes not.
> 
> Mikael
> _______________________________________________
> maker-devel mailing list
> maker-devel at box290.bluehost.com
> http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org




More information about the maker-devel mailing list