扩增子分析解读2提取barcode 质控及样品拆分 切除扩增引物

  • A+
所属分类:生物信息学

This is cutadapt 1.14 with Python 3.6.1Command line parameters: -g AACMGGATTAGATACCCKG -a GGAAGGTGGGGATGACGT -e 0.15 -m 300 --discard-untrimmed temp/PE250_split/seqs.fna -o temp/PE250_P5.faTrimming 2 adapters with at most 15.0% errors in single-end mode ...Finished in 73.83 s (58 us/read; 1.04 M reads/minute).
=== Summary ===
Total reads processed: 1,277,436Reads with adapters: 1,277,194 (100.0%)Reads that were too short: 8,849 (0.7%)Reads written (passing filters): 1,268,345 (99.3%)
Total basepairs processed: 522,379,897 bpTotal written (filtered): 495,607,409 bp (94.9%)
=== Adapter 1 ===
Sequence: GGAAGGTGGGGATGACGT; Type: regular 3'; Length: 18; Trimmed: 202757 times.
No. of allowed errors:0-5 bp: 0; 6-12 bp: 1; 13-18 bp: 2
Bases preceding removed adapters: A: 96.3% C: 1.5% G: 0.8% T: 1.3% none/other: 0.0%WARNING: The adapter is preceded by "A" extremely often. The provided adapter sequence may be incomplete. To fix the problem, add "A" to the beginning of the adapter sequence.
Overview of removed sequenceslength count expect max.err error counts3 3 19959.9 0 34 4 4990.0 0 46 2 311.9 0 28 1 19.5 1 111 1 0.3 1 113 1 0.0 1 115 9 0.0 2 917 42 0.0 2 4218 202626 0.0 2 20262619 56 0.0 2 5620 1 0.0 2 121 1 0.0 2 132 1 0.0 2 138 1 0.0 2 139 1 0.0 2 141 1 0.0 2 1309 2 0.0 2 2310 1 0.0 2 1311 3 0.0 2 3
=== Adapter 2 ===
Sequence: AACMGGATTAGATACCCKG; Type: regular 5'; Length: 19; Trimmed: 1074437 times.
No. of allowed errors:0-5 bp: 0; 6-12 bp: 1; 13-19 bp: 2
Overview of removed sequenceslength count expect max.err error counts3 2 19959.9 0 27 1 78.0 1 0 18 2 19.5 1 1 110 6 1.2 1 4 211 1 0.3 1 112 3 0.1 1 2 113 5 0.0 1 3 214 24 0.0 2 17 715 51 0.0 2 32 14 516 71 0.0 2 56 12 317 134 0.0 2 92 30 1218 327 0.0 2 189 117 2119 1059175 0.0 2 1056863 2069 24320 13846 0.0 2 1817 10955 107421 744 0.0 2 5 10 72922 1 0.0 2 123 2 0.0 2 224 1 0.0 2 125 2 0.0 2 227 5 0.0 2 528 2 0.0 2 229 2 0.0 2 230 1 0.0 2 131 2 0.0 2 232 10 0.0 2 1049 1 0.0 2 151 1 0.0 2 1166 1 0.0 2 1291 6 0.0 2 6401 2 0.0 2 2409 1 0.0 2 1443 1 0.0 2 1460 2 0.0 2 2465 2 0.0 2 2
WARNING: One or more of your adapter sequences may be incomplete. Please see the detailed output above.

  • 蛋白质组学业务咨询
  • 扫码咨询相关问题
  • weinxin
  • 代谢组学业务咨询
  • 扫码咨询相关问题
  • weinxin

发表评论

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen: