Statistics and figures for sequencing library: MM0581 (DefEndo)
Summary of input lanes for this library
Basic info is provided for each lane of data comprising this library
|Flowcell-Lane||Read Length||Quality Reads||Source Dir|
|61YEMAAXX_Lane7||2 x 51||68,611,871||/projects/analysis/analysis5/MM0581/61YEMAAXX_7/basecalls/|
|61YEMAAXX_Lane6||2 x 51||68,688,556||/projects/analysis/analysis5/MM0581/61YEMAAXX_6/basecalls/|
|61YEMAAXX_Lane5||2 x 51||68,962,188||/projects/analysis/analysis5/MM0581/61YEMAAXX_5/basecalls/|
|61YEMAAXX_Lane8||2 x 51||67,906,452||/projects/analysis/analysis5/MM0581/61YEMAAXX_8/basecalls/|
|61U0GAAXX_Lane7||2 x 51||65,495,773||/projects/analysis/analysis5/MM0581/61U0GAAXX_7/basecalls/|
Summary of sequence statistics for each lane
Basic statistics for each lane of sequence data are summarized below. 'Low quality' refers to reads with too many ambiguous bases (N's). 'Low complexity' refers to homopolymeric (e.g. polyA tails) or other sequences of low complexity identified by 'mdust'. 'R1' and 'R2' refers to the first and second read of a paired-end read.
|Flowcell-Lane||Paired Reads||Low Quality|
(R1 | R2)
(R1 | R2)
(R1 | R2)
|Total Bases||Passing Bases|
|61YETAAXX_Lane6||35,758,874||1.06% | 0.89%||0.00% | 0.00%||98.94% | 99.10%||2,074,014,692||98.96%|
|61YETAAXX_Lane5||36,612,692||1.03% | 0.87%||0.00% | 0.00%||98.97% | 99.13%||2,123,536,136||98.95%|
|61U0GAAXX_Lane7||33,267,242||0.95% | 1.09%||0.50% | 0.57%||98.54% | 98.33%||3,393,258,684||98.41%|
|61YEMAAXX_Lane8||35,760,773||4.36% | 4.84%||0.45% | 0.46%||95.19% | 94.70%||3,647,598,846||95.17%|
|61YEMAAXX_Lane5||36,299,625||4.30% | 4.74%||0.46% | 0.50%||95.23% | 94.75%||3,702,561,750||95.18%|
|61YEMAAXX_Lane6||36,145,905||4.28% | 4.73%||0.46% | 0.49%||95.26% | 94.78%||3,686,882,310||95.21%|
|61YEMAAXX_Lane7||36,124,222||4.34% | 4.80%||0.45% | 0.47%||95.21% | 94.72%||3,684,670,644||95.20%|
Summary of read assignments by assignment class
The assignment of reads to each read class is summarized below for the entire library. 'Unassigned' reads could not be assigned with high confidence to any known or predicted transcriptome or genome sequence. This does not mean that they had no significant similarity, only that they could not be assigned with high confidence. In this table, reads are summarized as individual reads (not paired reads)
|Read Class||Read Count||Percent of Total|
|Total read count||355,195,544||100%|
|Low quality (>1 N)||13,811,971||3.89%|
|Ensembl transcript (ambiguous)||17,001,030||4.79%|
|Novel exon junction||594,674||0.17%|
|Novel exon junction (ambiguous)||6,286||0.00%|
|Novel exon boundary extension||1,637,390||0.46%|
|Novel exon boundary extension (ambiguous)||49,149||0.01%|
|Repeat element (ambiguous)||542,657||0.15%|
Summary of mapping results for each lane
The mapping of reads to genome and transcriptome sequences is summarized below on a lane-by-lane basis. Only reads mapped with high confidence are summarized here. Reads mapped to each sequence type were assigned unambiguously. Reads mapping ambiguously (i.e. map equally well to multiple places) are summarized in the last column.
|Flowcell-Lane||Total||Transcript||Novel junction||Novel boundary||Intron||Intergenic||Repeat element||Ambiguous|
Summary of average coverage values by feature type
The grand average coverage observed for each type of sequence feature is summarized below. Average coverage is calculated as the cumulative number of mapped bases divided by the total number of base positions in the genome (for that type of sequence feature).
|Feature Type||Average Coverage||Total Base Count||Cumulative Coverage|
Summary of expressed events by feature type
The number of sequence features of each type detected above the level of background noise are summarized below. For each feature type, the total number of features is listed, followed by the subset of these features that are expressed versus not expressed above background.
|Feature Type||Feature Count||Expressed (%)||Not Expressed (%)|
|Gene||31,805||12,761 (40.12%)||19,044 (59.88%)|
|Transcript||49,040||13,073 (26.66%)||35,967 (73.34%)|
|ExonRegion||250,436||144,633 (57.75%)||105,803 (42.25%)|
|Junction||2,030,108||124,821 (6.15%)||1,905,287 (93.85%)|
|KnownJunction||201,049||123,913 (61.63%)||77,136 (38.37%)|
|NovelJunction||1,829,060||908 (0.05%)||1,828,152 (99.95%)|
|Boundary||398,581||11,860 (2.98%)||386,721 (97.02%)|
|KnownBoundary||14,969||5,995 (40.05%)||8,974 (59.95%)|
|NovelBoundary||383,613||5,865 (1.53%)||377,748 (98.47%)|
|Intron||193,082||2,283 (1.18%)||190,799 (98.82%)|
|ActiveIntronRegion||149,447||4,082 (2.73%)||145,365 (97.27%)|
|SilentIntronRegion||250,315||1,959 (0.78%)||248,356 (99.22%)|
|Intergenic||26,558||1,455 (5.48%)||25,103 (94.52%)|
|ActiveIntergenic||78,492||7,453 (9.50%)||71,039 (90.50%)|
|SilentIntergenic||72,440||2,291 (3.16%)||70,149 (96.84%)|
|3,530,304||326,671 (9.25%)||3,203,633 (90.75%)|
Estimates of signal-to-noise ratio
(Average coverage of exon regions / Average coverage of silent intron regions) = 108.38
(Average coverage of exon regions / Average coverage of silent intergenic regions) = 1448.84
Estimates of intronic and intergenic noise levels (95th percentiles of silent intron and intergenic regions)
95th percentile of silent intron regions for library (MM0581) is: 12.62 (log2 = 3.66)
95th percentile of silent intergenic regions for library (MM0581) is: 6.92 (log2 = 2.79)
Summary of library complexity - estimated by tag redundancy per million reads and compared to other libraries
Library complexity is calculated for the sequence library by randomly sampling 1 million reads and determining the number of unique and redundant reads within the pool. This sampling is repeated (with replacement) at least 3 times and average values across these samples is used. A 'redundant' read is one where both reads of a pair have either identical sequence (first figure) or mapping location (second figure) to at least one other read in the sampled pool. In each panel below, the value for the current library is indicated by a red dot and compared to other libraries, summarized as box and whisker plots. In the first panel (green), the percent of unique reads per million reads sampled is summarized. In the second panel (orange), the percent of redundant reads per million is summarized. In the third panel (yellow), the redundant reads are further examined to determine the number of distinct redundant reads. For example, if all redudant reads corresponded to a few sequences occuring many times this would result in a low number of distinct redundant reads. A 'good' library with high complexity should have high, low, and high values for panels one, two and three respectively.
Library complexity assessed by read sequence identity
Library complexity assessed by read sequence mapping positions
Distribution of relative positions of reads mapping within known transcripts (i.e. position bias test)
The frequency of read positions are plotted against the relative position within each transcript, where 0% is the 5' end of the transcript and 100% is the 3' end. Read positions were binned according to the size of transcripts they map within and the plot is produced for each bin (i.e. reads mapping to transcripts of 0-500 bp, 500-1000 bp, ..., 15000-20000 bp, and > 20000 bp).
Distribution of fragment size for paired-end mappings to transcripts
The distribution of fragment sizes is plotted for all paired-end reads in which both pairs could be aligned to the transcriptome. Pairs with reads aligning to different gene loci or different chromosomes are excluded. Fragment sizes with less than 0.01% of total reads are not plotted.
Distribution % of gene bases covered for each expressed gene (at various minimum depth levels)
The distribution of percent coverage levels for each gene is summarized as a box and whisker plot. These plots are produced for six minimum depth requirements. For example, in the first plot (>=1X), the percentage of bases covered to a depth of 1X or greater is determined for each gene. If a gene is completely covered by reads from the first base to the last, at a depth of 1X or greater, then this gene is given a value of 100%. The distribution of these percent values for all genes detected above background is summarized by the box and whisker plot.
Distribution of log2 raw expression values for each feature type
Box-and-whisker plots for log2 expression values for each feature type.
Distribution of log2 normalized expression values for each feature type
Box-and-whisker plots for normalized log2 expression values for each feature type.
Density scatter plot of exon region versus gene expression values
Density scatter plot of log2 expression values for exon regions versus corresponding gene expression values.
Density scatter plot of silent intron region versus gene expression values
Density scatter plot of log2 expression values for silent intron regions versus corresponding gene expression values.
Distribution of gene-by-gene expression cutoff values
Histogram depicting the distribution of expression cutoff values used for each gene. In order to be considered 'expressed above background', all features (genes, transcripts, exons, junctions, etc.) must be expressed above the level of INTERGENIC noise. This INTERGENIC cutoff is the 95th percentile of expression values for all Silent Intergenic Regions in this library and is depicted below as a dotted red line. The number of genes for which only INTERGENIC noise is considered is provided in the legend. For features within the boundaries of highly expressed genes, additional noise is expected due to the presence of un-processed RNA contamination. For this reason, a higher INTRAGENIC cutoff is determined. These are calculated by fitting a linear model to the 95th percentile of expression values for silent intronic regions. The INTRAGENIC cutoff for a gene is then determined by using the gene expression level and the coefficients of the model fit. The distribution of the resulting gene-by-gene cutoffs is depicted below as a histogram. The number of genes that required an INTRAGENIC cutoff is also indicated in the legend.
Histograms of expression values for each feature type
Histograms depicting distribution of log2 expression values for individual feature types. The 95th percentile of expression values for silent intergenic region is depicted as a dotted line on all plots.
Percentiles plot for expression of exon regions, silent intron regions and silent intergenic regions
Percentiles plot for exon region, silent intron region and silent intergenic region expression values. The 95th percentiles of intronic and intergenic distributions are depicted as colored, dotted lines
Cumulative distribution of mean Phred score for all reads
CDF plot depicting cumulative distribution of mean Phred score for all reads broken down by R1/R2 and sequencing lane