Question: RNAseq read mapping statistics from sam file
gravatar for Bioinfonext
4 months ago by
Bioinfonext200 wrote:


I used below script to run Hisat2 for multiple pair end read files to map on reference genome separately. I got the sam file but it did not generate the statistics of read mapping, is it still possible to get read mapping statistics from sam files?

#SBATCH --job-name=mpi_job_test      # Job name
#SBATCH --mail-type=END,FAIL         # Mail events (NONE, BEGIN, END, FAIL, ALL)     
#SBATCH --ntasks=8                  # Number of MPI ranks
#SBATCH --cpus-per-task=1            # Number of cores per MPI rank 
#SBATCH --nodes=1                    # Number of nodes
#SBATCH --time=200:05:00              # Time limit hrs:min:sec

#SBATCH --partition=lowpri


module add hisat2/2.1.0

hisat2 -x chrX_tran -1 Leaf_T1_F_R10_S1_L001_R1_paired.fq.gz -2 Leaf_T1_F_R10_S1_L001_R2_paired.fq.gz -S Leaf_T1_F_R10_S1_L001.sam

hisat2 -x chrX_tran -1 Leaf_T1_F_R2_S5_L001_R1_paired.fq.gz -2 Leaf_T1_F_R2_S5_L001_R2_paired.fq.gz -S Leaf_T1_F_R2_S5_L001.sam

hisat2 -x chrX_tran -1 Leaf_T1_F_R3_S4_L001_R1_paired.fq.gz -2 Leaf_T1_F_R3_S4_L001_R2_paired.fq.gz -S Leaf_T1_F_R3_S4_L001.sam

Many thanks nabiyogesh

bash rna-seq hisat2 • 142 views
ADD COMMENTlink written 4 months ago by Bioinfonext200

Please use the search function:

How To Get Mapping Statistics On Bam File

Rna-Seq Read Mapping Stats/Metrics

ADD REPLYlink written 4 months ago by ATpoint31k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1751 users visited in the last hour