Hi friends,
Got a small doubt about a genome file. I don't know if my chromosomes are already ordered in my fasta file
Actually I think it is because it's not 1 to 17.
Just want to confirm with you before to extract the region I wanted after a blast. My indexed genome order:
Advance thanks
1 answer
Ordering is whatever you need it to be. Some common sorting schemes (sort-bed, for instance) use "lexicographic" ordering of chromosome names, so as to take advantage of that structure. In that scheme, the chromosome label "chr11" is positioned before "chr2", for example. Other schemes can use numeric sorting, where "chr2" would come before "chr11". It just depends on where your data come from and what you need to do.
•
written
by
Alex Reynolds
Log in to answer this question.
-
Interpreting HOMER getDifferentialPeaks.pl output columns
written by Ronin •
-
BAM 2 FASTQ using Picard and Samtools
written by kcarey •
-
How to use samtools faidx
written by florian.grosmaire •
-
Extract CSV and SAM
written by mustafa_aljadi •
-
How to extract certain CDS from a GenBank file in linux terminal
written by Anisur Rahman •
-
featureCounts log: "WARNING contig is not found in the provided genome file!"
written by Andrés Ribone •
-
Read counts from bam file without available annotation file
written by E •
-
abyss-pe file prob double /1
written by Siva •
-
correct fastq file with SED
written by Siva •
-
How to generate a .sam/.bam file of a particular chromosome or even a particular region?
written by Chenglin
I would like to help you but I don't really understand your question.
Do you want to sort the chromosomes in your fasta file? The file you have attached it is sorted alphabetically, not by chromosome. What are you trying to do exactly?
I search a way to know if chromosomes of my genome are in this order.
Actually I got a genome file, and wanted to locate some genes and retrieve the part of the genome, but before that I should know if my genome has already ordered chromosomes.
I would suggest you to do that using
bedtools getfasta. Read here for more info: http://bedtools.readthedocs.org/en/latest/content/tools/getfasta.html.Nope airan I can't make a bedtools because haven't got bed file, and i don't want to retrieve sequences.
The file I give before is already from samtools faidx so no use to take that to make a getfasta.
I just want to know if the genome I have is already a contiguous one so with a good order of his chromosomes.
Your .fasta file doesn't seen to be "karyotypically" ordered.
You should read this question -> Karyotypically Ordered Hg19
A simple way to check the order in your fasta file would be:
Yeah when I did a grep I got the same order than in this file:
The question is: so in my genome the order of the chromosomes is what it's given by the grep result?
It could not be an other order ? In the real genome we'll find it in that order, am I right?
grep will give you the same order as chromosomes appear in your fasta file. If you want to change the order of chromosomes and create a new file, you should be able to do it using Picard. Search biostars for relevant posts. You can use keywords such as "change order chromosomes fasta".
OK thanks guys.
I think chromosomes are conserved in the genome, I'll continue with that.
P.S.: I don't want to change the order, just wanted to confirm the real order in the real genome.