MIG statisticsΒΆ
Description
A script to generate consensus coverage statistics, i.e. molecular identifier group (MIG) size distribution.
Usage
General:
java -jar migec.jar Histogram checkout/ histogram/
Running this script will generate several files in histogram folder, the one important for basic data processing is overseq.txt. The header of table contains MIG sizes (in log2 scale), while each row corresponds to a de-multiplexed sample contains the number of reads in MIGs of a given size (cumulative abundance).
For a decent dataset the plot of cumulative abundance display a small peak at MIG size of 1 that could be attributed to erroneous MIGs and has an exponential decline, and a clear peak at MIG size of 10+ containing amplified MIGs. Those erroneous MIGs could arise as experimental artifacts, however the most common reason for their presence is an error event in UMI sequence itself. Note that the latter is only valid when number of distinct UMIs is far lower than theoretically possible UMI diversity (e.g. 4^12 for 12-letter UMI regions)!
MIG size cutoff in Assemble should be set to dissect erroneous MIGs while retaining amplified ones. If peaks overlap collision filtering should be considered.
A simple plotting routine written in R can facilitate visualization of MIG size distributions, available here.
Important
The --only-first-read
option should be used if it is also
specified in Assemble routine for consensus coverage estimates
to be concordant.