...
LoadGenomeIntervals.sh
- TASSEL requires 50G for this assembly, add Xmx=50G to configSQLite.txt
LoadAssemblyAnchors.sh - do not use, it is being rewritten
CreateHaplotypes.sh
- change gatk-launch to gatk
- fix arguments
- index documents - https://gatkforums.broadinstitute.org/gatk/discussion/1601/how-can-i-prepare-a-fasta-file-to-use-as-reference
java -jar picard.jar CreateSequenceDictionary R=/data2/sra/161010_Chinese_Spring_v1.0_pseudomolecules.fasta
java -jar picard.jar CreateSequenceDictionary R=tmpFileDir/data/reference/161010_Chinese_Spring_v1.0_pseudomolecules.fasta
samtools faidx 161010_Chinese_Spring_v1.0_pseudomolecules.fasta
- intervals.bed - change contigs from 1A to chr1A and so on
...