Aggregated Frequencies Of Transcription Initiations Observed In Fantom5 Cage Data On Grch38, Including Alignments With Low Mapping Qualities

<strong>Overview</strong> Aligned reads of the FANTOM5 CAGE data have been used after filtering (ones with mapping quality less than 20 or percent identity less than 85% were discarded) for general purpose, resulting in the data set consisting of only the reads aligned with confidence. The filtering process made possible to interpret the data without ambiguity, however it also limited interpretation of paralogous or duplicated regions within the genome. Here all of the 5'-ends of the CAGE read alignments, including the ones with low mapping quality, were counted. The counts in the individual profiles were aggregated and summed up. <strong>Special usage note</strong> As noted above, this data derived from the alignments with low mapping qualities, as well as the ones with high mapping qualities. The result has to be examined very carefully: observations on the genome does not support transcription initiation with confidence, and even absence of such observation does not support silence of transcription with confidence. For example, file size on the forward strand is substantially larger than the one on the reverse strand, which is likely caused by an arbitrary preference of the alignment process. It does not mean transcription happens more frequently on the forward strand. Interpretation has to be made always in comparison with the standard data (BED files under http://fantom.gsc.riken.jp/5/datafiles/reprocessed/hg38_v4/basic/ or bigWig files under http://fantom.gsc.riken.jp/5/datahub/hg38/reads/). <strong>Data files</strong> The resulting data files are formatted as bigWig (https://genome.ucsc.edu/FAQ/FAQformat.html#format6.1). '*.fwd.bw' and '*.rev.bw' represent forward and reverse strand on the genome, respectively. <strong>Methods</strong> The BAM files under http://fantom.gsc.riken.jp/5/datafiles/reprocessed/hg38_v4/basic/ were subjected to 5'-end counting by bedtools v2.27.1 (https://github.com/arq5x/bedtools2), followed by conversion into bigWig with jksrc v357 (http://hgdownload.cse.ucsc.edu/admin/).