Filter hits that aligned concordantly exactly 1 times base on mpq>40
Hello,
I need to get all reads that aligned concordantly exactly 1 times base on mpq>40 from hisat2.
samtools view -b -q 40 -o output.bam alignments.bam
Will this command make a job?
• 2,205 views
•
link
updated
by
Brian Bushnell
•
written
by
User000
•
written
by
User000
1 answer
I assume "make a job" means "work". Yes, that will work.
• 0 views
•
link
written
by
Devon Ryan
Log in to answer this question.
More posts like this
-
featureCounts output summary assigned value higher than uniquely mapped reads from HISAT2
written by Prawesh •
-
Bowtie2 and Hisat2 read mapping comparisons
written by Allan •
-
sam file error
written by Chironex •
-
multiqc hisat2 transfer
written by smithkthedale •
-
Should we keep only the uniquely mapped reads for gene expression
written by User000
-
How to plot alignment statistics; how many reads are mapped to the genome?
written by anamaria
-
Why does Hisat2 report a different number of aligned pairs than samtools view?
written by adam.faranda
-
Align paired and unpaired reads simultaneously using Bowtie2?
written by xiaozhongzhiping •
-
Low fragment assignment RNAseq
written by said3427
-
Hisat2 versus Bowtie2
written by nalpas.nicolas •
