I am trying to analyze RNA-seq data in DESeq in Galaxy and wonder if anyone has a detailed instructions or work flow how DEseq can be used after alignment in galaxy. Any pointer should be helpful. Since it is galaxy question I have also posted similar question on galxay but though this area may have better coverage. Thanks
Either that species lacks UTRs or it's unknown where they are. So try using
CDS
instead ofexon
. Note that I'm not sure how well htseq-count handles GFF files, so if that doesn't work then just post a comment here and I'll make a proper GTF file for you.I am using HTseq on galaxy and it seems it is not accepting option CDS in feature type box. It will be a great help if you could make a proper GTF file for me. Thanks in advance :)
GFF file for brassica looks like this....
Seqname Source Feature Start End Score Strand Frame Group Scaffold02034 glean gene 5106 6434 . + 0 ID=BjuA000715; Scaffold02034 glean mRNA 5106 6434 . + 0 ID=BjuA000715;Parent=BjuA000715; Scaffold02034 glean CDS 5106 5174 0.99 + 0 ID=BjuA000715.cds;Parent=BjuA000715; Scaffold02034 glean CDS 5275 5425 0.97 + 0 ID=BjuA000715.cds;Parent=BjuA000715; Scaffold02034 glean CDS 6307 6434 0.61 + 2 ID=BjuA000715.cds;Parent=BjuA000715; Scaffold00650 glean gene 12497 13500 459.56 + . ID=BjuB006329; Scaffold00650 glean mRNA 12497 13500 459.56 + . ID=BjuB006329;Parent=BjuB006329;
Here you go. Make sure you explicitly label the file type as "gtf" when you upload it to Galaxy.
Thank you very much Sir ...