Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • FASTQ to SAM conversion

    What is the best program to use for converting fastq (or eland extended) files to SAM format? Thanks!

  • #2
    Originally posted by albrown415 View Post
    What is the best program to use for converting fastq (or eland extended) files to SAM format? Thanks!
    Do you want to align the data first, or just represent the FASTQ data in the SAM format? Is the data paired end (mate-pair)?

    Comment


    • #3
      You're right. This was a silly question. Most people map the reads at this point and choose an alignment program that outputs the data into the proper format. I'm working to get Bowtie running on my computer, which I believe should be able to input fastq and output SAM.

      Comment


      • #4
        If you did want to convert the FASTQ to an unaligned SAM or BAM file, try this:

        Comment


        • #5
          Thanks. That's great to know.

          Comment


          • #6
            Hi,

            I've received results from my first NGS run and like you posted, I'd like to convert my fastq file to a SAM file in order to upload and retrieve data from Galaxy. I'll need to map the reads and align them to the Rice genome...but is this something I could do on my MacOSX? I'm at a loss as far as how to retrieve the sequencing results! Any help would be greatly appreciated!!

            Comment


            • #7
              Originally posted by ebatis2 View Post
              Hi,

              I've received results from my first NGS run and like you posted, I'd like to convert my fastq file to a SAM file in order to upload and retrieve data from Galaxy. I'll need to map the reads and align them to the Rice genome...but is this something I could do on my MacOSX? I'm at a loss as far as how to retrieve the sequencing results! Any help would be greatly appreciated!!
              FASTQ is the raw reads with qualities. SAM is format to describe reads and their alignment. This was hinted in the previous respond above. It seems you've got to align the reads first if you just received the raw FASTQ. Perhaps you should be asking how to align your reads. I can't help much because I am not sure what your goal is. Galaxy has a tutorial on how to align your reads and produce a SAM file. Check it out.

              Comment


              • #8
                Originally posted by husamia View Post
                FASTQ is the raw reads with qualities. SAM is format to describe reads and their alignment. This was hinted in the previous respond above. It seems you've got to align the reads first if you just received the raw FASTQ. Perhaps you should be asking how to align your reads. I can't help much because I am not sure what your goal is. Galaxy has a tutorial on how to align your reads and produce a SAM file. Check it out.
                I have the raw FASTQ reads and in order to perform de novo assembly using transAbySS, I need to feed the input in the form of bam or sam. Can you please shed some light on this.

                Comment


                • #9
                  Hi Amolkote,

                  It is rather unusual for an assembly program to accept SAM/BAM input but not FastQ. I suspect it accepts FastQ, but I don't have any experience with transAbyss. Anyways, the only tool i know that will do the job (short of an awk or perl script, which can be dangerous) was mentioned above by maubp: FastqToSam.

                  Why are you using transAbyss exactly?

                  Comment


                  • #10
                    Hi,
                    I think you cannot run transAbyss on its own. Taking a short look at the manual, (http://www.bcgsc.ca/downloads/trans-...v1.2.0.doc.pdf) I figured you probably need to run Abyss first (see "Data Preparation" in the Workflow on page 7). Plus you might have to install all the external software mentioned in "Installation, 2. External Software" (page 5). Abyss will produce contigs (.fa) and the other aligners will produce the SAM/BAM files for you, so you don't have to convert them, if I got that right. Hope it helps.

                    Cheers,
                    Stroehli
                    MSc Bioinformatics student at the Free University Berlin , Germany

                    Comment


                    • #11
                      Thanks TheRob and Stroehli !!

                      TheRob - I was using transbyss for de-novo assembly, since I don't have a concrete reference to begin with.

                      Stroehli - I have used abyss to assemble the contigs. thank you.

                      Comment


                      • #12
                        fastq to sam

                        i have a raw reads dataset in format fastq, and i want to use it to find SNPs of the transcriptome data we have. after i searched some material i found that i can do it by using Samtools and SOAPsnp softwares, am i right? but before i use them i need to convert my raw reads fastq format to SAM format, right?
                        so i installed java, samtools and picard tools on my ubuntu 12.04(why i mention these here is because i am new at linux, so any suggestion would be appreciated). and then i write this commend in the terminal :
                        java -Xmx2g -jar FastqToSam.jar FASTQ=CD_ATGTCA_L007_R1_001.fastq.gz FASTQ2=CD_ATGTCA_L007_R2_001.fastq.gz OUTPUT=outputfile.sam PREDICTED_INSERT_SIZE=null QUALITY_FORMAT=Solexa SAMPLE_NAME=file4

                        then i got this :
                        Error: Unable to access jarfile FastqToSam.jar

                        i do not know what is going on.
                        i guess many people here may done these before ,so please anyone could share your knowledge ?!

                        Comment


                        • #13
                          To access the jar, you can try giving complete path of the jar where it is installed. That should work.

                          Comment


                          • #14
                            Originally posted by kurban910 View Post
                            i have a raw reads dataset in format fastq, and i want to use it to find SNPs of the transcriptome data we have. after i searched some material i found that i can do it by using Samtools and SOAPsnp softwares, am i right? but before i use them i need to convert my raw reads fastq format to SAM format, right?
                            so i installed java, samtools and picard tools on my ubuntu 12.04(why i mention these here is because i am new at linux, so any suggestion would be appreciated). and then i write this commend in the terminal :
                            java -Xmx2g -jar FastqToSam.jar FASTQ=CD_ATGTCA_L007_R1_001.fastq.gz FASTQ2=CD_ATGTCA_L007_R2_001.fastq.gz OUTPUT=outputfile.sam PREDICTED_INSERT_SIZE=null QUALITY_FORMAT=Solexa SAMPLE_NAME=file4

                            then i got this :
                            Error: Unable to access jarfile FastqToSam.jar

                            i do not know what is going on.
                            i guess many people here may done these before ,so please anyone could share your knowledge ?!
                            SAM is the abbreviation for Sequence Alignment/Map format, which tells you that it should contain aligned/mapped reads. Though it is possible to create a kind of unmapped SAM file from fastq, this will be useless to address your question.

                            My suggestion: Make yourself familiar with read alignment via tophat (the software is here: http://ccb.jhu.edu/software/tophat/tutorial.shtml; the paper is here: http://www.nature.com/nprot/journal/....2012.016.html) and samtools in general (I suggest Dave Tang's brief wiki: http://davetang.org/wiki/tiki-index.php?page=SAMTools) and samtools mpileup in particular (http://samtools.sourceforge.net/mpileup.shtml)

                            Comment


                            • #15
                              Since you mention SNP calling, you'll want to use a tools like BWA or bowtie2 rather than tophat for alignment. Aside from that, I'm in agreement with WhatsOEver.

                              Comment

                              Latest Articles

                              Collapse

                              • seqadmin
                                Techniques and Challenges in Conservation Genomics
                                by seqadmin



                                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                                Avian Conservation
                                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                                03-08-2024, 10:41 AM
                              • seqadmin
                                The Impact of AI in Genomic Medicine
                                by seqadmin



                                Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                                02-26-2024, 02:07 PM

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by seqadmin, 03-14-2024, 06:13 AM
                              0 responses
                              34 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-08-2024, 08:03 AM
                              0 responses
                              72 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-07-2024, 08:13 AM
                              0 responses
                              81 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-06-2024, 09:51 AM
                              0 responses
                              68 views
                              0 likes
                              Last Post seqadmin  
                              Working...
                              X