Hello,
I have a very peculiar problem. I used samtools to get SNPs using the mpileup command. The commands that i used are as follows:
samtools mpileup -uf cpp34_contigs_nr_prinseq_good_NAF4.fasta aln-sorted.bam | bcftools view -bvcg - > var.raw.bc
bcftools view var.raw.bcf > var.flt.vcf
The .vcf file that i got had the follwing results:
#CHROM POS ID REF ALT QUAL
PT_34_209 159 . T C 7.8
PT_34_209 164 . C T 7.8
PT_34_209 655 . A G 7.8
PT_34_209 658 . G C 7.8
PT_34_209 659 . C G 7.8
PT_34_238 533 . T C 7.8
PT_34_238 600 . T C 7.8
PT_34_238 786 . A T 7.8
PT_34_238 789 . A C 7.8
PT_34_238 790 . T C 7.8
PT_34_238 791 . T G 7.8
PT_34_238 792 . T C 7.8
PT_34_238 793 . G T 7.8
PT_34_238 795 . T G 7.8
PT_34_258 561 . A G 7.8
PT_34_549 520 . T C 7.8
Could anybody please help me with this? and let me know how i can correct the problem between the reference column and the quality column?
I have a very peculiar problem. I used samtools to get SNPs using the mpileup command. The commands that i used are as follows:
samtools mpileup -uf cpp34_contigs_nr_prinseq_good_NAF4.fasta aln-sorted.bam | bcftools view -bvcg - > var.raw.bc
bcftools view var.raw.bcf > var.flt.vcf
The .vcf file that i got had the follwing results:
#CHROM POS ID REF ALT QUAL
PT_34_209 159 . T C 7.8
PT_34_209 164 . C T 7.8
PT_34_209 655 . A G 7.8
PT_34_209 658 . G C 7.8
PT_34_209 659 . C G 7.8
PT_34_238 533 . T C 7.8
PT_34_238 600 . T C 7.8
PT_34_238 786 . A T 7.8
PT_34_238 789 . A C 7.8
PT_34_238 790 . T C 7.8
PT_34_238 791 . T G 7.8
PT_34_238 792 . T C 7.8
PT_34_238 793 . G T 7.8
PT_34_238 795 . T G 7.8
PT_34_258 561 . A G 7.8
PT_34_549 520 . T C 7.8
Could anybody please help me with this? and let me know how i can correct the problem between the reference column and the quality column?
Comment