Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #16
    Originally posted by michfish86 View Post
    Edit: Could the reason for the lack of matching to the reference transcriptome just be the presence of novel transcripts in my samples that were not present in the individuals in the other study that generated the reference? For instance, because the reference may not have included the specific type of tissue that I'm evaluating?
    While that is tempting do you have reasons to suspect that may be case. Is the genome odd (large with multiple chromosomes, ploidy difference)?

    Even though you may not have enough data you could try a de novo assembly and see what you get. Be sure to map the reads back to your own assembly to see what percentage of them map?

    Comment


    • #17
      Originally posted by GenoMax View Post
      Where do you see the 20% map?
      I think Macspider's phrasing a just a little confusing. I think s/he means that more than 20% of my reads are not mapping.

      Comment


      • #18
        Originally posted by GenoMax View Post
        While that is tempting do you have reasons to suspect that may be case. Is the genome odd (large with multiple chromosomes, ploidy difference)?
        Looking back at the paper that accompanies the reference transcriptome, they used several different tissue types, but did not include the exact tissue that I sequenced. The genome isn't particularly huge (C-value 4.75), but the chromosomal arraignments have been tough to nail down. One study that looked into it found a diploid number of 110 with a range of 94-185.

        Even though you may not have enough data you could try a de novo assembly and see what you get. Be sure to map the reads back to your own assembly to see what percentage of them map?
        I will give it a shot.

        Thanks again!!

        Comment


        • #19
          BBMap suite has seal.sh and tadpole.sh. You may want to look at them as substitute assemblers (though trinity has been the assembler of choice for eukaryotes of late).

          Comment


          • #20
            While my de novo assembly runs, I wanted to get clarification on a what could be causing the large difference in alignment percentage between Bowtie and Bowtie2 and BBMap. Reminder of alignment scores below:

            Bowtie: 43%
            Bowtie2: 71%
            BBMap: 78%

            I keep seeing that Bowtie2 and BBMap are "more sensitive" or "more flexible" but I can't figure out what that means and I'm skeptical of the higher alignment scores because I'm afraid of including a bunch of false positives.

            Two potential reasons that I'm aware of for Bowtie being substantially lower is because Bowtie doesn't allow ambiguous characters (e.g., N) in alignments and it disqualifies discordant reads, whereas Bowtie2 still seeks alignments despite these imperfections. I'd appreciate any additional insight in what could be causing these disparities since I don't want to blindly accept the higher alignment scores just because they're higher!

            Edit: Update: Running Bowtie2 with the --no-mixed option turned on reduces the alignment percent to 54%. It appears that Bowtie is "--no-mixed" by default, therefore this could be a major contributor to the difference between my Bowtie and Bowtie2 percentages. In "--no-mixed" mode, alignments must involve both pairs of reads, therefore, in Bowtie I may be losing some reads just because only one of the reads constituting a "pair" aligns. Again, I'd appreciate anyone else's insights.
            Last edited by michfish86; 08-25-2016, 12:12 PM.

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Quality Control Essentials for Next-Generation Sequencing Workflows
              by seqadmin




              Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

              Nucleic Acid Quality Control
              Preparing for NGS starts with isolating the...
              02-10-2025, 01:58 PM
            • seqadmin
              An Introduction to the Technologies Transforming Precision Medicine
              by seqadmin


              In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
              01-27-2025, 07:46 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 02-07-2025, 09:30 AM
            0 responses
            54 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 02-05-2025, 10:34 AM
            0 responses
            88 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 02-03-2025, 09:07 AM
            0 responses
            70 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 01-31-2025, 08:31 AM
            0 responses
            44 views
            0 likes
            Last Post seqadmin  
            Working...
            X