Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • GenoMax
    replied
    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).

    Leave a comment:


  • Kulvait
    replied
    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.

    Leave a comment:


  • Brian Bushnell
    replied
    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.

    Leave a comment:


  • Kulvait
    started a topic Adapter trimming with cutadapt

    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.

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
9 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