[maker-devel] Extract CDS sequences without UTRs

[Carson Holt]

fasta_tool with the —trim_maker_utr option.

—Carson



> On Jun 27, 2021, at 7:55 PM, Sarah Baker <vireofeathers at gmail.com> wrote:
> 
> Hello,
> 
> I have annotated several genomes using the Maker2 pipeline with the goal of estimating dN/dS ratios for many genes. 
> 
> I have been using the fasta_merge script to extract the coding sequences, but I just noticed that the nucleotide sequences that it outputs (in *.all.maker.transcripts.fasta) sometimes include the 5' and 3' UTRs and they are not always in the correct reading frame. Is there a way to output CDS sequences in-frame and without UTRs (eg. so that the contents of *.all.maker.transcripts.fasta could be directly translated to the *.all.maker.proteins.fasta output by fasta_merge)?
> 
> Thank you for this great program!
> _______________________________________________
> maker-devel mailing list
> maker-devel at yandell-lab.org
> http://yandell-lab.org/mailman/listinfo/maker-devel_yandell-lab.org

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://yandell-lab.org/pipermail/maker-devel_yandell-lab.org/attachments/20210719/aed9c6db/attachment-0002.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 1376 bytes
Desc: not available
URL: <http://yandell-lab.org/pipermail/maker-devel_yandell-lab.org/attachments/20210719/aed9c6db/attachment-0002.p7s>