Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • A question about BWA index

    I have use bwa index Homo_sapiens.GRCh37.59.dna.toplevel.fa (download from Ensembl) ,the file is 4.9G

    and when I index it,


    [bwa_index] Pack FASTA... 121.12 sec
    [bwa_index] Reverse the packed sequence... 24.93 sec
    [bwa_index] Construct BWT for the packed sequence...
    TextLengthFromBytePacked(): text length > 2^32!
    louis@bio001:~/program/bwa-0.5.8c$
    louis@bio001:~/program/bwa-0.5.8c$ ./bwa bwasw genomes/Homo_sapiens.GRCh37.59.dna.toplevel.fa reads/1_2_RT.fastq > hgoutput.sam
    [bwt_restore_bwt] fail to open file 'genomes/Homo_sapiens.GRCh37.59.dna.toplevel.fa.bwt'. Abort!
    Aborted

    I dont't know what happened ,if someone know,please let me know


    Best Regard!!!

  • #2
    I think something is wrong with your FASTA file - the index failed, apparently your sequence is too long to index (2 to the power of 32 bases is very big, 4.2 billion!).

    What URL did you download the FASTA file from?

    Comment


    • #3
      Originally posted by maubp View Post
      I think something is wrong with your FASTA file - the index failed, apparently your sequence is too long to index (2 to the power of 32 bases is very big, 4.2 billion!).

      What URL did you download the FASTA file from?

      hi,it is my download fils 's url ftp://ftp.ensembl.org/pub/current/fa...toplevel.fa.gz

      The file's sorce is from ensembl.

      Would you help me find the error thanks!!!

      Comment


      • #4
        Did you decompress it properly? e.g. try:

        head Homo_sapiens.GRCh37.59.dna.toplevel.fa

        Comment


        • #5
          Originally posted by maubp View Post
          Did you decompress it properly? e.g. try:

          head Homo_sapiens.GRCh37.59.dna.toplevel.fa


          I use command "gunzip Filename.gz" to decompress this file.

          sorry I don;t understand "head" What is this command?

          Thanks!
          Last edited by louis7781x; 11-01-2010, 07:46 AM.

          Comment


          • #6
            head is a Unix command to see the start of a text file (short for header I think), tail shows you the end of a text file (head and tail being the opposite ends of an animal).

            Comment


            • #7
              To he best of my understanding you can't use the top level files as the size exceeds the maximum supported by the BWT used in BWA. This is because the top level files include entire duplicate chromosomes for the different haplotypes. Most people are using the 1000 genomes version of GRCh37.

              Comment


              • #8
                Originally posted by Jon_Keats View Post
                To he best of my understanding you can't use the top level files as the size exceeds the maximum supported by the BWT used in BWA. This is because the top level files include entire duplicate chromosomes for the different haplotypes. Most people are using the 1000 genomes version of GRCh37.
                hi Jon ,my research is to find gene fusions in brain tumor's cDNA library generated from 454.
                I read many papers,and they usually use 454 data align against to hg19 and refseq of cDNA ,and just extract "non-mapping reads". then,using non-mapping reads to find the read where can align across to two exon.

                I don't understand .In my research,Is 1000 genomes useful ?

                Thanks!

                Comment


                • #9
                  Hi Louis,

                  The 1000 genomes version of human genome build GRCh37 should be as useful as the top level file at ensembl, maybe more so as you can the get BWA running. I'm assuming your analysis strategy is to map all reads to human genome, take all those that don't map, and map against human transcriptome, then take those that still do not map and blast against genome to look for novel hybrid junctions? The only difference is they curated this version to get rid of the redundant duplication that is not necessary and is likely to cause problems in your analysis. If you want a bit more explaination see my thread (http://seqanswers.com/forums/showthread.php?t=4589)

                  Comment


                  • #10
                    Hi,
                    size limit in last versions of bwa (use bwa 0.6 )

                    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
                    12 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, Yesterday, 06:07 PM
                    0 responses
                    10 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-22-2024, 10:03 AM
                    0 responses
                    51 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-21-2024, 07:32 AM
                    0 responses
                    68 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X