BAM index via samtools
In order to use genomic data in IGV I converted SAM files to BAM and them executed the index via samtools. I observed that when generating the index from the same BAM file the output could have different sizes. To test I inputed the same genome three times, with diferent index in IGV and verified slighth diferences as shown in the image. Does anyone knows if i have done something wrong, or if that is a normal thing? Thanks in advance!
• 1,442 views
•
link
written
by
carolina.santiago.t •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Discrepancy between positions on IGV and NCBI's PrimerBLAST
written by Rozita •
-
No variants found ( have problem in the handbook)
written by quanyu •
-
GATK workflow for de novo transcriptome assembly of non-model organism
written by poecile.pal •
-
Read count difference between STAR and Bowtie
written by hothriananya
-
IGV error: Error loading SAM header: For sequence 34711, text and binary have different lengths in …
written by EVR
-
SNP calling through bcftools
written by palak.siddhu •
-
Are the identified SNPs in the given plot paralogs or positive variants (alleles)?
written by kirannbishwa01
-
IGV not reading BAM file
written by novice
-
I loaded a BAM (RNA-seq) file into IGV but cant see anything!
written by jyalda2 •
-
454 Sequencing With Bwa Mapping
written by Gabriel Wajnberg
Size of files is never a good criteria to judge validity of a command/action. If IGV is not complaining about the indexes then they are fine.
Please use: How to add images to a Biostars post
Thank you for your help! I can't understand why, when using the same input file and the same command, the size of the output can be different. That´s messing with my head. Do you know the reason why this happens?
Can you show us what kind of difference are we looking at? Many NGS data analysis program produce "non-deterministic" output (i.e you may get slightly different results each time you run a dataset). In practice this is not a big problem. Some of the programs will allow you to set a
seedso you get exactly the same result each time. I don't know ifsamtools indexis one of those operations.Thank you, once again, for your help! I am really sorry but my question was the result of a mistake, that i only realysed now. I was using the computer at distance, i don't know if it was some kind of delay but the truth is that now i'm looking at the same files at the computer and they do have the exact same size! I'm really sorry! Thank you once again!