我用seqtk包将.fasta文件转换成.fastq后,用下列语句进行比对
bowtie2 -p 16 - hap_genome -U .fq -S .sam
后形成的sam文件用samtools view转换成.bam文件
出现的问题,如图
我换了fasta_to_fastq.pl进行文件转化再进行比对,而且比对不同的基因组也是同样的结果。我的sam的头文件如图
整个比对是在linux中进行的,在miniconda的环境中,samtools还是bowtie2都是用conda安装的。是我的原始测序数据有问题吗?网上找不到解决办法,恳请大家帮我想想问题出在哪里?非常感谢!