I am also trying to check the coverage at each position of my reference using the CoverageBySample tool (with and without the –L argument):
java -Xmx30g -jar GenomeAnalysisTK.jar \ -T UnifiedGenotyper \ –T CoverageBySample \ –R ref.fasta \ -I input.bam \ -o output.cov\
The output (below) is giving the right coverage but without the positions on the reference and also skipping all positions with no coverage. Is there any way to get these positions in the output file?
eo78 10 eo78 10 eo78 10 eo78 10 eo78 10 eo78 11 eo78 12 eo78 12 eo78 12
UnifiedGenotyper- Can it be used to call SNP's in bacterial genomes