Thank you. It works.
• 0 views
•
link
written
by
Elizabeth •
Command for counting gene model after each round of MAKER
Hello, I am looking for a command that can tell me the number of gene models after each round of MAKER. I have tried the following:
less Maker2.all.gff |awk '$3=="mRNA"' |grep "mRNA-1" |awk '{print $5-$4}' |summary.sh
cat <roundN.full.gff> | awk '{ if ($3 == "gene") print $0 }' | awk '{ sum += ($5 - $4) } END { print NR, sum / NR }'
But, my terminal seems to hang, every time I run the command. Has anyone encountered a similar problem ?
• 2,061 views
•
link
updated
by
ringourquia •
•
written
by
Elizabeth •
•
written
by
Elizabeth •
1 answer
The "fathom -gene-stats" command used in the path to generate the SNAP learning-file displays several stats of your gff including the number of genes. Otherwise, I use:
cut -f3 your_merged_gff.gff | grep "gene" | wc -l
• 0 views
•
link
written
by
ringourquia •
Log in to answer this question.
More posts like this
-
Average length of mRNA on NCBI annotation seems too high
written by yh1126 •
-
How to filter and keep longest isoform per gene in gff3 file ?
written by Sony Nguyen •
-
Using samtools to determine Rh blood type
written by esiv •
-
Multiple commands in a for loop
written by lamia_203
-
Getting error with awk when using parallel processing in bedtools
written by nastaran.esfahani •
-
SNAP error in MAKER exon:out_of_bounds
written by preethibagopi5 •
-
Forum: Awk in Bioinformatics
written by Shicheng Guo
-
AWK: if value falls within a range, print sum of values
written by heso •
-
Bash script error on reading coordinates in a file
written by Ali HEBRA •
-
Linux command line use to rename a gff annotation file "Name" attribute
written by rob234king
What does
summary.sh?