mutations file
 

Hello all,

I just finished generating BAM file of resequenced tomato line from paired end (reads size 100bp).
and view it through IGV (aligned to reference genome). Do you familiar with any way to load all mutations that were found in the resequenced line to one file instead of looking for them in IGV?

Thanks in advance,

Itay

When you refer to "mutations file" I assume you are referring to SNP's? You can use samtools to call SNP (and INDEL's too). Here is a simple guide (skip to the "Basic Variant Calling" section): https://wikis.utexas.edu/display/bio...using+SAMtools

https://www.ebi.ac.uk/training/sites...np_calling.pdf