Hello, a quick question for anybody familiar with bowtie. In my reference genome files, which have names like SVID_XXX, the first line of text, before the actual sequence data is a line like >Chromosome, or >ContigXXXX, .... In the bowtie output data, Chromosome or ContigXXXX shows up as the file name it mapped to instead of SVID_XXX, the actual file name. My question is this, is bowtie taking the > symbol in the first line of text as a command to redirect to this file name or whatever comes after it in that line? I know that > redirects output to a file, but why does this happen in bowtie since the file already has a file name, and > only shows up in the first line of a text file? I hope this makes sense, any thoughts?
-
-
Your reference genome files are FASTA files.
The FASTA file format is simple.
There is one line, starting with the symbol ">" which describes the sequence.
The first word after ">" is the name of the sequence.
Following the sequence description is the sequence itself.
Bowtie and all programs dealing with raw genomic data are written to be able to parse FASTA files.
The FASTA format has nothing to do with the Unix redirection operator, and dates back to 1985.Last edited by blancha; 11-13-2015, 07:50 PM. -
Thanks, that helps a lot.Comment
-
Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.
We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing...
Comment