This is a test version of Biostars. For the public version, visit https://www.biostars.org.
remove scaffold and other unplaced sequence before mapping ?

Hi,
I downloaded reference genomes from Ensembl (fasta format). But there are lots of sequences with name "dna:scaffold": https://github.com/CTLife/TEMP/tree/master/RefGenomes

Such as Mouse_GRCm38 (mm10), except chromosome 1-19, Mt, X and Y; others should be removed before mapping ? https://github.com/CTLife/TEMP/blob/master/RefGenomes/Mouse_GRCm38.p4.txt

Such as Human_GRCh38.p5 (hg38), https://github.com/CTLife/TEMP/blob/master/RefGenomes/Human_GRCh38.p5.txt, there are 516 sequences. In addition to chromosome 1-22, Mt, X and Y; others (such as CHR_HG2241_PATCH and KI270728.1) should be removed before mapping ?

rna-seq chip-seq genome sequencing next-gen
• 5,503 views • link
updated by abascalfederico abascalfederico • written by yongxpeng • yongxpeng

1 answer

The latest release of the human genome (don't know about mice) contains alternative contigs. You will need an alternative-contig aware algorithm like BWA: https://github.com/lh3/bwa/blob/master/README-alt.md

If you are not using one of this kind of algorithms it would be better to remove the alternative contigs. That's because a read may map to multiple alternative contigs and be (incorrectly) considered a non-uniquely mapped read.

HTH

• 1 views • link
written by abascalfederico abascalfederico

If I understood well that link is about repetitive sequences, not about alternative contigs

For RNA-seq... it depends. For example, if you want to analyse HLA genes, which are highly diverse, you would need the alternative contigs. I guess most people just ignore alternative contigs because of the increase in complexity.

• 1 views • link
written by abascalfederico abascalfederico

Log in to answer this question.

More posts like this