Skip to content
Snippets Groups Projects
Commit 53df9825 authored by tcarver's avatar tcarver
Browse files

add indexed GFF support

parent f93db0b1
No related branches found
No related tags found
No related merge requests found
...@@ -44,11 +44,34 @@ This function only reads the feature section of the input file - the sequence ...@@ -44,11 +44,34 @@ This function only reads the feature section of the input file - the sequence
</LISTITEM> </LISTITEM>
<LISTITEM> <LISTITEM>
<PARA> <PARA>
FASTA files and indexed FASTA files can be read. The files are indexed FASTA files
using <ULINK URL="http://samtools.sourceforge.net/">SAMtools</ULINK> (e.g. </PARA>
<COMMAND>samtools faidx ref.fasta</COMMAND>). A drop down menu of the contigs or chromosomes </LISTITEM>
is provided in the Entry toolbar to select the sequence. Using indexed FASTA files <LISTITEM>
improves the memory management and enables large genomes to be viewed. <PARA>
Indexed FASTA files can be read in. The files are indexed
using <ULINK URL="http://samtools.sourceforge.net/">SAMtools</ULINK>:
</PARA>
<SYNOPSIS>
samtools faidx ref.fasta
</SYNOPSIS>
</LISTITEM>
<LISTITEM>
<PARA>
The indexed FASTA can be used with an indexed GFF to overlay annotation on the sequence.
To index the GFF first sort and bgzip the file and then use tabix with "-p gff" option (see the
<ULINK URL="http://samtools.sourceforge.net/tabix.shtml">tabix manual</ULINK>):
</PARA>
<SYNOPSIS>
(grep ^"#" in.gff; grep -v ^"#" in.gff | sort -k1,1 -k4,4n) | bgzip > sorted.gff.gz;
tabix -p gff sorted.gff.gz
</SYNOPSIS>
<PARA>
A drop down menu of the contigs or chromosomes sequences is provided in the Entry toolbar
to select the sequence. Using indexed FASTA and indexed GFF files improves the memory management
and enables large genomes to be viewed. Note that as it is indexed the sequence and annotation are
read-only and cannot be edited. When there are many contigs to select from it can be easier
to display the one of interest by typing the name into the drop down list.
</PARA> </PARA>
</LISTITEM> </LISTITEM>
<LISTITEM> <LISTITEM>
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment