...
Expand | ||||
---|---|---|---|---|
| ||||
Use a visualization tool, and if it’s IGV, using “load from server -> Annotations -> Genes -> Ensemble genes” will make this task much easier. |
Here is a basic command useful for parsing/sorting the gene_exp.diff
or isoform_exp.diff
files:
Code Block | ||
---|---|---|
| ||
cat isoform_exp.diff | awk '{print $10 "\t" $4}' | sort -n -r | head
|
Note that you can sort by a different column, like 12 (p-value) instead of 10; sort
can also sort on more than one column.
Now go to the GABBR1 gene in either IGV or UCSC browsers. Hint: you can type a gene name directly into the coordinate box on either browser and if it can find a match it will take you there.
...