Reviewing a fastq_mergepairs report to check for problems 信息搬運(yùn)

信息來(lái)源https://www.drive5.com/usearch/manual/merge_report.html


Below is an example report produced by the -report option of fastq_mergepairs. This information is also shown on the terminal (standard error output stream). The options -fastq_minmergelen 230 -fastq_maxmergelen 270 were used because these are 2 x 250 reads of amplicons generated by three different primer pairs including V4, V3-V4 and V4-V5. Using the?length range?230 to 270 selects the V4 reads.

For each parameter used in pre-processing, alignment, merging and filtering the report shows how many pairs were successfully processed or discarded. The parameter value is also shown. For example, 7.9M reads (58%) were discarded because the alignment had >5 mismatches (parameter set by the fastq_maxdiffs option).

Here we have long overlaps, shown by the mean alignment length of 248. Mis-alignments are therefore very unlikely, and it would be reasonable to increase the -fastq_maxdiffs and -fastq_maxdiffspct values to increase the number of merged pairs. Quality filtering will take care of discarding reads where many mismatches induce a large number of expected errors. This doesn't necessarily happen -- e.g., if low quality base calls in R2 are mismatches against high-quality base calls in R1 then the merged Q scores can still be high.

If you think that too many reads are being discarded, then you can?use the tabbedout file to investigate further.


?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書(shū)系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

  • Introduction What is Bowtie 2? Bowtie 2 is an ultrafast a...
    wzz閱讀 6,186評(píng)論 0 5
  • rljs by sennchi Timeline of History Part One The Cognitiv...
    sennchi閱讀 7,871評(píng)論 0 10
  • All the following information come from www.drive5.com, I...
    代號(hào)北極能閱讀 1,168評(píng)論 0 0
  • 揮灑汗水的2018終于過(guò)去了,在廈門(mén)迎來(lái)了2019。難得的元旦三天假期,節(jié)前兩天才確認(rèn)行程,又是說(shuō)走就走。看著公司...
    Helen煥煥閱讀 236評(píng)論 0 0
  • 睡前。 寶寶指指右邊的爸爸 爸爸是男生 指指左邊的媽媽 媽媽是女生 那你呢? 我是單身! 天哪,誰(shuí)教你的?!這……
    大樂(lè)之野閱讀 343評(píng)論 0 0

友情鏈接更多精彩內(nèi)容