Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Understanding paired-end SOLiD

    Hello everyone,
    I am working with paired-end SOLiD data, and I am not sure how to handle it.

    1) What exactly is the difference between the F3 and F5 reads? Is the F5 reversed?

    2) I'm using BWA for the analysis. How do I use it with paired-end? In the documentation it says 'file 1' and 'file 2', but I don't know which is which.
    I have tried to run it both ways and in both times I got very few variants - 3000, 5000 while expecting 50000.
    I have a color-space reference index, and I used the -c parameter in the alignment. Am I missing something?

    3) Does the csfasta conversion work differently with paired-end data?

    Thanks a lot,
    Amit

    =]

  • #2
    1) What exactly is the difference between the F3 and F5 reads? Is the F5 reversed?
    They are strand-specific reads. For the same RNA sequence, they will map in different directions to each other.

    I think the 5' read is reverse and done first, so if you're doing an analysis taking reads in the usual order you need to specify RF orientation in the mapping program (I'm not very sure on this point).

    3) Does the csfasta conversion work differently with paired-end data?
    Nope. You shouldn't do this conversion with either. Map in colour-space and let the mapping tool give you the equivalent base-space sequence.

    Comment


    • #3
      Thank you very much!

      Do you run such analyses with BWA?

      Comment


      • #4
        No, I haven't [yet] bothered with looking at another program. While I might get a higher proportion of mapped reads, I've got enough problems working with error in the reads I have from Bowtie.

        Comment


        • #5
          But... Bowtie does not handle color-space... or does it?

          Comment


          • #6
            Bowtie can do color-space, Bowtie2 cannot (and is not likely to in the future).

            Comment


            • #7
              ~bump~
              Does anyone work with BWA on SOLiD paired-end data?

              Comment


              • #8
                BWA does not support paired end reads with SOLiD.

                Comment


                • #9
                  I know the new version does not work with color-space.
                  Are you saying that it's true for the older versions too?

                  Comment


                  • #10
                    The new version does not work with color space. The older version works with mate pairs and color space, but not paired end and color space.

                    Comment


                    • #11
                      How do you know it's meant for mate-pairs?

                      On the documentation it says:
                      Generate alignments in the SAM format given paired-end reads.

                      Comment


                      • #12
                        Originally posted by AmitL View Post
                        How do you know it's meant for mate-pairs?

                        On the documentation it says:
                        Check out the source code "bwape.c" and search for SOLID. It uses the F3/R3 "pairing", which is mate pair.

                        Comment


                        • #13
                          Alright. Thanks a lot!

                          Comment


                          • #14
                            Have you consider to map your csfasta reads with SHRiMP2?

                            Comment


                            • #15
                              Yes.
                              Shrimp2 collapses because of non-existing input problems. I have checked the input manually and it was ok. It works with all other tools.

                              I preferred not to deal with such a buggy tool.

                              But thanks!

                              Comment

                              Latest Articles

                              Collapse

                              • seqadmin
                                Exploring the Dynamics of the Tumor Microenvironment
                                by seqadmin




                                The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
                                07-08-2024, 03:19 PM
                              • seqadmin
                                Exploring Human Diversity Through Large-Scale Omics
                                by seqadmin


                                In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
                                06-25-2024, 06:43 AM

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by seqadmin, 07-10-2024, 07:30 AM
                              0 responses
                              30 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 07-03-2024, 09:45 AM
                              0 responses
                              202 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 07-03-2024, 08:54 AM
                              0 responses
                              212 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 07-02-2024, 03:00 PM
                              0 responses
                              194 views
                              0 likes
                              Last Post seqadmin  
                              Working...
                              X