Question: Extracting DE genes from annotation file
0
gravatar for clementl
2.4 years 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 • 680 views
ADD COMMENTlink modified 2.4 years ago • written 2.4 years ago by clementl0
1
gravatar for clementl
2.4 years 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 2.4 years ago by genomax87k • written 2.4 years 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: 1421 users visited in the last hour