Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Best tool for the job (noob question)

    What are the best tools for aligning a query against a target with the following characteristics? I had previously tried bowtie but only got 6% alignment when over 90% was expected:

    Query file:
    > Format: multisequence fasta
    > Number of sequences: 180,099
    > Minimum sequence length: 69 bases
    > Maximum sequence length: 1,045 bases
    > Average sequence length: 825 bases
    > Number of Ns in file: 158,231
    > Longest N stretch: 363
    > Number of sequences containing N: 59,998
    Target file:
    > Number of sequences: 1,168
    > Minimum sequence length: 1,028 bases
    > Maximum sequence length: 63,086,992 bases
    > Average sequence length: 843,760 bases
    > Number of Ns in file: 18,289,543
    > Longest N stretch: 1,185,460
    > Number of sequences containing N: 771

    If anyone could point me in the right direction it would be greatly appreciated. Thanks!
    Last edited by bre; 01-05-2010, 10:43 AM.

  • #2
    Originally posted by bre View Post
    What are the best tools for aligning a query against a target with the following characteristics? I had previously tried bowtie but only got 6% alignment when over 90% was expected:

    Query file:
    > Format: multisequence fasta
    > Number of sequences: 180,099
    > Minimum sequence length: 69 bases
    > Maximum sequence length: 1,045 bases
    > Average sequence length: 825 bases
    > Number of Ns in file: 158,231
    > Longest N stretch: 363
    > Number of sequences containing N: 59,998
    Target file:
    > Number of sequences: 1,168
    > Minimum sequence length: 1,028 bases
    > Maximum sequence length: 63,086,992 bases
    > Average sequence length: 843,760 bases
    > Number of Ns in file: 18,289,543
    > Longest N stretch: 1,185,460
    > Number of sequences containing N: 771

    If anyone could point me in the right direction it would be greatly appreciated. Thanks!
    How about trying BLAT? It may be slow, but in this situation it may give you proper sensitivity.

    Comment


    • #3
      Would bwa work? I need gapped alignments which bowtie cannot do, and I'd also like to get SAM-format output.

      Comment


      • #4
        Originally posted by bre View Post
        Would bwa work? I need gapped alignments which bowtie cannot do, and I'd also like to get SAM-format output.
        Since the query is up to 1Kb I would try the "bwa bwasw" command, though I have not tried it. You might have to split your queries into long and short ones. Anyone else have any advice?

        Comment


        • #5
          It looks right now like bfast may be my best bet. I noticed bwa doesn't allow paired-end reads for sequences longer than ~200bp.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Non-Coding RNA Research and Technologies
            by seqadmin




            Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

            Nobel Prize for MicroRNA Discovery
            This week,...
            Yesterday, 08:07 AM
          • seqadmin
            Recent Developments in Metagenomics
            by seqadmin





            Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
            09-23-2024, 06:35 AM
          • seqadmin
            Understanding Genetic Influence on Infectious Disease
            by seqadmin




            During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.

            Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...
            09-09-2024, 10:59 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 10-02-2024, 04:51 AM
          0 responses
          95 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 10-01-2024, 07:10 AM
          0 responses
          106 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 09-30-2024, 08:33 AM
          1 response
          106 views
          0 likes
          Last Post EmiTom
          by EmiTom
           
          Started by seqadmin, 09-26-2024, 12:57 PM
          0 responses
          20 views
          0 likes
          Last Post seqadmin  
          Working...
          X