Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • building a Mosaik reference for Mouse

    Can anybody point me in the right direction for building a single, full genome, reference for the mouse that i can then align my illumina read data to? I have FASTA reference files, one per chromosome, for the mouse which I downloaded from UCSC. I can pass one of those at a time to MosaikBuild to produce .dat files for each chromosome but that seems a little crazy because that means I'd have to run a single lane of data against each chromosome, 1 at a time.

    If this is how other people do it then that's totally fine - it just seems like I should be able to build a single reference file for the entire genome.
    /* Shawn Driscoll, Gene Expression Laboratory, Pfaff
    Salk Institute for Biological Studies, La Jolla, CA, USA */

  • #2
    Hi there,

    All you have to do is create a concatenated FASTA file and you'll be all set with MOSAIK.

    For example, if I wanted to combine the first four mouse chromosomes into one file, I could type:

    cat mm_ref_chr1.fa >> mouse_ref.fa
    cat mm_ref_chr2.fa >> mouse_ref.fa
    cat mm_ref_chr3.fa >> mouse_ref.fa
    cat mm_ref_chr4.fa >> mouse_ref.fa

    You could keep doing this for all of the mouse chromosomes or if you're savvy at creating bash scripts, you could pretty much automate the above in a small script.

    Cheers,

    // Michael

    Comment


    • #3
      cool thanks. it just wasn't clear in the documentation that you could just cat files together to make one larger reference. now i just need to figure out this jump database thing and i'll be off and running. Mosaik chews up some serious RAM and i've only got 16 GB on the system I'm running it on. looks like a jump database will help for running full genome alignments on this system.
      /* Shawn Driscoll, Gene Expression Laboratory, Pfaff
      Salk Institute for Biological Studies, La Jolla, CA, USA */

      Comment


      • #4
        so i made this cat'd reference file (2.6 GB) and compiled it down. then i made a jump database and started a run with MosaikAligner using the jump database of this full genome reference. i ran pretty much the default settings listed in the manual except with only 4 cpu cores. it looks like it munched up about 19GB of RAM to load the jump database files into memory but once the alignment actually started it wasn't using all 4 cores - it was only using about 4% of the CPU and it was processing only 3.5 reads per second with an ETA of 53 DAYS. what could be going wrong?
        /* Shawn Driscoll, Gene Expression Laboratory, Pfaff
        Salk Institute for Biological Studies, La Jolla, CA, USA */

        Comment


        • #5
          cat all files

          maybe concatting different files might get little tedious. i tried:
          cat chr*.fa >> human_ref.fa

          it worked well.

          Comment


          • #6
            Are you still using your system with 16GB of RAM?
            You are probably swapping if it's using 19GB...

            I found reading the manual to get the right parameters made a huge difference to the reads per second. The magic parameters I found were
            -bw 13 -act 20 -mm 4 -mhp 100
            That took me from a few reads per second to 700 - 800 per second.

            I also needed around 20GB of RAM for the jump database.

            Comment


            • #7
              thanks for sharing the magic.
              /* Shawn Driscoll, Gene Expression Laboratory, Pfaff
              Salk Institute for Biological Studies, La Jolla, CA, USA */

              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, Yesterday, 06:37 PM
              0 responses
              10 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              10 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              51 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-21-2024, 07:32 AM
              0 responses
              67 views
              0 likes
              Last Post seqadmin  
              Working...
              X