Hi Everyone,
Sorry if this is a dup, but I can't seem to find a satisfactory answer on this site or others.
I'm wondering what, if any, pre-processing I should perform on a reference genome fasta/gff file prior to mapping using BWA or Bowtie. For example, if I wanted to map something to the orangutan genome, should I remove entries that are labeled as "unplaced/unlocalized genomic scaffold" from the gff and fasta files--i.e. only map to canonical chromosomes.
I notice that even these scaffolds have "BestRefSeq" categories in the gff file for genes, indicating that they still have useful information on them.
The reason I ask is because I was told by someone who no doubt knows much more than me about this stuff that I SHOULD remove these chromosomes. I'm wondering, however if this person is wrong.
Thanks!
1 answer
That is your choice.
If the sequence is known to belong to the genome (but is unplaced at the moment) then it should remain. You can ignore reads that align there if you are not interested in that region. Omitting the region from the reference may force an aligner to align those reads elsewhere, which you would probably do not want.
Log in to answer this question.
-
how to split a Bam by chromosome and unplaced scaffolds.
written by Maxine •
-
Is it a problem that the reference genome is not at the chromosome level?
written by beausoleilmo
-
What's the right way to submit a genome assembly to NCBI?
written by liorglic
-
samtools: splitting a bam file putting all scaffolds together
written by Fedster •
-
gffread query for gff to gtf conversion
written by prasundutta87
-
Forum: Filtering out chromosomes from reference genome
written by Bastien Hervé
-
Using satsuma output when creating manhattan plots
written by ash.t.cook •
-
scaffold in reference genome
written by ft.mastouri •
-
GFF and Genome Fasta file
written by aj123
-
Discrepancy In The Number Of Scaffolds Of A Genome Assembly
written by nobrainw •
