Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • comparaison to the reference genome

    Hello,
    I generated reads a file (.fasta) I want to test the effectiveness of my file to the reference genome (before assembly).
    how to do?
    and what tool I have to use it to make the comparison of my file with the reference genome?
    Thank you

  • #2
    I suggest using Quast; it will compare your assembly to the reference. If you mean that you want to compare reads to an assembly, then you should use mapping, and a tool like ALE.

    Comment


    • #3
      Originally posted by Brian Bushnell View Post
      I suggest using Quast; it will compare your assembly to the reference. If you mean that you want to compare reads to an assembly, then you should use mapping, and a tool like ALE.
      to compare read file to the reference genome i should use mapping tools?
      like ALE?
      thanks

      Comment


      • #4
        If you want to compare reads, then map them, and then use ALE.

        Comment


        • #5
          hello,
          i used BBMAP for mapping and this is the results:
          Code:
             ------------------   Results   ------------------
          
          Genome:                 1
          Key Length:             13
          Max Indel:              16000
          Minimum Score Ratio:    0.56
          Mapping Mode:           normal
          Reads Used:             1626223 (793501894 bases)
          
          Mapping:                21520.101 seconds.
          Reads/sec:              75.57
          kBases/sec:             36.87
          
          
          Read 1 data:            pct reads       num reads       pct bases          num bases
          
          mapped:                  84.1693%         1368781        84.7943%          672844516
          unambiguous:             75.6937%         1230949        76.2771%          605259839
          ambiguous:                8.4756%          137832         8.5173%           67584677
          low-Q discards:           0.0000%               0         0.0000%                  0
          
          perfect best site:       72.2714%         1175294        73.8835%          586267084
          semiperfect site:        72.2780%         1175402        73.8902%          586320010
          
          Match Rate:                   NA               NA        94.7978%          668713525
          Error Rate:              13.8973%          193385         5.1976%           36664236
          Sub Rate:                12.5303%          174364         0.2647%            1867550
          Del Rate:                 6.0770%           84563         4.6165%           32565546
          Ins Rate:                 5.8459%           81348         0.3163%            2231140
          N Rate:                   0.0630%             877         0.0046%              32301
          do i have good results or not?
          and can explain these lines?
          thanks

          Comment


          • #6
            any response?

            Comment


            • #7
              First off, it is very important to clearly explain what you did. Mapping reads to their assembly is different than mapping to a reference; which did you do? And if you have a reference, what are you making reads and assemblies for? Is this a different strain, different species, or what? What kind of data do you have, from which platform?

              Those numbers are very useful when trying to determine which is a better assembly of two assemblies from the same reads, or which is a better set of reads when you have two libraries made from the same organism. However, it's hard to determine the absolute quality of an assembly just from the mapping results, since you don't know whether the mismatches are caused by misassemblies or low-quality data, or how much contamination there is, etc. The fact that 72.27% of the reads match the assembly perfectly is nice. But the fact that only 84.17% of the reads map to the assembly indicates that it is probably incomplete.

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Strategies for Sequencing Challenging Samples
                by seqadmin


                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                03-22-2024, 06:39 AM
              • seqadmin
                Techniques and Challenges in Conservation Genomics
                by seqadmin



                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                Avian Conservation
                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                03-08-2024, 10:41 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 03-27-2024, 06:37 PM
              0 responses
              13 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-27-2024, 06:07 PM
              0 responses
              11 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              53 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-21-2024, 07:32 AM
              0 responses
              69 views
              0 likes
              Last Post seqadmin  
              Working...
              X