Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Demultiplexing combinatorial barcodes from paired end data

    Hey, I've been trying to figure out how to demultiplex paired-end 250 miseq data that has three different barcodes.

    The samples all come from 96 well plates. I have two 5' barcodes representing plate number as well as plate column. These two barcodes are linked via a common adapter sequence. I then have one 3' barcode that represents the plate column.

    I am fairly sure I would be able to demultiplexing the data in 2 rounds in which I first demultiplex based upon each plate/column combination, and then demultiplex the result based upon the 3' barcode.

    This seems fairly cumbersome as I have 2880 unique samples. Is there some sort of software in which I can provide simply the three unique barcodes, and then demultiplex in one step?

    Also, I usually use the A5Pipeline (https://code.google.com/p/ngopt/wiki/A5PipelineREADME) to error correct sequence reads, remove Illumina adapter sequences, and assemble the paired-reads. Would I be able to assemble the reads using A5 and then demultiplex the data, or do I have to run the pipeline on each individual sample after first demultiplexing. Thanks!

  • #2
    Does anyone have experience with combinatorial barcodes?

    Comment


    • #3
      How long are these barcodes? It may be possible to write some custom code to do this in one step but without seeing what you reads look like (or a good schematic) no one is going to be able to help. It may be best to find someone locally who can work with you interactively.

      If the A5Pipeline changes the context of the barcodes it will make things complicated. You may have to run the pipeline post demultiplexing.

      Comment


      • #4
        This will do it https://code.google.com/p/jmhc/

        Comment


        • #5
          Originally posted by JackieBadger View Post
          Can you import fastq files? Description states that it operates on fasta files (internally?). I don't think this will do 3 barcodes in a single pass.

          Comment


          • #6
            use skewer

            Hi ddelalca,
            You may use skewer for this purpose. An example command line is as follows:
            Code:
            skewer -x V3-forward.txt -y V4-reverse.txt --mode ap -t 4 ampliconReads-pair1.fastq ampliconReads-pair2.fastq assigned/
            If you only use some of the combinations, you may specify the matrix by -M as well.
            Good luck!
            Attached Files

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