Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Illumina HiSeq BclConverter

    Is it possible to do one lane at a time Bcl -> QSEQ.txt conversion using BclConveter? Can someone show the exact command ?

    Here is what I have been using

    ~bin/setupBclToQseq.py -i FOLDER_CONTAINING_LANE_DIRS -p POSITION_DIR -o OUT_DIR --overwrite

    followed by make in OUT_DIR

    My command stopped at lane 6 because some file was missing , but I'd like to continue with lanes 7 and 8.
    Should I be using make -j 8 so all lanes are processed simultaneously?

    Is there a '--lane=N" option for one lane at a time conversion?

  • #2
    These are the commands I usually use:

    ~/BclConverter-1.7.1/bin/setupBclToQseq.py -i /solexa/*62P5DAAXX/Data/Intensities/BaseCalls -p /solexa/*62P5DAAXX/Data/Intensities -o /solexa/*62P5DAAXX/Data/Intensities/BaseCalls --in-place --overwrite
    nohup make -j 8 > ./make.out 2> make.err

    If I got an error like that, I think I'd follow up with [email protected]. I've had quite good experiences with them.

    Comment


    • #3
      Hi,

      I asked this to illumina techsupport a while ago, and their answer was :
      The BCL converter doesn't allow you to specify lanes. I would strongly recommend not to split a RunFolder apart until after you have the qseq files. You could then run a couple of different GERALD jobs to split lanes apart. I wouldn't do it before that stage.
      So, there is not "--lane" option as you would like.

      Comment


      • #4
        Digging a little deeper, I found that "stats" file is missing for "28" and "43", which might have led to unsuccessful conversion. (See below)

        Any ideas?

        Thanks,


        Data/Intensities/BaseCalls/L006/C9.1 $ ls -1

        s_6_1.bcl
        s_6_1.stats
        s_6_21.bcl
        s_6_21.stats
        s_6_22.bcl
        s_6_22.stats
        s_6_23.bcl
        s_6_23.stats
        s_6_24.bcl
        s_6_24.stats
        s_6_25.bcl
        s_6_25.stats
        s_6_26.bcl
        s_6_26.stats
        s_6_27.bcl
        s_6_27.stats
        s_6_28.bcl
        s_6_2.bcl
        s_6_2.stats
        s_6_3.bcl
        s_6_3.stats
        s_6_41.bcl
        s_6_41.stats
        s_6_42.bcl
        s_6_42.stats
        s_6_43.bcl
        s_6_44.bcl
        s_6_44.stats
        s_6_45.bcl
        s_6_45.stats
        s_6_46.bcl
        s_6_46.stats
        s_6_47.bcl
        s_6_47.stats
        s_6_48.bcl
        s_6_48.stats
        s_6_4.bcl
        s_6_4.stats
        s_6_5.bcl
        s_6_5.stats
        s_6_61.bcl
        s_6_61.stats
        s_6_62.bcl
        s_6_62.stats
        s_6_63.bcl
        s_6_63.stats
        s_6_64.bcl
        s_6_64.stats
        s_6_65.bcl
        s_6_65.stats
        s_6_66.bcl
        s_6_66.stats
        s_6_67.bcl
        s_6_67.stats
        s_6_68.bcl
        s_6_68.stats
        s_6_6.bcl
        s_6_6.stats
        s_6_7.bcl
        s_6_7.stats
        s_6_8.bcl
        s_6_8.stats

        Comment


        • #5
          generate missing qseq files with Bustard

          rerun setupBclToQseq.py
          go into your output directory
          edit Makefile.config file
          from the line : lanes := s_1 s_2 s_3 s_4 s_5 s_6 s_7 s_8, delete the lanes for which you don't want the qseq and delete the corresponding s_?_tiles lines

          In the s_?_tiles you are interested in, delete the tiles for which you have your problem.

          for your case, you would have :

          Code:
          lanes :=  s_6 s_7 s_8
          s_6_tiles := 0001 0002 0003 0004 0005 0006 0007 0008 0021 0022 0023 0024 0025 0026 0027 0041 0042 0044 0045 0046 0047 0048 0061 0062 0063 0064 0065 0066 0067 0068
          s_7_tiles := 0001 0002 0003 0004 0005 0006 0007 0008 0021 0022 0023 0024 0025 0026 0027 0028 0041 0042 0043 0044 0045 0046 0047 0048 0061 0062 0063 0064 0065 0066 0067 0068
          s_8_tiles := 0001 0002 0003 0004 0005 0006 0007 0008 0021 0022 0023 0024 0025 0026 0027 0028 0041 0042 0043 0044 0045 0046 0047 0048 0061 0062 0063 0064 0065 0066 0067 0068

          then run the make or qmake command to generate the wanted qseq files.

          You can then get the missing qseq files for tiles s_6_28 and s_6_43 using Bustard. Please have a look à the manual as you can specify only the tiles you are interested in.

          Regards

          Gérald
          Last edited by gerald2545; 11-30-2010, 12:28 PM.

          Comment


          • #6
            Thanks, Gerald. This is very helpful.

            Comment


            • #7
              BCL Converter error

              BCL Conversion for this run stopped at lane 4 with the following error:

              Error: /storage/BclConverter-1.7.1/c++/basecalling/BclToQseqConverter.cpp: unsigned int ga::basecalling::BclReader::readNumberOfClusters(): line 122: 2010-Nov-30 13:39:36: Invalid argument: Incorrect number of clusters in /storage/101116_ILLUMINA-E84855_00027_FC_BLACKH348pM/Data/Intensities/BaseCalls/L004/C99.1/s_4_39.bcl: expected 478507: got 0

              # Why there is no cluster or incorrect cluster?
              # How would I solve this problem

              Thank you!

              Comment


              • #8
                As per a suggestion on this forum, I contacted the Illumina tech-support at [email protected] and they got back with me very quickly.

                It is worth trying for BclConversion issues.

                good luck.

                Comment


                • #9
                  To skip a lane, just create in the make directory an empty file s_N_finished.txt
                  where N stands for the lane number. In order to redo some lane, delete the corresponding "finished.txt" file and run make again.

                  Comment


                  • #10
                    Hi wanfahmi,

                    I got the same problem as you. Did you manage to solve it in the end?
                    Thanks very much!
                    JJ

                    Comment


                    • #11
                      BCL converter vs. OLB vs. upcoming CASAVA 1.8

                      We have a new HiSeq1000 and cannot get .bcl converted to fastq or qseq. Were told (by someone from illumina) that BCL converter standalone package won't work with HiSeq; OLB won't install readily (fftw issues), CASAVA 1.7 doesn't include converter and CASAVA 1.8 isn't available yet. Why is there such a disconnect between the RTA and alignment/assembly steps?

                      Comment

                      Latest Articles

                      Collapse

                      • seqadmin
                        Strategies for Sequencing Challenging Samples
                        by seqadmin


                        Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                        03-22-2024, 06:39 AM
                      • 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

                      ad_right_rmr

                      Collapse

                      News

                      Collapse

                      Topics Statistics Last Post
                      Started by seqadmin, Yesterday, 06:37 PM
                      0 responses
                      7 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, Yesterday, 06:07 PM
                      0 responses
                      7 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 03-22-2024, 10:03 AM
                      0 responses
                      49 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 03-21-2024, 07:32 AM
                      0 responses
                      66 views
                      0 likes
                      Last Post seqadmin  
                      Working...
                      X