Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • maria_mari
    replied
    Hi,
    size limit in last versions of bwa (use bwa 0.6 )

    Leave a comment:


  • Jon_Keats
    replied
    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)

    Leave a comment:


  • louis7781x
    replied
    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!

    Leave a comment:


  • Jon_Keats
    replied
    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.

    Leave a comment:


  • maubp
    replied
    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).

    Leave a comment:


  • louis7781x
    replied
    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.

    Leave a comment:


  • maubp
    replied
    Did you decompress it properly? e.g. try:

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

    Leave a comment:


  • louis7781x
    replied
    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!!!

    Leave a comment:


  • maubp
    replied
    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?

    Leave a comment:


  • louis7781x
    started a topic A question about BWA index

    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!!!

Latest Articles

Collapse

  • seqadmin
    Quality Control Essentials for Next-Generation Sequencing Workflows
    by seqadmin




    Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

    Nucleic Acid Quality Control
    Preparing for NGS starts with isolating the...
    Yesterday, 01:58 PM
  • seqadmin
    An Introduction to the Technologies Transforming Precision Medicine
    by seqadmin


    In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
    01-27-2025, 07:46 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 02-07-2025, 09:30 AM
0 responses
19 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-05-2025, 10:34 AM
0 responses
37 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-03-2025, 09:07 AM
0 responses
37 views
0 likes
Last Post seqadmin  
Started by seqadmin, 01-31-2025, 08:31 AM
0 responses
41 views
0 likes
Last Post seqadmin  
Working...
X