How to address "*" as RNAME in SAM file
I'm making a program in java for alignment visulaization. The input is a ref genome and a sam file with the alignment. I'm using RNAME field in the SAM file to see to which sequence does the alingment belong to. How to address when there is "" instead of a RNAME. In my SAM file for Ecoli alignemnt all RNAME fields are "" so I've added into my code that if there is only one sequence in the referente genome file the program takes the sequence name from the header as RNAME. I don't know if this aproach is correct and I don't know what to do when the referente genome file has more than one sequence. Any advice/explanation would be appreciated.
• 773 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
Matej •
•
written
by
Matej •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to utilize the bai index file while parsing a BAM file
written by matejm469 •
-
Parsing SAM files with java
written by Matej •
-
Find and replace the RNAME field in BAM/SAM file
written by Begonia_pavonina •
-
When can you called a genome closed/no longer draft?
written by kristensen1234 •
-
variant calling ecoli
written by lima80in •
-
bwa XA field: What is the procedure for the binning of all but one alternative alignments into XA?
written by izaak.coleman1 •
-
Redirecting the Log.txt file from a STAR indexing run
written by biomagician
-
htseq-count error RNAME == '*' although flag bit &0x0004 cleared"
written by super •
-
Tutorial: SAM files processing and variants calling in bacterial genomes
written by shilparaopradeep •
-
Tutorial: Sam File Format - Lesser Known Tips And Tricks
written by Ya •
don't reinvent the wheel and use htsjdk. https://github.com/samtools/htsjdk/