Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bowtie2 index error help - beginner

    Hi guys,

    I am a newbie to bioinformatics and having some trouble getting a grip on a Bowtie2 index.

    I have my 2 files (wu_0_A_wgs.fastq and wu_0.v7.fas) located in /home/guest

    when I run the following command.

    bowtie2 -p 4 -x /home/guest/wu_0 wu_0_A_wgs.fastq -S wu.bt2.sam

    I get the following error:
    Could not locate a Bowtie index corresponding to basename "/home/guest/wu_0"
    Error: Encountered internal Bowtie 2 exception (#1)
    Command: /usr/local/bin/bowtie2-align-s --wrapper basic-0 -p 4 -x /home/guest/wu_0 -S wu.bt2.sam wu_0_A_wgs.fastq
    (ERR): bowtie2-align exited with value 1

    I am trying to see how my input is not allowing the index to take place.
    The help shows an input:
    bowtie2 [options]* -x <bt2-idx> {-1 <m1> -2 <m2> | -U <r>} [-S <sam>]
    Just can't see where i'm going wrong.
    Any advice would be amazing!
    Thank you

  • #2
    Is "wu_0.v7.fas" the reference fasta file and the other file actual reads?

    You would need to create an index separately with the reference fasta first. Something like

    Code:
    $ bowtie2-build wu_0.v7.fas wu_indx
    This will create multiple files with "wu_indx" name as prefix. These files constitute the bowtie index.

    You will then align your fastq file against this reference genome index:

    Code:
    $ bowtie2 -p 4 -x wu_indx -U /home/guest/wu_0_A_wgs.fastq -S wu_aligned.sam
    Last edited by GenoMax; 03-30-2016, 05:10 PM. Reason: Correction

    Comment


    • #3
      Thanks for your help...
      yes wu_0.v7.fas is the reference fasta.

      Running that exact code:
      $ bowtie2 -p 4 -x wu_indx -U /home/guest/wu_0 wu_0_A_wgs.fastq -S wu.bt2.sam I get the following error:

      Extra parameter(s) specified: "wu_0_A_wgs.fastq"
      Error: Encountered internal Bowtie 2 exception (#1)
      Command: /usr/local/bin/bowtie2-align-s --wrapper basic-0 -p 4 -x wu_indx -S wu.bt2.sam -U /home/guest/wu_0 wu_0_A_wgs.fastq
      (ERR): bowtie2-align exited with value 1

      Quite confused as to why..
      Thanks again

      Comment


      • #4
        Sorry. My mistake. I had copied and pasted your original command and forgot to take out two extra words that were there. Corrected version is in post #2. Please try that.

        Comment


        • #5
          Thank you heaps for your help!

          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
          28 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 12-13-2024, 08:24 AM
          0 responses
          43 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 12-12-2024, 07:41 AM
          0 responses
          29 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 12-11-2024, 07:45 AM
          0 responses
          42 views
          0 likes
          Last Post seqadmin  
          Working...
          X