Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • carmeyeii
    replied
    Hi!

    I'm analyzing a "second-hand" dataset generated using SOLiD 4. It is a transcriptome mate pair library that is 52 x 37 nt, and I cannot for the sake of me find the protocol that was used to generate those specific read lengths. I have F3 and R3 reads, so I am assuming it is a circularization protocol, but I do not know what the size selection parameters were, or how the circles were cut to produce the final fragments. This info would be very valuable for a more accurate mapping.

    Any knowledge would be greatly appreciated!

    Thanks a lot,

    Carmen

    Leave a comment:


  • carmeyeii
    replied
    I'd like to bump this up as I'm trying to map some SOLID 4 reads. Thanks!

    Leave a comment:


  • mrawlins
    replied
    SOLiD Paired-end with Bowtie

    This post is really helpful, addressing Illumina mate-pair and paired-end as well as the SOLiD mate-pair. The newest SOLiD machines also have a paired-end, which as near as I can tell is of the configuration:

    Code:
    read1                           read2
    -------------->                 <--------------
    ===============================================
    Now, since the reverse compliment of color-space is the just the reverse sequence, it seems to me this is equivalent to the -fr from bowtie. Has anybody run Bowtie with the SOLiD Paired-end data?

    Leave a comment:


  • days369
    replied
    kmcarr, thanks so much. your reply is really helpful.

    I guess that my reads should be produced using the Illumina paired end protocol, so most of them are in --fr orientation.

    But a small portion of the paired-end reads are in --rf/--ff orientation, because they can be mapped with those two options. Some are even in an orientation as illustrated below. I guess they are different from the -ff orientation, because read2 is 5' to read1.

    Do you have any idea why the same illumina paired-end protocol can generate some reads that is not in -fr orientation? or Do you meet this problem before?

    read2 read1
    --------------> -------------->
    ===============================================

    OR

    ===============================================
    <-------------- <--------------
    read1 read2

    Leave a comment:


  • kmcarr
    replied
    Dezhi,

    There is no --rr because that is the same as --ff, molecularly speaking. The 'f' and 'r' do not refer to the strand which the read aligns to, but the relative orientation of read1 and read2. The relative orientation of the reads 1 and 2 is determined by the method used to generate the paired reads; it should be identical for all read pairs in your data set and you should now what the expected orientation is based on the method used to construct the paired end library. (Let's ignore for now that unexpected pair orientations can signal structural variations.)

    The --fr orientation means that the two reads are aligned to opposite strands and are pointed toward each other (i.e. the 3' ends of the reads are closer together than the 5' ends).

    Code:
    -------------->
    ===============================================
                                    <--------------
    The --rf orientation means that the two reads are aligned to opposite strands and are pointed away from each other (i.e. the 5' ends of the reads are closer together than the 3' ends).

    Code:
                                    -------------->
    ===============================================
    <--------------
    The --ff orientation means that both reads are aligned to the same strand and are pointed in the same direction. Further, read1 should always be 5' to read2, relative to the strand they are aligned to. In the two cases above you should have an equal distribution of read1 being the read "on the left".

    Code:
    read1                           read2
    -------------->                 -------------->
    ===============================================
    
    OR
    
    ===============================================
    <--------------                 <--------------
              read2                            read1
    Read pairs in the --fr orientation are produced using the Illumina paired end protocol. Read pairs in the --rf orientation are produced using the Illumina mate-pair protocol. Read pairs in the --ff orientation are produced in using the SOLiD mate-pair protocol.

    Leave a comment:


  • days369
    started a topic Bowtie Paired-end options

    Bowtie Paired-end options

    Hi everyone,

    I just began to use bowtie to map my solexa paired-end reads to a genome. I found that bowtie has three options to specify the orientation of the paired-end reads, --fr, --rf, and --ff. I tried all of them, and found:

    --fr, as the default one, can map most of my data onto the genome. For the left ones, --rf and --ff options can map some, but still leave some reads, which can be separately mapped in non-paired-end mode(should I call it single-end mode), but cannot be recognized as Paired-end reads with any of the three options.

    Carefully examining those remaining reads, I found most of them are in an orientation of rev-rev. In another word, both ends of each pair are mapped to the reverse strand. Why does bowtie have no --rr option? How can I process those reads in paired-end mode?

    THanks a lot!
    Dezhi

Latest Articles

Collapse

  • seqadmin
    Best Practices for Single-Cell Sequencing Analysis
    by seqadmin



    While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
    06-06-2024, 07:15 AM
  • seqadmin
    Latest Developments in Precision Medicine
    by seqadmin



    Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

    Somatic Genomics
    “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
    05-24-2024, 01:16 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 07:24 AM
0 responses
10 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-13-2024, 08:58 AM
0 responses
11 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-12-2024, 02:20 PM
0 responses
16 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-07-2024, 06:58 AM
0 responses
184 views
0 likes
Last Post seqadmin  
Working...
X