b3833e993f98922e22e1278e632716e729ec54ad hiram Fri Jul 19 15:24:34 2019 -0700 adding gapOverlap and tandemDups tracks to the asmHub build refs #23734 diff --git src/hg/utils/automation/asmHubTrackDb.sh src/hg/utils/automation/asmHubTrackDb.sh index c3fd82a..2a21dce 100755 --- src/hg/utils/automation/asmHubTrackDb.sh +++ src/hg/utils/automation/asmHubTrackDb.sh @@ -75,30 +75,78 @@ visibility full autoScale Off maxHeightPixels 128:36:16 graphTypeDefault Bar gridDefault OFF windowingFunction Mean color 0,0,0 altColor 128,128,128 viewLimits 30:70 type bigWig 0 100 bigDataUrl bbi/%s.gc5Base.bw html html/%s.gc5Base\n\n" "${asmId}" "${asmId}" $scriptDir/asmHubGc5Percent.pl $asmId $buildDir/html/$asmId.names.tab $buildDir > $buildDir/html/$asmId.gc5Base.html fi +# see if there are gapOverlap or tandemDup bb files +export gapOverlapCount=0 +export tanDupCount=0 +if [ -s $buildDir/trackData/gapOverlap/${asmId}.gapOverlap.bb ]; then + gapOverlapCount=`zcat $buildDir/trackData/gapOverlap/${asmId}.gapOverlap.bed.gz | wc -l` + rm -f $buildDir/bbi/${asmId}.gapOverlap.bb + ln -s $buildDir/trackData/gapOverlap/${asmId}.gapOverlap.bb $buildDir/bbi/${asmId}.gapOverlap.bb +fi +if [ -s $buildDir/trackData/tandemDups/${asmId}.tandemDups.bb ]; then + tanDupCount=`zcat $buildDir/trackData/tandemDups/${asmId}.tandemDups.bed.gz | wc -l` + rm -f $buildDir/bbi/${asmId}.tandemDups.bb + ln -s $buildDir/trackData/tandemDups/${asmId}.tandemDups.bb $buildDir/bbi/${asmId}.tandemDups.bb +fi + +if [ "${gapOverlapCount}" -gt 0 -o "${tanDupCount}" -gt 0 ]; then + + printf "track tanDups +shortLabel Tandem Dups +longLabel Paired identical sequences +compositeTrack on +visibility hide +type bigBed 12 +group map +html html/%s.tanDups\n\n" "${asmId}" + + if [ "${gapOverlapCount}" -gt 0 ]; then + printf " track gapOverlap + parent tanDups on + shortLabel Gap Overlaps + longLabel Paired exactly identical sequence on each side of a gap + bigDataUrl bbi/%s.gapOverlap.bb\n\n" "${asmId}" + fi + + if [ "${tanDupCount}" -gt 0 ]; then + printf " track tandemDups + parent tanDups on + shortLabel Tandem Dups + longLabel Paired exactly identical sequence survey over entire genome assembly + bigDataUrl bbi/%s.tandemDups.bb\n\n" "${asmId}" + fi + + $scriptDir/asmHubTanDups.pl $asmId $buildDir/html/$asmId.names.tab $buildDir/trackData > $buildDir/html/$asmId.tanDups.html + +fi + +export rmskCount=`(ls $buildDir/trackData/repeatMasker/bbi/${asmId}.rmsk.*.bb | wc -l) || true` + + # see if there are repeatMasker bb files export rmskCount=`(ls $buildDir/trackData/repeatMasker/bbi/${asmId}.rmsk.*.bb | wc -l) || true` if [ "${rmskCount}" -gt 0 ]; then printf "track repeatMasker compositeTrack on shortLabel RepeatMasker longLabel Repeating Elements by RepeatMasker group varRep visibility dense type bed 3 . noInherit on html html/%s.repeatMasker\n\n" "${asmId}" $scriptDir/asmHubRmsk.pl $asmId $buildDir/html/$asmId.names.tab $buildDir/trackData/repeatMasker/$asmId.rmsk.class.profile.txt > $buildDir/html/$asmId.repeatMasker.html fi