samtools view for selecting the unmapped reads
Hi there,
I used Interactive Genome Viewer (IGV), to visualization of produced sam file from bowtie2, after loading of sorted-sam in IGV..from where and how I can select and identify the unmapped reads? Because by rolling mouse cursor around the each read, in box saying that mapped..
Don't be silent please and help me
Thank you in advance
1 answer
You won't see unmapped reads in IGV. If you want to extract them from a SAM/BAM/CRAM file, then use samtools (samtools view -f 4 foo.bam)
• 0 views
•
link
updated
by
Ram
•
written
by
Devon Ryan
•
written
by
Devon Ryan
Log in to answer this question.
More posts like this
-
Detect Transgene Insertion Site using SAM format
written by kspata
-
Error of SAMSequenceDictionary
written by niu.shengyong •
-
peak calling for medip seq analysis using macs
written by subhoshree.ghose •
-
IGV error: Error loading SAM header: For sequence 34711, text and binary have different lengths in …
written by EVR
-
How can I Annotate my miRNA data against MirBase?
written by k.kathirvel93
-
Bowtie2 --un-conc output empy
written by t2
-
--u option in bowtie2
written by zizigolu
-
Why sequence is unmapped by bwa but position is given in Sam file and also can be aligned by blat ?
written by ZheFrench
-
I loaded a BAM (RNA-seq) file into IGV but cant see anything!
written by jyalda2 •
-
Method used for reducing contigs IonTorrant 400bp [unpaired] (critique/ways to improve would be app…
written by bobo •
Genome browser is to view mapped reads against genome and other annotations. You can filter reads based on various flags. https://broadinstitute.github.io/picard/explain-flags.html
I went across the link you suggested, then what I should type for flag?
As Devon suggested below.
You mean, in cmd I should type the command Devon suggested?
Yes.
It would be better to switch to linux or something like Cygwin ( if you intended to use Windows only ) rather than trying out to do it on windows,
Yes you are right but I have internet connection problem in linux, that's why I'm using Windows
Anyway I typed this command but told failed to open foo.bam
I copied the
eg1.samproduced in bowtie2 folder in samtools folder and typed the command but got errorIf you have sam file, specify that input is sam using
-Soption. By default, samtools expects bam fileThank you, the same error
Can you give exact command, the folder path where the file is and error?
Sorry,
I extracted
samtools.ziptosamtools, copiedeg1.samproduced in bowtie2 folder to samtools folder and I typed in cmd:It's saying that:
Does it actually say
foo samor insteadfoo.sam? In the former case try specifying the correct name. In the latter case check the file permissions.You are right, saying
foo.samSorry, what is file permissions?
Given that you're in Berlin, perhaps you're more familiar with "Datei-Zugriffsrechte". Normally one right clicks on a file/directory and goes to "Preferences" (something along the lines of "Einstellungen" in the German version of Windows, if I remember correctly).
thank you, I am going to check