Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #31
    Originally posted by GenoMax View Post
    This is not working for me:

    Code:
    tar: This does not look like a tar archive
    tar: Skipping to next header
    tar: Read 4461 bytes from ./sorted-reference.xml
    tar: Error exit delayed from previous errors
    make: *** [Temp/sorted-reference.xml] Error 2
    My mistake. Apologies. It is not the sorted-reference.xml but the tarball created by 'isaac-pack-reference':

    Code:
    rm -rf Temp
    isaac-unpack-reference -j 1 -w 6 -i packed-reference.tar.gz

    Comment


    • #32
      Commands used for the final steps in a nutshell.

      Code:
      $ isaac-pack-reference -j 1 -r ./sorted-reference.xml -o ./packed-reference.tar.gz
      
      $ isaac-unpack-reference -j 1 -w 6 -i ./packed-reference.tar.gz
      The end result was a set of 64 files

      kmer-positions-32-00.dat through kmer-positions-32-63.dat
      And one

      Code:
      2uniqueness.16bpb.gz
      file.

      I have started a new isaac2 genome creation job for the MM9 genome with -w 6 option upfront.
      Last edited by GenoMax; 07-29-2015, 03:11 AM.

      Comment


      • #33
        Got the same just 5 minutes ago :-)

        So the default for isaac-sort-reference should be changed or, alternatively, it should always be called with '--mask-width 6'.

        Comment


        • #34
          I had started an isaac2 index creation job for mm9 genome (with -w 6). It has been running for a week and still making files in Temp directory.

          Comment


          • #35
            In an attempt to make it easier to use Isaac2, we will make the packed index reference for commonly used genomes on BaseSpace. At the moment, the only 2 genomes available are hg19 and mm9. Feel free to request other genomes.

            Also, the issues and recommendations around indexing genomes are summarized on the isaac2 github wiki page "Reference Indexes".

            The link to the already indexed genomes in basespace might change in the future, please refer to the wiki page on github for updates.

            Hopefully, this will help.

            Come

            Comment


            • #36
              Hallo again ;-)

              we are now with Isaac3. Cool .. ;-)

              Creating indices for grch38 and grcm38 leaves some open questions:

              I have run index creation as follows (mask-width 0 is the default, I just put it there as a "reminder" for future index creation runs):

              Code:
              isaac-sort-reference \
                --output-directory iSAACindex \
                --jobs 1 \
                --mask-width 0 \
                --genome-file genome.fa
              That left me with exact 3 files and a 1.1TiB Temp folder:

              Code:
              -rw-rw-r-- 1 klages klages 618M 2016.08.26 01:05:08 2repeatness.8bpb.gz
              -rw-rw-r-- 1 klages klages 678M 2016.08.25 22:19:13 2uniqueness.8bpb.gz
              -rw-rw-r-- 1 klages klages 108K 2016.08.26 01:05:09 sorted-reference.xml
              drwxrwxr-x 2 klages klages 8.0K 2016.08.26 01:05:09 Temp
              make reported
              Code:
              [all]    INFO: All done!
              At least it is "packable" by isaac-pack-reference.

              hg19-packed-reference.tar.gz from BaseSpace (btw, would be fine to have some grch38/grcm38 though) shows:

              Code:
              -rwxr-x--- rpetrovski/aladdin 644685308 2014-11-19 21:38 2uniqueness.16bpb.gz
              -rw-r--r-- rpetrovski/aladdin 386961748 2014-11-20 13:03 neighbors-1or2-16.1bpb
              -rw-r--r-- rpetrovski/aladdin 386961748 2014-11-20 13:06 neighbors-1or2-32.1bpb
              -rwxr-xr-- rpetrovski/aladdin 3157608038 2014-11-20 12:53 genome.fa
              -rw-r--r-- rpetrovski/aladdin      48044 2014-11-20 12:54 sorted-reference.xml
              * Is that a complete and valid index??
              * Do I still need Temp for any task after index creation?
              * What are the differences compared to isaac2 indices?

              best,
              Sven

              Comment


              • #37
                Dear all,
                Can iSAAC work on mac os platform?

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Strategies for Sequencing Challenging Samples
                  by seqadmin


                  Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                  03-22-2024, 06:39 AM
                • seqadmin
                  Techniques and Challenges in Conservation Genomics
                  by seqadmin



                  The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                  Avian Conservation
                  Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                  03-08-2024, 10:41 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 03-27-2024, 06:37 PM
                0 responses
                13 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-27-2024, 06:07 PM
                0 responses
                11 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-22-2024, 10:03 AM
                0 responses
                53 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-21-2024, 07:32 AM
                0 responses
                69 views
                0 likes
                Last Post seqadmin  
                Working...
                X