hi, I wonder what’s the difference when metaPhlan deal with paired-end reads and single reads? I run my paired-end data with wrong command (“metaphlan fastq1 fastq2 -o”) ,which end with a re-write bowtie2out to my raw data(a bowtie2out named fastq2 rewrite the raw fastq2), and a result which actually result from only one file in paired end. my question is that, if it’s necessary to use only one file in pared-end reads for metaPhlan to qualified the tax abundance instead of two?
Hi @jiazhong28
For running metaphlan with more than 1 FASTq file, you should specify the files separated by comma (wo spaces) e.g. metaphlan fastq1,fast2 ...
thanks for your reply, i’m still wondering that if it necessary for me to simply regared paired-end file as doubled results from just fastq1 only, since that metaphlan do not use paired information for abundabce caculation?
Even if MetaPhlAn does not use the paired-end information, markers not mapping to one of the ends could still be mapping to the other read pair when the mapping region is one of the ends of the marker. So I will still recommend using both pairs when executing metaphlan