This is a test version of Biostars. For the public version, visit https://www.biostars.org.
why hisat2 mapping sam file has many sam flags

Hi, I used hisat2 to map my reads to hg19 genome, and it seems there's something wrong with this sam file. The @PG header shows this:

@PG     ID:hisat2       PN:hisat2       VN:     CL:"envs/rna/bin/hisat2-align-s --wrapper basic-0 -p 8 --dta -x GRCh37.p13.genome.clean.fa -S 1.sam -1 1.non_rRNA_paired.fq -2 2.non_rRNA_paired.fq"

and the samtools flag in this sam file has many different numbers: 113 ,129,133,137,145 ,147,153,161,163,177,321,329,337,339,345,353,355,369,385,393,401,403,409,417,419,433,65,69,73,81,83,89,97,99. But usually the double-end sequencing files only has 99,147 or something else, but not so many samtools flag numbers, so what's wrong here?

sam rna-seq sequencing hisat2 genome

0 answers

No answers yet.

Log in to answer this question.