Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Demultiplexing problem

    Hi All,

    I have been trying to demultiplex some Illumina (101bp PE, dual barcoded) libraries using Picard - e.g.:

    Code:
    java -Xmx2g -jar /seq/software/picard/1.626/bin/ExtractIlluminaBarcodes.jar BASECALLS_DIR=/analysis/140615/_temp/140615.H9FT9ADXX/Data/Intensities/BaseCalls/ LANE=2 READ_STRUCTURE=101T8B101T BARCODE_FILE=/analysis/140615/barcodeData.2 METRICS_FILE=/analysis/140615/_logs/barcode.metrics.H9FT9ADXX.2.txt MAX_MISMATCHES=0 MINIMUM_BASE_QUALITY=25 NUM_PROCESSORS=1
    I have done this a million times before, but with these new libraries I keep getting the following error messages:

    Code:
    Barcode AGGCAGAA specified more than once in /analysis/140615/barcodeData.2
    Barcode CAGAGAGG specified more than once in /analysis/140615/barcodeData.2
    Barcode CAGAGAGG specified more than once in /analysis/140615/barcodeData.2
    .
    .
    .
    I have never seen this happen before and I have tried everything (I can think of...) to fix the problem. Basically it appears as if Picard isn’t accepting the dual barcodes and just treats them individually - which is strange, since I have demultiplexed many runs that were dual barcoded. If I just to try to demultiplex samples where both barcodes are different, then everything goes through.

    My barcode data files look like this:
    Code:
    barcode_name	library_name	barcode_sequence_1	barcode_sequence_2
    H9FT9ADXX-EM-096_NuGEN	EM-096_NuGEN	AGGCAGAA	TATCCTCT
    H9FT9ADXX-EM-098_NuGEN	EM-098_NuGEN	AGGCAGAA	AGAGTAGA
    H9FT9ADXX-EM-095B_NuGEN	EM-095B_NuGEN	CAGAGAGG	TAGATCGC
    H9FT9ADXX-G3691-1_NuGEN	G3691-1_NuGEN	CAGAGAGG	CTCTCTAT
    H9FT9ADXX-K562B_NuGEN	K562B_NuGEN	CAGAGAGG	CTAAGCCT
    H9FT9ADXX-EM-095_NuGEN	EM-095_NuGEN	CTCTCTAC	CTAAGCCT
    H9FT9ADXX-G3670-1_NuGEN	G3670-1_NuGEN	CTCTCTAC	TAGATCGC
    H9FT9ADXX-G3679-1_NuGEN	G3679-1_NuGEN	CTCTCTAC	TATCCTCT
    H9FT9ADXX-G3680-1_NuGEN	G3680-1_NuGEN	CTCTCTAC	AGAGTAGA
    H9FT9ADXX-G3681-1_NuGEN	G3681-1_NuGEN	CTCTCTAC	GTAAGGAG
    H9FT9ADXX-G3682-1_NuGEN	G3682-1_NuGEN	CTCTCTAC	ACTGCATA
    H9FT9ADXX-G3683-1_NuGEN	G3683-1_NuGEN	CTCTCTAC	AAGGAGTA
    H9FT9ADXX-G3686-1_NuGEN	G3686-1_NuGEN	GGACTCCT	ACTGCATA
    H9FT9ADXX-G3676-1_NuGEN	G3676-1_NuGEN	TAGGCATG	AAGGAGTA
    H9FT9ADXX-G3676-2_NuGEN	G3676-2_NuGEN	TAGGCATG	CTAAGCCT
    H9FT9ADXX-G3677-1_NuGEN	G3677-1_NuGEN	TAGGCATG	GTAAGGAG
    H9FT9ADXX-G3677-2_NuGEN	G3677-2_NuGEN	TAGGCATG	ACTGCATA
    H9FT9ADXX-G3687-1_NuGEN	G3687-1_NuGEN	TAGGCATG	TAGATCGC
    Any idea how to fix this? It's driving me potty.

  • #2
    ..... aaaaand after 3.5h I figured it out - my good old friend 'read structure'! I forgot to change 101T8B101T to 101T8B8B101T

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Essential Discoveries and Tools in Epitranscriptomics
      by seqadmin




      The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
      04-22-2024, 07:01 AM
    • seqadmin
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 08:47 AM
    0 responses
    16 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    60 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    60 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 09:21 AM
    0 responses
    54 views
    0 likes
    Last Post seqadmin  
    Working...
    X