Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • mido1951
    Senior Member
    • May 2014
    • 123

    #76
    Code:
    for example we have two fragments.
    S1: ATCGTTGAGCAGACT and S2: TGAGCAGACTTAAGTAGTTTT they have an overlap zone"TGAGCAGACT". so the consensus is:ATCGTTGAGCAGACTTAAGTAGTTTT
    
    we sequenced this two fragments S1 and S2 in paired end reads R1 and R2.
    S1:R1: ATCGTTGAGCA
    S1:R2: AGTCTGCTCAA (illumina read the fragment in the other way so I think that to write in reverse complement. correct me if I'm wrong)
    S2:R1:TGAGCAGACTT
    S2:R2:AAAACTACTTA (reverse complement)I think this way is built the R2.
    
    how to make assembly in this case?
    I speak of here. How the assembly in this case?
    in both files (paired end), is what we must do the merging of the two files?

    Comment

    • OTU
      Member
      • May 2013
      • 44

      #77
      You realize that you are not doing this manually, right? You usually use a program, which creates overlapping/merging reads and gives you the contig.

      Comment

      • westerman
        Rick Westerman
        • Jun 2008
        • 1104

        #78
        Use Panda or Flash (or probably a number of other programs) to do merging.

        Comment

        • mido1951
          Senior Member
          • May 2014
          • 123

          #79
          yes I know.
          I'm talking of phase before getting contigs.
          you know to make an assembly of reads that it must overlap.
          In the case of paired end reads, how to find overlaps between the two files?
          we F1.fq that contains reads RX.1 and F2.fq that contains reads RX.2.
          how to find the overlap between these reads? is what we should do a merging between RX.1 and RX.2?
          So the merging is necessary?
          thanks

          Comment

          • GenoMax
            Senior Member
            • Feb 2008
            • 7142

            #80
            @Rick/@OTU: In order to avoid a re-hash of things that have been already discussed in other threads I am posting two below.

            Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

            Bridged amplification & clustering followed by sequencing by synthesis. (Genome Analyzer / HiSeq / MiSeq)


            @mido1951 is either mis-understanding some basic concepts about sequencing/assembly or I am not able to understand what @mid1951 wants to know.

            Comment

            • Brian Bushnell
              Super Moderator
              • Jan 2014
              • 2709

              #81
              Originally posted by mido1951 View Post
              yes I know.
              I'm talking of phase before getting contigs.
              you know to make an assembly of reads that it must overlap.
              In the case of paired end reads, how to find overlaps between the two files?
              we F1.fq that contains reads RX.1 and F2.fq that contains reads RX.2.
              how to find the overlap between these reads? is what we should do a merging between RX.1 and RX.2?
              So the merging is necessary?
              thanks
              It's very simple. You do this:

              bbmerge.sh in1=F1.fq in2=F2.fq out=merged.fq outu=unmerged.fq

              Comment

              • mido1951
                Senior Member
                • May 2014
                • 123

                #82
                thanks for your answer.
                So, the merging is necessary ??
                thanks

                Comment

                • Brian Bushnell
                  Super Moderator
                  • Jan 2014
                  • 2709

                  #83
                  No, merging is optional.

                  Comment

                  • mido1951
                    Senior Member
                    • May 2014
                    • 123

                    #84
                    if the merging is optional,
                    how to do overlap reads without merging the two files into one file??

                    Comment

                    Latest Articles

                    Collapse

                    • GATTACAT
                      Reply to Nine Things a Sample Prep Scientist Thinks About Before Sequencing
                      by GATTACAT
                      Love this - good data definitely starts from good input, and poor input can only give relatively poor data. I particularly like the mention of Nanodrop/absorbance based methods for quantification. It's such a toss up if you'll get an accurate reading or what amounts to a randomly generated number, and a lot of library/sequencing related issues can be traced back to poor quant.
                      07-01-2026, 11:43 AM
                    • SEQadmin2
                      Nine Things a Sample Prep Scientist Thinks About Before Sequencing
                      by SEQadmin2


                      I’m not a sequencing expert. I’m a purification scientist who uses NGS to evaluate workflows my group develops. With this perspective, we think about the sample first and the NGS workflow second. The sequencer is an exceptionally honest reporter, but it can only report on what you give it, so whether you get clean, interpretable data from an NGS workflow is largely determined before you begin.

                      Here are nine questions we think about, in roughly the order they matter, before...
                      06-18-2026, 07:11 AM

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by SEQadmin2, 07-02-2026, 11:08 AM
                    0 responses
                    25 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 06-30-2026, 05:37 AM
                    0 responses
                    23 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 06-26-2026, 11:10 AM
                    0 responses
                    23 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 06-17-2026, 06:09 AM
                    0 responses
                    55 views
                    0 reactions
                    Last Post SEQadmin2  
                    Working...