It makes sense to implement the extensive functionality offered by IVar to get from aligned BAM to consensus sequence. Varscan 2 and BCFTools can still be used for metagenomics data too. We are currently only using it to trim primer sequences so would need to add the other commands as described in the MANUAL.
Also, need to make sure that the parameters we use are applicable to Illumina data too.
It makes sense to implement the extensive functionality offered by
IVar
to get from aligned BAM to consensus sequence.Varscan 2
andBCFTools
can still be used for metagenomics data too. We are currently only using it to trim primer sequences so would need to add the other commands as described in the MANUAL.Also, need to make sure that the parameters we use are applicable to Illumina data too.