Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Adapter trimming with cutadapt

    Hi,
    can you please give me an advice. I have paired end sequencing data from amplicon sequencing using TruSeq adaptors (and some barcoding using Myeloid panel Illumina kit).

    I have noticed using FastQC that there are many reads that have contamination by Illumina universal adapter in R2 but not in R1. I quite do not understand the fact that there is such disproportion. So maybe I am missing something such as reverse complement the reads from R1 or adapter sequence.

    However knowing that adapter sequence starts with AGATCGGAAGAG (I guess for both adapters) is it sufficient to remove adapters by this command?

    cutadapt -a AGATCGGAAGAG -A AGATCGGAAGAG -o out.1.fastq -p out.2.fastq IN_R1_001.fastq IN_R2_001.fastq

    Thank you,
    Vojtech.
    Data analysis www.persmed.eu

  • #2
    Adapter sequence for fragment libraries should always appear at the exact same location in R1 as in R2 for a given pair. Thus adapter-trimmed paired reads should come out with R1 and R2 being the same length.

    If you are not sure what your adapter sequence is, you can find out with BBMerge (post #48). 12bp is kind of short for adapter removal, particularly if mismatches are allowed.

    Comment


    • #3
      Thank you,
      you pointed me out in a good direction.

      Acording to extensive analysis of supplied files I have found out adapter sequences:
      Read1 :GCGAATTTCGACGATCGTTGCATTAACTCGCGA
      Read2 :AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT

      resulting in a command
      cutadapt -a GCGAATTTCGACGATCGTTGCATTAACTCGCGA -A AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT -o /dev/null -p /dev/null R1.fastq R2.fastq

      I am also considering using trimmomatic.

      I have barcoded sequences. And the automatic detection would be good. The problem with my data is that Illumina (or my sequencing facility in their name) refused to confirm adapter sequences. I know that the data are coming from http://www.illumina.com/products/trusight-myeloid.html and are barcoded using Illumina barcodes. Adaptor contamination is present only in a proportion of my files. In some files they are not present.

      I promise to test my data with your tool if the detection is correct. The problem is that I do not understand how to supply first pair reads and how second pair reads. And how your tool will figure out where to put N's for barcode when I supply only one sample barcoded with only one barcode. I have approximately 96 samples sequenced in one run, multiplexed. All samples (same library) were sequenced on the two lanes as paired end experiment. Thus I have 4 files for each sample. If you are interested in testing your program on these data I am interested. PM me.

      Another Issue will be to remove primer sequences from the reads and then it would be nice to use FASTX or your software to cover whole amplicon with a read. But I doubt that using such tools is of no benefit when primer sequences are present in the data.

      Vojtech.
      Data analysis www.persmed.eu

      Comment


      • #4
        Vojtech: Both BBDuk (part of BBMap) and trimmomatic include sequence files for standard illumina adapters.

        With BBDuk if you're not sure which adapters are used, you can add "ref=truseq.fa.gz,truseq_rna.fa.gz,nextera.fa.gz" to your command line and get them all (this will increase the amount of overtrimming, though it should still be negligible).

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Best Practices for Single-Cell Sequencing Analysis
          by seqadmin



          While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
          06-06-2024, 07:15 AM
        • seqadmin
          Latest Developments in Precision Medicine
          by seqadmin



          Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

          Somatic Genomics
          “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
          05-24-2024, 01:16 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Yesterday, 07:49 AM
        0 responses
        12 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 06-20-2024, 07:23 AM
        0 responses
        14 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 06-17-2024, 06:54 AM
        0 responses
        16 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 06-14-2024, 07:24 AM
        0 responses
        25 views
        0 likes
        Last Post seqadmin  
        Working...
        X