Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • vl80
    Member
    • Nov 2012
    • 34

    #1

    Bowtie on Galaxy - multiple amplicon sequences instead of reference genome

    I am trying to use Bowtie in Galaxy via the web GUI to align Miseq generated amplicons. Our amplicons are from a multiplex PCR and thus instead of a single reference sequence we have a number of (27 to be exact) short sequences (250 - 400 bp) as the references. Now the products can align against any of these, so in the ideal situation I need to be able to upload a sigle file with all the different sequenes and then after alignment generate a table to say how many reads were aligned agaist which sequence. I tried to do this but get a failure message in Galaxy. Is it possible for this to be done or do I need to creat a file where all the sequences are strined together?
    Thanks a advance.
  • GenoMax
    Senior Member
    • Feb 2008
    • 7142

    #2
    Did you do it this way? http://wiki.galaxyproject.org/Learn/CustomGenomes

    Comment

    • vl80
      Member
      • Nov 2012
      • 34

      #3
      My ref sequences are in FASTA formates with compatible lengths per line
      However the file is something like

      >Ref_amplicon1
      AAAACCGGTTTAA
      >Ref_amplicon2
      GGTTCCAAAAAA
      >Ref_Amplicon3
      CCCTTGGAAA

      I am not sure if this would be an acceptable form

      Comment

      • GenoMax
        Senior Member
        • Feb 2008
        • 7142

        #4
        That format should be fine.

        The failure message is for what step? Is the alignment failing or something else?

        Comment

        • vl80
          Member
          • Nov 2012
          • 34

          #5
          Actually, alignment went OK, my mistake in posting.
          When I wanted to convert SAM-BAM an error came. Also, I tried to filter the alined SAM data, the number of sequences remaining in the end were very low by just using the paird reads flag only. Our service provider used bow tie script based method and gave summary tables to us, which showed a high percentage of aligned reads!, Not sure where things go wrong.

          Comment

          • vl80
            Member
            • Nov 2012
            • 34

            #6
            This is the error I get on trying to get the coverage stats...

            An error occurred with this dataset: Picked up _JAVA_OPTIONS: -Djava.io.tmpdir=/space/g2main [Thu Apr 04 12:10:28 EDT 2013] net.sf.picard.sam.CreateSequenceDictionary REFERENCE=/space/g2main/tmp-gatk-dbNNN9/gatk_input.fasta OUTPUT=/space/g2main/tmp-gatk-dbNNN9/dict7890443155236506554.tmp

            Comment

            • GenoMax
              Senior Member
              • Feb 2008
              • 7142

              #7
              Are you using the main galaxy site or a local mirror?

              Comment

              • vl80
                Member
                • Nov 2012
                • 34

                #8
                I am using the Galaxy site at https://main.g2.bx.psu.edu/ Not a locally downloaded version..

                Comment

                • GenoMax
                  Senior Member
                  • Feb 2008
                  • 7142

                  #9
                  Have you posted to the galaxy help list ([email protected])? They are very helpful and prompt to respond.

                  Comment

                  • vl80
                    Member
                    • Nov 2012
                    • 34

                    #10
                    Yes I did, unfortunately no reply yet. I tried to send an error message too. Still no feed back.

                    Comment

                    Latest Articles

                    Collapse

                    • SEQadmin2
                      Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
                      by SEQadmin2



                      Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

                      Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
                      ...
                      Yesterday, 09:55 AM
                    • 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

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by SEQadmin2, Today, 10:52 AM
                    0 responses
                    8 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, Yesterday, 05:39 AM
                    0 responses
                    17 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 10-13-2026, 09:56 AM
                    0 responses
                    15 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 10-09-2026, 11:14 AM
                    0 responses
                    27 views
                    0 reactions
                    Last Post SEQadmin2  
                    Working...