Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Allpath-LG stuck on CloseUnipathsGaps?

    Hello all,

    I am working on a de novo assembly with Allpaths-LG. I have a 300bp paired-end Illumina MiSeq run with ~400-500 bp fragment size for the overlapping fragment library (~50X), and several 100bp Illumina Hiseq jumping libraries (250bp, 500bp, 2kb, 5kb, 10kb). Allpaths overlaps the Miseq fragment library well, with 99% overlap (I also verified overlaps with PEAR program).

    I am running allpaths with 12 threads on a machine with 2 Intel Xeon X5660, and 192 GB of RAM available - so the machine is decent, and I'm the only user utilizing the machine.

    When Allpaths gets to the CloseUnipathGaps module, it is seeming to get hung up on the "Building those extenders" step. Allpaths ran on this building extenders step for 6 days before I decided to stop it, and try again. The second run is also stuck here for a day now. In other peoples' allpaths logs I have seen, this step only takes a couple minutes to complete.

    Code:
    --------------------------------------------------------------------------------
    Mon Nov 24 18:20:08 2014 run on d017, pid=25701 [Jun  5 2014 15:37:46 R49856 ]
    CloseUnipathGaps NUM_THREADS=12                                                \
                     DIR=/data/FG-183/pulicaria/allpaths-assembly2/run1            \
                     IN_HEAD=frag_reads_corr_cpd                                   \
                     UNIBASES=filled_reads.unibases.k96 UNIBASES_K=96              \
                     USE_LINKS=False OUT_HEAD=filled_reads.gap_closed WORKDIR=tmp  \
                     MM=True _MM_INTERVAL=10 _MM_SUMMARY=False                     \
                     _MM_OUT=/data/FG-183/pulicaria/allpaths-assembly2/run1/makein \
                     fo/filled_reads.gap_closed.CloseUnipathGaps.log.mm.CloseUnipa \
                     thGaps
    --------------------------------------------------------------------------------
    Mon Nov 24 18:20:23 2014: Finding seeds for unipath gaps
    Mon Nov 24 18:20:23 2014: Creating KmerParcel files for FindUnipathGapSeeds
    Mon Nov 24 18:20:23 2014: n_reads = 22919962
    Mon Nov 24 18:21:56 2014: Creating seeds...
    Mon Nov 24 18:21:56 2014: compute seeds[query_ID].size() for each query_ID.
    Mon Nov 24 18:22:03 2014: reserve space for seeds[query_ID].
    Mon Nov 24 18:22:04 2014: store seeds in seeds[query_ID].
    Mon Nov 24 18:22:12 2014: 22919962 number of reads processed.
    Mon Nov 24 18:22:12 2014: 46415402 seeds created.
    Mon Nov 24 18:22:33 2014: Creating KmerParcel files for FindUnipathGapSeeds
    Mon Nov 24 18:22:33 2014: n_reads = 22919962
    Mon Nov 24 18:24:00 2014: Creating seeds...
    Mon Nov 24 18:24:00 2014: compute seeds[query_ID].size() for each query_ID.
    Mon Nov 24 18:24:06 2014: reserve space for seeds[query_ID].
    Mon Nov 24 18:24:07 2014: store seeds in seeds[query_ID].
    Mon Nov 24 18:24:15 2014: 22919962 number of reads processed.
    Mon Nov 24 18:24:15 2014: 46685568 seeds created.
    Mon Nov 24 18:24:37 2014: Loading read quality scores so we can build extenders
    Mon Nov 24 18:24:45 2014: Building those extenders
    
    Sun Nov 30 15:11:45 2014.  Interrupt received (perhaps a ctrl-c).  Stopping.
    I had previously run Allpaths on the Illumina Hiseq jumping library data only, without the MiSeq library, and this problem was not encountered during that assembly. Only with the addition of the MiSeq library did this issue appear.

    Does anyone have experience with Allpaths and have any suggestions? Do you think the 2x300 MiSeq library is too much for Allpaths to handle? (I know the fragment library code is designed for 2x100 reads with 180bp)

    I'd like to be able to try out Allpaths for this assembly. But if it can't complete, I do have the merged overlapping MiSeq reads from PEAR paired-end merger - so I can try out other assemblers, using the merged MiSeq reads as single-end reads. Any recommendations on this front are also appreciated.
    Last edited by CraigJ; 12-01-2014, 11:31 AM.

  • #2
    Try 'CLOSE_UNIPATH_GAPS=False'.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Exploring the Dynamics of the Tumor Microenvironment
      by seqadmin




      The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
      07-08-2024, 03:19 PM
    • seqadmin
      Exploring Human Diversity Through Large-Scale Omics
      by seqadmin


      In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
      06-25-2024, 06:43 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Today, 06:53 AM
    0 responses
    6 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-10-2024, 07:30 AM
    0 responses
    30 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-03-2024, 09:45 AM
    0 responses
    202 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-03-2024, 08:54 AM
    0 responses
    212 views
    0 likes
    Last Post seqadmin  
    Working...
    X