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
          Exploring the Dynamics of the Tumor Microenvironment
          by seqadmin




          The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
          07-08-2024, 03:19 PM
        • seqadmin
          Exploring Human Diversity Through Large-Scale Omics
          by seqadmin


          In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
          06-25-2024, 06:43 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Yesterday, 05:49 AM
        0 responses
        12 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-15-2024, 06:53 AM
        0 responses
        23 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-10-2024, 07:30 AM
        0 responses
        36 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-03-2024, 09:45 AM
        0 responses
        204 views
        0 likes
        Last Post seqadmin  
        Working...
        X