Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • 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

  • #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

    • seqadmin
      Advancing Precision Medicine for Rare Diseases in Children
      by seqadmin




      Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...
      12-16-2024, 07:57 AM
    • seqadmin
      Recent Advances in Sequencing Technologies
      by seqadmin



      Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.

      Long-Read Sequencing
      Long-read sequencing has seen remarkable advancements,...
      12-02-2024, 01:49 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 12-17-2024, 10:28 AM
    0 responses
    33 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 12-13-2024, 08:24 AM
    0 responses
    48 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 12-12-2024, 07:41 AM
    0 responses
    34 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 12-11-2024, 07:45 AM
    0 responses
    46 views
    0 likes
    Last Post seqadmin  
    Working...
    X