Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Kaskere
    Junior Member
    • Jul 2018
    • 7

    #1

    bwa alignement

    Hello,
    I am trying to align viral RNA sequences against reference fasta file.
    It looks like it is not working properly

    my index code looks like that

    ./bwa index reference.fasta

    but I keep getting output like

    [M::bwa_idx_load_from_disk] read 0 ALT contigs
    [M:rocess] read 136386 sequences (10000131 bp)...
    [M:rocess] read 136392 sequences (10000087 bp)...
    [M::mem_pestat] # candidate unique pairs for (FF, FR, RF, RR): (0, 24370, 7, 0)
    [M::mem_pestat] skip orientation FF as there are not enough pairs
    [M::mem_pestat] analyzing insert size distribution for orientation FR...
    [M::mem_pestat] (25, 50, 75) percentile: (138, 175, 224)
    [M::mem_pestat] low and high boundaries for computing mean and std.dev: (1, 396)
    [M::mem_pestat] mean and std.dev: (186.45, 62.40)
    [M::mem_pestat] low and high boundaries for proper pairs: (1, 482)
    [M::mem_pestat] skip orientation RF as there are not enough pairs
    [M::mem_pestat] skip orientation RR as there are not enough pairs
    [M::mem_process_seqs] Processed 136386 reads in 3.278 CPU sec, 3.070 real sec

    Should I include options like -p and -a? I have no idea what exactly should be written in "Prefix of the output database" after index command
  • GenoMax
    Senior Member
    • Feb 2008
    • 7142

    #2
    What is not working properly? Indexing genome and alignment are two independent steps. You completed the indexing of the genome. Now you need to align using `bwa mem`+ this index + your reads => Aligned SAM file => Aligned BAM file.

    You could give your database a name that is what the "prefix" refers to. If you don't provide a name then your index is going to use the same "basename" as your reference file. In this case, "reference.fasta".

    Comment

    Latest Articles

    Collapse

    • SEQadmin2
      New Genomics Technologies Take Aim at Long-Standing Limits
      by SEQadmin2


      Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

      We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
      ...
      09-28-2026, 10:25 AM
    • SEQadmin2
      How Immunogenomics Decodes Immunity’s Genetic Blueprint
      by SEQadmin2




      The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.

      This convergence of genetics, immunology, and computation...
      09-01-2026, 05:41 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, Yesterday, 09:51 AM
    0 responses
    13 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 09-25-2026, 09:06 AM
    0 responses
    37 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 09-23-2026, 11:05 AM
    0 responses
    28 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 09-18-2026, 11:37 AM
    1 response
    48 views
    0 reactions
    Last Post pekgio
    by pekgio
     
    Working...