thank you for your replying！
• 0 views
•
link
written
by
minghuiguo448 •
Question: how to get low quality read from sam file
Hello, l am a starter.I used hisat2 to generate a sam file and have extracted unmapping reads. l want to extract reads with a mapping quality lower than 10.My instructor reminded me to think about how to get the mapping quality.but I don't know how to extract it. Thanks for your help.
• 1,679 views
•
link
updated
by
ATpoint
•
written
by
minghuiguo448 •
•
written
by
minghuiguo448 •
2 answers
For example use sambamba:
sambamba view -c -F "mapping_quality < 10" --sam-input -f sam -o filtered.sam in.sam
https://lomereiter.github.io/sambamba/docs/sambamba-view.html
samtools view is able to filter for "quality greater than". But for "lower than", you could use awk:
samtools view -h in.bam | awk -F '\t' '$0 ~ /^@/ || $5< 10'
• 0 views
•
link
written
by
Pierre Lindenbaum
thank you for your replying！
• 0 views
•
link
written
by
minghuiguo448 •
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 •
-
Adding @RG header line to bam and sam files in hisat2
written by Nemo •
-
Question about a paper called Reference flow: reducing reference bias using multiple population ge…
written by minghuiguo448 •
-
featureCounts: extremely low rate of 'Successfully assigned alignments'
written by Yi Chu •
-
How to extract uniquely mapped reads in Bowtie2
written by Ankit
-
Extract unique reads from SAM file
written by niu.shengyong •
-
Mapping RNA-Seq data on genome from another species
written by tlorin
-
Puzzle on size of Bam file
written by devinliao0918 •
-
Filtration Of Reads With Length Lower Than 30 From Bam
written by filipzembol
-
Extract Regions From A Bam To A Bed File Where The Average Mapping Quality Is Below 30.
written by William
samfile is a text file and one of the columns codes for map qualities. Filter the data by map quality column and validate with standard tools such as samtools.