# working folder
	/var/www/jbrowse/Martijn_Sornchai

	

# format refseq 
	copy all fna to /var/www/jbrowse/Tonia_Sporulation/docs/fna
	copy all gff to /var/www/jbrowse/Tonia_Sporulation/docs/gff

# in this option only BS168 is added because all strains are mapped on BS168
bin/prepare-refseqs.pl --fasta data/dataset/virgibacillus_with_plasmids.fna

# add annotation
# if not exists:  mkdir /data/tracks
bin/flatfile-to-json.pl --gff data/dataset/virgibacillus_with_plasmids.gff --trackLabel Genbank_Annotation --trackType CanvasFeatures



# make feature searchable
# bin/prepare-refseqs.pl -gff docs/gff/xxxxx.gff
bin/prepare-refseqs.pl -gff data/dataset/virgibacillus_with_plasmids.gff


# active the names
	bin/generate-names.pl

# make link to the /data folder
	# First data set of 4 genomes T01..T56
	ln -s /home/anne/PROJECTS/tonia/2014-09_bsubt168 dataset01
	# second data set
	ln -s /home/anne/PROJECTS/yi_tonia/2015_08 dataset02



# make data/tracks.conf file
# usefull tool is to convert a table direct to a tracks.conf file.
# goto /data and run:
/usr/molgentools/tools/jbrowse_table_2_tracks.pl -i tracks_conf.table > tracks.conf

# Change default spacing between tracks open jbrowse.conf and add the line "view.trackPadding = 7"

# If BAM is empty, check the genome header name
/usr/molgentools/tools/find_and_replace_single.pl -i MG1363_A.sorted.bam -f 'gi|125622882|ref|NC_009004.1|' -r 'NC_009004.1' -o MG1363_A.sorted.new.bam