samtools having to use -A input option

0

Entering edit mode

4.7 years ago

williamsbrian5064 ▴ 510

Hi,

I am using samtools to make a bcf file and then a vcf file. I am using the following code

samtools mpileup --redo-BAQ --min-BQ 30 --per-sample-mF -C 50 --output-tags DP,DV,DPR,INFO/DPR,DP4,SP -A -u -f Felis_catus.Felis_catus_9.0.dna_sm.toplevel.fa --BCF C04893_L001.sorted.dedup.bam > C04893.bcf

I have used samtools on a number of samples and never had to use the -A input option. For some reason, one of my samples wouldn't start making a the bcf file without using the -A input option. I was curious why what happened all of the sudden? I was wondering if anyone had a similar issue? Is there a way to see if there is something weird about the bam file? Would samtools flagstat tell me anything?

assembly next-gen genome alignment • 712 views

ADD COMMENT • link 4.7 years ago by williamsbrian5064 ▴ 510

0

Entering edit mode

Hello,

could you please show the output of samtools flagstat for a "normal" sample andbthe "weird" one?

fin swimmer

ADD REPLY • link 4.7 years ago by finswimmer 16k

0

Entering edit mode

The simple explanation is that there is something really different about that fastq file. Did you look up what -A meant before you asked this question?

ADD REPLY • link 4.7 years ago by swbarnes2 14k

Login before adding your answer.