Question: How to generate custom fasta and gtf files for RNA spike-ins
gravatar for 2405592M
6 months ago by
2405592M30 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 • 293 views
ADD COMMENTlink modified 6 months ago by michael.ante3.3k • written 6 months ago by 2405592M30
gravatar for michael.ante
6 months ago by
michael.ante3.3k 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 6 months ago by michael.ante3.3k
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: 1502 users visited in the last hour