This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
espinos:results.primary.analysis [2015/08/06 19:43] cespinos |
espinos:results.primary.analysis [2015/09/27 19:34] (current) cespinos |
||
|---|---|---|---|
| Line 6: | Line 6: | ||
| ==== BIER's pipeline ==== | ==== BIER's pipeline ==== | ||
| - | In this pipeline, there are four stages where the number of reads decreases (mapping, filter by mapping quality, remove duplicates and intervals realignment). This {{:espinos:count_reads_bier_eng.txt|table}} shows reads remaining after these stages. | + | In this pipeline, there are four stages where the number of reads decreases (mapping, filter by mapping quality, remove duplicates and intervals realignment). This table shows reads remaining after these stages. |
| - | * **N_reads_forward and reverse**: initial number of reads forward and reverse obtained in the exome sequencing process. | + | * **N_reads_forward and reverse**: initial number of reads forward and reverse obtained in the exome sequencing process |
| - | * **N_mapped_read_pairs**: number of read pairs mapped to the human genome reference. | + | * **N_mapped_read_pairs**: number of read pairs mapped to the human genome reference |
| - | * **%_mapped_read_pairs**: percentage of initial read pairs mapped to the human genome reference. | + | * **%_mapped_read_pairs**: percentage of initial read pairs mapped to the human genome reference |
| - | * **N_mapped_reads_mapq>10**: number of mapped reads whose mapping quality (mapq) is higher than 10. | + | * **N_mapped_reads_mapq>10**: number of mapped reads whose mapping quality (mapq) is higher than 10 |
| - | * **%_mapped_reads_mapq>10**: percentage of initial mapped reads whose mapping quality (mapq) is higher than 10. | + | * **%_mapped_reads_mapq>10**: percentage of initial mapped reads whose mapping quality (mapq) is higher than 10 |
| - | * **N_reads_single_hit**: number of reads uniquely mapped to the human genome reference. | + | * **N_reads_single_hit**: number of reads uniquely mapped to the human genome reference without duplicates |
| - | * **%_reads_single_hit**: percentage of initial reads uniquely mapped to the human genome reference. | + | * **%_reads_single_hit**: percentage of initial reads uniquely mapped to the human genome reference without duplicates |
| - | * **N_reads_single_hit_realigned**: number of reads realigned in the exome capture kit targets. | + | * **N_reads_single_hit_realigned**: number of reads located in the exome capture kit targets who had been realigned |
| - | * **%_reads_single_hit_realigned**: percentage of initial reads realigned in the exome capture kit targets. | + | * **%_reads_single_hit_realigned**: percentage of initial reads located in the exome capture kit targets who had been realigned |
| + | |||
| + | ^ Sample ^ N_reads forward ^ N_reads reverse ^ N_mapped read_pairs ^ %_mapped read_pairs ^ N_mapped reads mapq>10 ^ %_mapped reads mapq>10 ^ N_reads single_hit ^ %_reads single_hit ^ N_reads single_hit realigned ^ %_reads single_hit realigned ^ | ||
| + | | SGT038 | 31471997 | 31471997 | 31304563 | 99.47 | 32844926 | 52.18 | 28115216 | 44.67 | 22724178 | 36.10 | | ||
| + | | SGT077 | 27308034 | 27308034 | 27166800 | 99.48 | 28443473 | 52.08 | 24680760 | 45.19 | 20024427 | 36.66 | | ||
| + | | SGT161 | 27566691 | 27566691 | 27429566 | 99.50 | 28775620 | 52.19 | 24997317 | 45.34 | 20213242 | 36.66 | | ||
| + | | SGT187 | 29730857 | 29730857 | 29593092 | 99.54 | 31001122 | 52.14 | 26579044 | 44.70 | 21530503 | 36.21 | | ||
| + | | SGT230 | 30415770 | 30415770 | 30296046 | 99.61 | 31712462 | 52.13 | 27084895 | 44.52 | 22377265 | 36.79 | | ||
| + | | SGT238 | 29472514 | 29472514 | 29333890 | 99.53 | 30770469 | 52.20 | 26406037 | 44.80 | 21543680 | 36.55 | | ||
| + | | SGT241 | 29513223 | 29513223 | 29365739 | 99.50 | 30803708 | 52.19 | 26174909 | 44.34 | 21579194 | 36.56 | | ||
| + | | SGT274 | 30394832 | 30394832 | 30268625 | 99.58 | 31699515 | 52.15 | 27122226 | 44.62 | 22387721 | 36.83 | | ||
| + | |||
| Line 23: | Line 34: | ||
| In contrast with BIER's pipeline, there are only two stages where the number of reads decreases (mapping and remove duplicates). | In contrast with BIER's pipeline, there are only two stages where the number of reads decreases (mapping and remove duplicates). | ||
| - | * **N_reads_forward and reverse**: initial number of reads forward and reverse obtained in the exome sequencing process. | + | * **N_reads_forward and reverse**: initial number of reads forward and reverse obtained in the exome sequencing process |
| - | * **N_mapped_read_pairs**: number of read pairs mapped to the human genome reference. | + | * **N_mapped_read_pairs**: number of read pairs mapped to the human genome reference |
| - | * **%_mapped_read_pairs**: percentage of initial read pairs mapped to the human genome reference. | + | * **%_mapped_read_pairs**: percentage of initial read pairs mapped to the human genome reference |
| - | * **N_read_pairs_single_hit**: number of read pairs uniquely mapped to the human genome reference. | + | * **N_read_pairs_single_hit**: number of read pairs uniquely mapped to the human genome reference remaining after removing duplicates |
| - | * **%_read_pairs_single_hit**: percentage of initial read pairs uniquely mapped to the human genome reference. | + | * **%_read_pairs_single_hit**: percentage of initial read pairs uniquely mapped to the human genome reference remaining after removing duplicates |
| - | ^Sample^N_reads_forward^N_reads_reverse^N_mapped_read_pairs^%_mapped_read_pairs^N_read_pairs_single_hit^%_reads_pairs_single_hit^ | + | ^ Sample ^ N_reads_forward ^ N_reads_reverse ^ N_mapped read_pairs ^ %_mapped read_pairs ^ N_read_pairs single_hit ^ %_reads_pairs single_hit ^ |
| | SGT038 | 31471997 | 31471997 | 27098380 | 86.10 | 26533415 | 84.31 | | | SGT038 | 31471997 | 31471997 | 27098380 | 86.10 | 26533415 | 84.31 | | ||
| | SGT077 | 27308034 | 27308034 | 23450991 | 85.88 | 22904031 | 83.87 | | | SGT077 | 27308034 | 27308034 | 23450991 | 85.88 | 22904031 | 83.87 | | ||