Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Conversion of qseq.txt format to fastq

    Hello Folks,

    Could any one please let me know,if their is a way to convert the qseq.txt format files of illumina to fastq format. I finished my run and I would like to align my reads to ref seq using MAQ. MAQ requires the read data to be in fastq. Currently I've my data in qseq.txt format.

    Any suggestions would be greatly appreciated.

    Thank you

    Rakesh

  • #2
    Using the forums search tool would have found these threads:

    Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

    Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

    Comment


    • #3
      You can try the following one liner in Bash: go to where those _qseq.txt files locate, and if those qseq.txt are in form of s_N_1_M_qseq.txt where N is lane number (1-8), M is tile number (1-120), then

      Code:
      $ for ((x=1;x<=8;x+=1)); do cat s_"$x"_1_*_qseq.txt | awk -F '\t' '{gsub(/\./,"N", $9); if ($11 > 0) printf("@%s_%04d:%s:%s:%s:%s#%s/%s\n%s\n+%s_%04d:%s:%s:%s:%s#%s/%s\n%s\n",$1,$2,$3,$4,$5,$6,$7,$8,$9,$1,$2,$3,$4,$5,$6,$7,$8,$10)}' > s_"$x"_sequence.fastq; done
      The resulting s_N_sequence.fastq (N is 1->8) are the 8 fastq files you will need.

      Comment


      • #4
        Hi Dukeven,

        On the Unix machine, I went to the directory where my qse.txt file were and typed in the command Bash and pasted your code

        $ for ((x=1;x<=8;x+=1)); do cat s_"$x"_1_*_qseq.txt | awk -F '\t' '{gsub(/\./,"N", $9); if ($11 > 0) printf("@%s_%04d:%s:%s:%s:%s#%s/%s\n%s\n+%s_%04d:%s:%s:%s:%s#%s/%s\n%s\n",$1,$2,$3,$4,$5,$6,$7,$8,$9,$1,$2,$3,$4,$5,$6,$7,$8,$10)}' > s_"$x"_sequence.fastq; done

        It shows:

        syntax error near unexpected token `('


        I'm new so....... not sure why it doesn't work

        Thanks

        Comment


        • #5
          It works just fine here. Dont know why you got that error. Just to be sure that the command should be in one line only in your terminal.

          Comment


          • #6
            rakeshponnala you had this $ 5 instead of $5

            Comment


            • #7
              I ran the script, and msg is "syntax error near unexpected token `('"

              Comment


              • #8
                Maybe I am missing something, but why don't you use the normal Illumina pipeline with Casava?

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Choosing Between NGS and qPCR
                  by seqadmin



                  Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
                  10-18-2024, 07:11 AM
                • seqadmin
                  Non-Coding RNA Research and Technologies
                  by seqadmin




                  Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

                  Nobel Prize for MicroRNA Discovery
                  This week,...
                  10-07-2024, 08:07 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 11-01-2024, 06:09 AM
                0 responses
                13 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 10-30-2024, 05:31 AM
                0 responses
                16 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 10-24-2024, 06:58 AM
                0 responses
                24 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 10-23-2024, 08:43 AM
                0 responses
                52 views
                0 likes
                Last Post seqadmin  
                Working...
                X