This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Would this level of base pair mismatch suggest the data is inaccurate?

On viewing a bam file there frequently appears to be a large amount of base-pair mismatch such as this. Should this data be disregarded or is it in any way possible for this to be accurate data?

enter image description here

enter image description here

accuracy data base-pair mismatch alignment
• 834 views • link
updated by GenoMax GenoMax • written by Lou11 • Lou11

Are you using the right reference (genome build/organism) to view the data in IGV? It is possible that your sample has several changes compared to the reference (windows above are 25kb).

• 0 views • link
written by GenoMax GenoMax

I believe it is the right reference. The alternative (hg19) gives a hugely worse-looking result.

• 0 views • link
written by Lou11 • Lou11

If you did not create the BAM yourself, check the header of the BAM for clues to see what genome build was used. Generally the alignment command should give you an indication. Use something like.

samtools view -H your.bam | grep "CL"
• 0 views • link
written by GenoMax GenoMax

0 answers

No answers yet.

Log in to answer this question.

More posts like this