Question: Extracting DE genes from annotation file
0
gravatar for clementl
17 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 • 464 views
ADD COMMENTlink modified 17 months ago • written 17 months ago by clementl0
1
gravatar for clementl
17 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 17 months ago by genomax71k • written 17 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: 1521 users visited in the last hour