Is there is Python module / library deals with the flag field in SAM files?
I am writing some python script for sam file , and I want to get the information from the flag file about the sequence, Is there any module or lib that can do it or I need to write mine?
if there is not previous code to do it what is the pseudo code that deals with it shall I convert from int to hex then bin and use some tabular to get the info or what?
Thanks
2 answers
pysam can do that for you. An AlignedSegment object has fields such as is_duplicate, is_proper_pair, is_read1, etc.. That makes things pretty convenient.
• 0 views
•
link
written
by
Devon Ryan
Log in to answer this question.
More posts like this
-
Write a file with CWL
written by hoseintoossi •
-
Tutorial: How to select reads with specific flag using samtools
written by imlituan •
-
best practices to use GATK with PacBio
written by Medhat
-
How to filter unique reads from sam file?
written by Medhat
-
How to filter unique reads from sam file?
written by Medhat
-
Best way to extract gene (CDS) sequences from genome assemblies.
written by bioinformaticsscientist •
-
how to exclude sequence from sam file?
written by Medhat
-
Trouble Compiling C++ Code With Bamtools
written by user
-
Attribute Error Using Biopython Bio.Searchio.Convert
written by charles.bridges
-
How To Edit A Bam File With Bio::Db::Bam
written by Owen S.
Have you tried pysam library? I think the sam/bam object has an attribute called
flag.Three replies of "try pysam" at once must be a new record.
is there any badge for it ? :-P
We should request one :)
It means, I need to make my homework with more concentration before I post a question
You just need more coffee :)
:D I agree