File tree Expand file tree Collapse file tree
data/genomics/homo_sapiens Expand file tree Collapse file tree Original file line number Diff line number Diff line change @@ -228,6 +228,21 @@ cat << EOF > GenomeSize.xml
228228EOF
229229```
230230
231+ Similarly, the GenomeSize.xml index for the entirety of chr22 (` homo_sapiens/genome/chr22/sequence/GenomeSize.xml ` ) was created accordingly:
232+
233+ ``` bash
234+ stripped_seq=$( mktemp)
235+ tail -n +2 hg38.chr22.fasta | tr -d ' \n' > $stripped_seq
236+ length=$( wc -c $stripped_seq | cut -f 1 -d ' ' )
237+ md5=$( md5sum $stripped_seq | cut -f 1 -d ' ' )
238+
239+ cat << EOF > GenomeSize.xml
240+ <sequenceSizes genomeName="Homo Sapiens (NCBI GRCh38)">
241+ <chromosome fileName="genome.fa" contigName="chr22" totalBases="${length} " build="GRCh38" isCircular="false" md5="${md5} " ploidy="2" species="Homo_sapiens" knownBases="${length} " type="Autosome" />
242+ </sequenceSizes>
243+ EOF
244+ ```
245+
231246### Genome map
232247
233248There is multiple type of genetic map depending on the softwares.
Original file line number Diff line number Diff line change 1+ <sequenceSizes genomeName =" Homo Sapiens (NCBI GRCh38)" >
2+ <chromosome fileName =" genome.fa" contigName =" chr22" totalBases =" 50818468" build =" GRCh38" isCircular =" false" md5 =" ac37ec46683600f808cdd41eac1d55cd" ploidy =" 2" species =" Homo_sapiens" knownBases =" 50818468" type =" Autosome" />
3+ </sequenceSizes >
4+
You can’t perform that action at this time.
0 commit comments