This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to use GSNAP/GMAP for mm10 genome

Hello,

I am not sure how to download the mm10 genome for genome guided alignment on GSNAP. In a separate project (for an Aiptasia alignment), I used an .fna file for building the gmap database and a .gff file for the intron store building, and I used Trinity subsequently (these files were provided for me and I didn't download them so I don't know where to get them). I wish to follow this protocol for my mm10 alignment for my mouse sequences if possible.

I appreciate your help!

assembly alignment genome rna-seq
• 2,627 views • link
updated by GenoMax GenoMax • written by Moneeb Bajwa • Moneeb Bajwa

1 answer

This is the GENCODE official data page for all things related to mouse genome. You will find the files you need there (fasta genome sequence, GTF annotation).

• 0 views • link
written by GenoMax GenoMax

Thank you. Does it matter which combination I choose? Could I choose the first in the gff list with the first in the fasta file list? I also see that it is .fa as opposed to .fna; does that matter?

• 0 views • link
written by Moneeb Bajwa • Moneeb Bajwa

.fa should be equivalent to .fna. Please read the descriptions for the files carefully before selecting. If you are not interested in patches and haplotypes then choose the files labelled as PRI (Primary) only.

• 0 views • link
written by GenoMax GenoMax

Also what if I did not unzip the files when supplying them to GSNAP and the other tools; would it work?

• 0 views • link
written by Moneeb Bajwa • Moneeb Bajwa

It would depend on the tool. Some tools can accept gzipped files others you will need to unzip them.

• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this