Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • lcmb
    Junior Member
    • May 2015
    • 6

    Comparing Merged Seqs to a Reference Seq

    Hello,

    I recently received paired-end reads from an amplicon library I sequenced with Illumina. I used BBMerge to merge the paired-ends and, now, I would like to compare my merged sequences to a reference. My goal is to analyze where the merged sequence may differ from the reference.

    I have tried to write a small Perl program, but I am continually running into problems. I study DNA repair and, essentially, I would like to compare my merged sequences to the reference so that I may better understand repair efficiency.

    Is anyone aware of any programs/scripts that could help me reach my goal of comparing my merges to a reference? They are about 200bp long. Thank you!
  • GenoMax
    Senior Member
    • Feb 2008
    • 7142

    #2
    BBMap.sh (where you got BBMerge from) will do the comparisons/alignments.

    Added: Is your reference 200 bp long?

    Comment

    • lcmb
      Junior Member
      • May 2015
      • 6

      #3
      Originally posted by GenoMax View Post
      BBMap.sh (where you got BBMerge from) will do the comparisons/alignments.

      Added: Is your reference 200 bp long?
      Yes, it is. Thank you, GenoMax! Will BBMap allow me to set any criteria for the alignment? For example, if I did not want any sequences that differ by more than 10% from my reference.

      Comment

      • GenoMax
        Senior Member
        • Feb 2008
        • 7142

        #4
        Run bbmap.sh on its own to see all the command line options.

        Here are some relevant options

        Post-Filtering Parameters:

        idfilter=0 Independant of minid; sets exact minimum identity
        allowed for alignments to be printed. Range 0 to 1.
        subfilter=-1 Ban alignments with more than this many substitutions.
        insfilter=-1 Ban alignments with more than this many insertions.
        delfilter=-1 Ban alignments with more than this many deletions.
        indelfilter=-1 Ban alignments with more than this many indels.
        editfilter=-1 Ban alignments with more than this many edits.
        inslenfilter=-1 Ban alignments with an insertion longer than this.
        dellenfilter=-1 Ban alignments with a deletion longer than this.
        Last edited by GenoMax; 05-01-2015, 10:37 AM.

        Comment

        • GenoMax
          Senior Member
          • Feb 2008
          • 7142

          #5
          You could also convert your reads to fasta format (use reformat.sh from BBMap) and then use blat for the alignments.

          Once you identify the reads that align I suppose you want to do a multiple sequence alignment?

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Pathogen Surveillance with Advanced Genomic Tools
            by seqadmin




            The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
            03-24-2025, 11:48 AM
          • seqadmin
            New Genomics Tools and Methods Shared at AGBT 2025
            by seqadmin


            This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

            The Headliner
            The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
            03-03-2025, 01:39 PM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 03-20-2025, 05:03 AM
          0 responses
          41 views
          0 reactions
          Last Post seqadmin  
          Started by seqadmin, 03-19-2025, 07:27 AM
          0 responses
          49 views
          0 reactions
          Last Post seqadmin  
          Started by seqadmin, 03-18-2025, 12:50 PM
          0 responses
          36 views
          0 reactions
          Last Post seqadmin  
          Started by seqadmin, 03-03-2025, 01:15 PM
          0 responses
          191 views
          0 reactions
          Last Post seqadmin  
          Working...