Kun:LabNotes/SingleCellExpr/2015-4-25: Difference between revisions
Jump to navigation
Jump to search
(Created page with "==Extract mapped reads from Bam files for IGV inspection== *I wrote a script for quick extraction of mapped sequencing reads near a list of genes from multiple data sets. This...") |
|||
Line 1: | Line 1: | ||
==Extract mapped reads from Bam files for IGV inspection== | ==Extract mapped reads from Bam files for IGV inspection== | ||
*I wrote a script for quick extraction of mapped sequencing reads near a list of genes from multiple data sets. This makes it easy to visually inspect the mapped reads with IGV, checking whether the genes/exons of interest have read coverage or not. | *I wrote a script for quick extraction of mapped sequencing reads near a list of genes from multiple data sets. This makes it easy to visually inspect the mapped reads with IGV, checking whether the genes/exons of interest have read coverage or not. | ||
extract_bam_for_selected_genes.pl gene_list data_set_folder output_folder | [[Media:extract_bam_for_selected_genes.txt|extract_bam_for_selected_genes.pl]] gene_list data_set_folder output_folder | ||
THREE ARGUMENTS: | THREE ARGUMENTS: | ||
gene_list: a plain text file containing the HUGO IDs for a list of genes, one gene per line | gene_list: a plain text file containing the HUGO IDs for a list of genes, one gene per line |
Revision as of 05:30, 26 April 2015
Extract mapped reads from Bam files for IGV inspection
- I wrote a script for quick extraction of mapped sequencing reads near a list of genes from multiple data sets. This makes it easy to visually inspect the mapped reads with IGV, checking whether the genes/exons of interest have read coverage or not.
extract_bam_for_selected_genes.pl gene_list data_set_folder output_folder THREE ARGUMENTS: gene_list: a plain text file containing the HUGO IDs for a list of genes, one gene per line data_set_folder: the folder that contains multiple mapped data set, each data set folder has a suffix "_mapped" output_folder: the folder where all extracted bam files will be placed.
Example: in /home/kunzhang/temp/brain_bulk_bams, run /home/kunzhang/RNAseq/SCAP/scripts/extract_bam_for_selected_genes.pl brain_selected_genes.txt \ /media/LTS_33T/RL_LTS33T/201404_201405_7Samples_BulkNucleiBatch1-20140623_Expt146/STAR .
- All the extracted bam files can be loaded into IGV for visual inspection.