Calling SNPs INDELs with SAMtools BCFtools employs SAMtools 1.7 mpileup and BCFtools 1.8 call to call SNPs and INDELs
Test data for this app appears directly in the Discovery Environment in the Data window under Community Data -> iplantcollaborative -> example_data -> Samtools_mpileup.
Use ex1.bam from the directory above as bam input and ex1.fa from the directory above as reference sequence.
When the app is run in the Discovery Environment, use the following parameters with the above input file(s) to get the output provided in the next section below.
- Use these parameters within the DE app interface:
- Set the mpileup output filename.
- Set the BCFtools call output filename and set the output file type to uncompressed VCF.
- Leave all other parameters as default.
The output file in the example data generated with these settings is bcftools_1.8_call_output.vcf as output.