Question: Extracting DE genes from annotation file
0
gravatar for clementl
12 months ago by
clementl0
clementl0 wrote:

I'm working with differentially expressed genes from different comparisions. I ran HISAT2, StringTie and DESeq2 and now I have txt files of all my DE genes from the 9 comparisons. I'm trying to extract out the locations and names of the DE genes from my merged annotation file (I have GTF and FASTA format of this file) and I can't draw them out.

TXT file ex of DE genes: MSTRG.1 MSTRG.1356 MSTRG.1345

I've used grep, samtools xargs faidx, and I still cannot figure it out. Any advice? Thank you!

de genes next-gen • 350 views
ADD COMMENTlink modified 12 months ago • written 12 months ago by clementl0
1
gravatar for clementl
12 months ago by
clementl0
clementl0 wrote:

save as .sh

 #!/bin/bash
while read p; do
grep -w -A1 "$p" mergedfasta.fa | head -n 2 >> comparison.de.fa
done < degenes.txt
ADD COMMENTlink modified 12 months ago by genomax65k • written 12 months ago by clementl0
Please log in to add an answer.

Help
Access

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