Statistics
| Metric | R941 | R1041 |
|---|---|---|
| Estimate (average) | ||
| Estimated read accuracy | 94.54 | 97.49 |
| Read length | 2412.33 | 1538.32 |
| Read GC content | 49.45 | 49.38 |
| Observed (average) | ||
| Observed read accuracy | 94.94 | 97.34 |
| Observed read identification | 98.00 | 99.00 |
| Substitution proportion | 1.91 | 1.21 |
| Insertion proportion | 1.15 | 0.54 |
| Deletion proportion | 2.00 | 0.91 |
| Homopolymer accuracy (A) | 84.16 | 90.62 |
| Homopolymer accuracy (T) | 84.34 | 90.72 |
| Homopolymer accuracy (G) | 76.12 | 85.54 |
| Homopolymer accuracy (C) | 76.12 | 85.70 |
Estimate
Note: If the scale of accuracy is not suitable, please use the giraffe_plot function to replot.
giraffe_plot estimate_acc --input Estimated_information.txt --x_min 95 --x_max 100 --x_gap 1
Observe
Note: If the scale of accuracy is not suitable, please use the giraffe_plot function to replot.
giraffe_plot observe_acc --input Observed_information.txt --x_min 95 --x_max 100 --x_gap 1
Note: If the scale of proportion is not suitable, please use the giraffe_plot function to replot.
giraffe_plot observe_mismatch --input Observed_information.txt --y_max 5 --y_gap 1
Note: If the scale of accuracy is not suitable, please use the giraffe_plot function to replot.
giraffe_plot observe_homo --input Homoploymer_summary.txt --y_min 90 --y_max 100 --y_gap 2
GC bias
Note: If the scale of GC content is not suitable, please use the renormalization_sequencing_bias for normalzation and giraffe_plot for plotting.
renormalization_sequencing_bias -i S1_distribution.txt -l 30 -r 60 -o S1.txt
giraffe_plot gcbias --input new_gcbias.txt --x_min 20 --x_max 50 --x_gap 2