Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Sequence Trimming without Knowing Primer Sequences

    Hi,

    I have an IonTorrent output but the provider doesn't want to give information about primer sequences. Instead it gives a bed file includes regions of interest. So if I do alignment and bedtools intersect than trimming sequences by outside parts which is not in the region, I'm not comfortable this is a good approach. We want to use the output as a data for finding SNP frequencies and for the best approach primer sequences must be trimmed.

    So, is there a way to do this job properly? I mean is there a way to trim sequences with no knowledge about primers?

  • #2
    What kind of experiment are you doing? Is it some kind of exon-capture, or an amplicon approach? Or, why do you expect to see primers in the sequence?

    Comment


    • #3
      Originally posted by Brian Bushnell View Post
      What kind of experiment are you doing? Is it some kind of exon-capture, or an amplicon approach? Or, why do you expect to see primers in the sequence?
      It's basically a variant analysis process. What I need is trimming sequences by primer sequences but we have no information a about primer sequences.

      Comment


      • #4
        Use fastq contaminant detection and removal approaches. Do a bit of de novo...

        If you would like to have some idea about the primers you can try:
        Converting to fastq format and having a look at fastqc at high kmer settings (-k8 or -k10).

        Up the memory to java wm (-Xmx) if you run out of memory for the fastqc application.

        Also try assembling a susbset (1k - 10k) of reads matching the primer-like kmers, and having a look at the de novo assembly results in consed or similar program.
        Newbler assembler has primer detection capability and may give you some idea about the primers present (check it's log).

        Local setup of WWWblast / STARblast or similar tool comes very handy for doing investigations.

        grep -c primer *.fasta against the sequences in fasta format may give you an idea about the abundance of the given kmer.

        Other options to investigate includes trimmomatic and similar tools.

        Also look for high quality unaligned region a the the ends of your reads when examining the mapping results.

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