GMAP or Bowtie index from FASTA of genes
0
0
Entering edit mode
5.9 years ago
gtasource ▴ 60

Hi all,

I'm looking to use either Bowtie2 or GMAP/GSNAP to try to map sequences to gene names. The FASTA file that I want to index contains all of the genes from a species. Problem is, I keep getting errors when I try to index this FASTA file, as GMAP wants a chromosomal assembly I believe. The problem with using the chromosome assembly as an index, is it won't tell me what that gene is (hence why I create my own FASTA with gene names.) Does anybody have any suggestions for this?

Cheers!

GMAP bowtie • 2.0k views
ADD COMMENT
1
Entering edit mode

You should align to the genome, and GMAP is perfectly capable of doing so. It is however unclear which sequences you are trying to map. If you have all genes from a species you most likely could also get a GTF/GFF annotation file from this species, which would tell you which gene is where on the genomic DNA.

For now it is unclear which problem you are trying to solve.

ADD REPLY
1
Entering edit mode

I keep getting errors

And what are the errors?

ADD REPLY

Login before adding your answer.

Traffic: 2819 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6