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/08 12:23] cespinos [BIER's pipeline] |
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 | ||
| Line 13: | Line 13: | ||
| * **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 located in the exome capture kit targets who had been realigned | * **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 located in the exome capture kit targets who had been realigned | * **%_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 ^ | + | ^ 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 | | | 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 | | | SGT077 | 27308034 | 27308034 | 27166800 | 99.48 | 28443473 | 52.08 | 24680760 | 45.19 | 20024427 | 36.66 | | ||
| Line 34: | 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 ^ | ||