Question: How to generate custom fasta and gtf files for RNA spike-ins
gravatar for 2405592M
2.0 years ago by
2405592M90 wrote:

Hi guys, I've been doing some RNA-seq on mice and I used 3x E. coli tRNA spike-ins. I'm assuming for analysis, I'll have to generate a fasta and GTF file for these spike-ins? If so, how do I do this? Also, I'm assuming this will have to be merged with the mouse fasta and GTF files. How do I do this as well? Once these steps are completed, can I move on to alignment? The spike-ins were custom made oligos and so I only have the 3 sequences. I don't know any gene loci/coordinates etc.

rna-seq spike-in • 834 views
ADD COMMENTlink modified 2.0 years ago by michael.ante3.6k • written 2.0 years ago by 2405592M90
gravatar for michael.ante
2.0 years ago by
michael.ante3.6k wrote:

Assuming the tRNAs are unspliced sequences, you can add them like the ERCC92 spike in annotation: Each spike in has its own contig and each gene annotation goes from start to the end.

Something like:


and the GTF:

mySpikeIn1   spikein  exon 1  9999 0.0  +  . gene_id "mySpikeIn1"; transcript_id "mySpikeIn1.1";

Please refer to the GTF file format documentation regarding spaces/tabs and necessary fields.



ADD COMMENTlink written 2.0 years ago by michael.ante3.6k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2121 users visited in the last hour