| |
Summary
We will visualize the output produced by SIM on the genomic MHC syntenic region between human and mouse. We will use a filter to convert its output into GFF, the C program ali2gff is included in the gff2aplot distribution.
Once we have generated the whole genomic region figure, we will zoom into a single gene annotation. This will illustrate how you can expand a given area from your pictures, to highlight single features on multigenic annotation plots for instance.
NOTE.- For the sake of clarity, we are going to use long names for the comand-line switches. See the command-line help if you prefer short names for those cases in which a short name is available.
Bitmaps for the examples were generated as PNGs (Portable Network Graphics). If your browser is not ready for such format yet, you can visualize the PDF or PS versions by clicking on the links below each snapshot. Links to customization files, log files, GFF input files and output PostScript figures, are also available on each command-line shown.
Contents
|
Parsing SIM Output
|
|
|
Raw Output
|
|
[PNG] [PS] [PDF]
| |
gff2aplot.pl \
--verbose \
--percent-box-label "Identity%" \
--show-percent-box \
--title "Human/Mouse MHCII Region" \
--subtitle "" \
--feature-var "seqbounds::alignment_color=grey" \
--feature-var "alignment::alignment_color=red" \
--feature-var "fragment::alignment_color=green" \
--strand-var "..::strand_layer=0" \
--strand-var "++::strand_layer=2" \
-- hs.gff \
mm.gff \
hs-mm.sim.gff \
> hs-mm.sim2.ps \
2> hs-mm.sim2.log
| |
|
Finishing Plot for MHC Region
|
|
|
Zooming Into the LMP2 Gene Region
|
|
[PNG] [PS] [PDF]
| |
gff2aplot.pl \
--verbose \
--percent-box-label "Identity%" \
--show-percent-box \
--title "Human/Mouse MHCII Region" \
--subtitle "Zooming into the LMP2 gene region" \
--zoom-area \
--end-x-sequence 6500 \
--end-y-sequence 6500 \
--x-label "Human X87344" \
--y-label "Mouse AF100956+AF027865" \
--custom-filename mhc2.rc \
--custom-filename mhc2labels.rc \
-- hs.gff \
mm.gff \
hs-mm.sim.gff \
> hs-mm.sim.zbox.ps \
2> hs-mm.sim.zbox.log
| |
|
|
