...的snp.unknown_multianno.vcf文件绘制时有报错信息，Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, :
line 46 did not have 209 elements
Calls: read.table -> scan，应该是生成的vcf格式有问题，但是用一开始过滤完的snp.clean.vcf.gz...
...k/tmp/tmpW7bSxCiAuQqVCErU3a2S.txt"
Traceback (most recent call last):
File "/biosoft/miniconda/envs/qiime1/bin/assign_taxonomy.py", line 417, in <module>
main()
File "/biosoft/miniconda/envs/qiime1/bin/assign_taxonomy.py", line 394, in main
log_path=log_path)
File "/...
not well-formed (invalid token) at line 177, column 44, byte 5589 at /biosoft/miniconda/envs/genefamilyenv/lib/meme-5.1.0/perl/CheckingSAX.pm line 46.Warning: meme_xml_to_html exited abnormally and may have failed to create HTML output.
...qc/cleandata/${i}.bowtie.map.log # Filter and output as SAM(or BAM),-f 4 to Remove reads from the alignment, # if you want to keep the -F on the alignment (12 for both ends), add -b before the output to output in BAM format samtools view --threads 6 \ -f 12 -h \ $workdir/1.host_del_and_qc/cleand...
...文献：
Carsten Lund Pedersen, This simple tool shows you how to choose your mentors, natureindex, 28 August 2019.
生物女学霸，一个自称学霸其实很渣的生物汪，努力将有趣的、有用的、有血有肉的科研那些事儿呈现给大家，关注一个好不啦~