Hi
I have some .sam files as output of bowtie run. i want to read them into Rstudio. i searched for it but all i can find is about importing .bam files. some suggest converting .sam into .bam. is it the only way? and is this convertion lossless?
thanks
•
written
by
Amirosein
1 answer
Sam file is a text file. So you can read it directly in R. But there is no point.
Best way is to convert sam to bam and then use a library like RSamtools to read the bam files.
You could read it. Its a text file. But to have an easy and quick access the information it has, better to use some libraries, rather than writing something from scratch. Converting to bam file has many other advantages, which is not the question here.
This might be useful, http://davetang.org/wiki/tiki-index.php?page=SAMTools
Log in to answer this question.
-
Which command line can run dockedtopdbq in autodock4 (including mgltools 1.5.7)
written by nonaddldy •
-
Is it okay to convert bam files to fastq and then getting seqkit results?
written by Baran •
-
samtools manipulation of bowtie output with suppressed fields
written by Sjannie Lefevre •
-
parallel bwa mem with compressed fastq.bz2 files
written by Amirosein
-
SAMTOOLS flagstat outputs
written by jt358 •
-
How to check new variants for functional effects?
written by Amirosein
-
Variant calling for sanger reads
written by Amirosein
-
Empty SAM files with sam-dump
written by adamlangenbucher •
-
Basic SAMtools tutorial
written by statfa
-
Best practice in vcf calling - From fastq to vcf ?
written by amir.jebelli •
The main question is "why would you want to do that?" Sam and bam contain the same information, bam is just binary of sam.