Take the 2-minute tour ×
Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. It's 100% free, no registration required.

I want to align WGS data from several mammals to one reference such as the human genome sequence. Since in most cases exons should be conserved and spliced in the same way and introns should vary, I want to use a tool that is aware of the annotation of exon and intron.

That is to say that matches in an exon region should have be given a higher weighting than matches in intronic regions.

Also, when a read pair maps to 2 exons, the insert size should be more flexible as the intron length can be vary among different mammals. Since the library insert is less than 2kb for short libraries, the more common situation should be one read on exon and another on intron.

This is different from RNA-seq as I am using WGS data of a whole genome, so both exons and introns are in the reads. And I want to align them to a human reference sequence.

share|improve this question

1 Answer 1

If you are not trying to assemble but just to align each read to the genome, you can use exonerate. On a Unix/Linux platform, once you have installed it run something like:

exonerate -m genome2genome WGS.fasta genome.fasta > out.txt 

From the exonerate manual:

          genome2genome
                 This  model  is  similar  to  the  cod‐
                 ing2coding  model,  except  introns are
                 modelled on both sequences.  (not work‐
                 ing well yet)

What I would recommend though, is to align against a reference cDNA dataset, not the whole genome. In that case, you should use this instead:

exonerate -m cdna2genome genome_cdna.fasta WGS.fasta > out.txt 

From the exonerate manual:

          cdna2genome
                 This   combines   properties   of   the
                 est2genome and coding2genome models, to
                 allow modeling of an whole cDNA where a
                 central coding region can be flanked by
                 non-coding UTRs.  When  the  CDS  start
                 and  end  is  known it may be specified
                 using  the  --annotation  option   (see
                 below)  to permit only the correct cod‐
                 ing region to appear in the alignemnt.
share|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.