Hi,
size limit in last versions of bwa (use bwa 0.6 )
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
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:
-
Originally posted by Jon_Keats View PostTo 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.
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:
-
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:
-
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:
-
Originally posted by maubp View PostDid 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:
-
Did you decompress it properly? e.g. try:
head Homo_sapiens.GRCh37.59.dna.toplevel.fa
Leave a comment:
-
Originally posted by maubp View PostI 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:
-
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:
-
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!!!Tags: None
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 08:47 AM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
54 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
Leave a comment: