# Quality Check before metaphlan analysis **URL:** \ **Category:** MetaPhlAn\ **Created:** [June 2, 2020, 3:28am UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539 "2020-06-02T03:28:07Z")\ **Posts on this page:** 15\ **Page:** 1
**Author:** ![DEEPCHANDA7](https://avatars.discourse-cdn.com/v4/letter/d/eb9ed0/32.png) [@DEEPCHANDA7](https://forum.biobakery.org/u/DEEPCHANDA7)\ **Post date:** [June 2, 2020, 3:28am UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/1 "2020-06-02T03:28:07Z")
Hi friends!!! I want to analyse already submitted data in SRA/EMBL by MetaPhlAn3. Should I do quality check each time before I run with MetaPhlAn? Thanks & Regards, DC7 ---
**Author:** ![fbeghini](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/fbeghini/32/81_2.png) [@fbeghini](https://forum.biobakery.org/u/fbeghini)\ **Post date:** [June 4, 2020, 12:38pm UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/2 "2020-06-04T12:38:12Z")
Hi, assuming that the submitted data was already QC’d, there’s no need to to QC. Internally, MetaPhlAn discards reads longer less than 70bp ---
**Author:** ![DEEPCHANDA7](https://avatars.discourse-cdn.com/v4/letter/d/eb9ed0/32.png) [@DEEPCHANDA7](https://forum.biobakery.org/u/DEEPCHANDA7)\ **Post date:** [June 4, 2020, 12:50pm UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/3 "2020-06-04T12:50:25Z")
Thanks @fbeghini. Does MetPhlAn discards sequences containing bases with low quality score? DC7 ---
**Author:** ![fbeghini](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/fbeghini/32/81_2.png) [@fbeghini](https://forum.biobakery.org/u/fbeghini)\ **Post date:** [June 4, 2020, 1:02pm UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/4 "2020-06-04T13:02:33Z")
Yes, marker hits with MAPQ below 5 are discarded ---
**Author:** ![DEEPCHANDA7](https://avatars.discourse-cdn.com/v4/letter/d/eb9ed0/32.png) [@DEEPCHANDA7](https://forum.biobakery.org/u/DEEPCHANDA7)\ **Post date:** [June 5, 2020, 11:09am UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/5 "2020-06-05T11:09:49Z")
But, Sir @fbeghini , how can I be sure that those submitted data were QC’d ? ---
**Author:** ![fbeghini](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/fbeghini/32/81_2.png) [@fbeghini](https://forum.biobakery.org/u/fbeghini)\ **Post date:** [June 5, 2020, 2:37pm UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/6 "2020-06-05T14:37:24Z")
If nothing is specified on SRA, the trimming + filtering of low MAPQ alignments done by MetaPhlAn is sufficient ---
**Author:** ![DEEPCHANDA7](https://avatars.discourse-cdn.com/v4/letter/d/eb9ed0/32.png) [@DEEPCHANDA7](https://forum.biobakery.org/u/DEEPCHANDA7)\ **Post date:** [June 7, 2020, 3:39pm UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/7 "2020-06-07T15:39:23Z")
I have a WGS file where average read length is 35bp. So, are you saying that MetPhlAn will be unable to profile it. Right? (PRJEB23292) ---
**Author:** ![fbeghini](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/fbeghini/32/81_2.png) [@fbeghini](https://forum.biobakery.org/u/fbeghini)\ **Post date:** [June 8, 2020, 7:43am UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/8 "2020-06-08T07:43:02Z")
No, you can set the minimum read length using `--read_min_len 34` (`metaphlan -h`), in this case all reads with length below 34 will be discarded ---
**Author:** ![DEEPCHANDA7](https://avatars.discourse-cdn.com/v4/letter/d/eb9ed0/32.png) [@DEEPCHANDA7](https://forum.biobakery.org/u/DEEPCHANDA7)\ **Post date:** [June 11, 2020, 9:32pm UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/9 "2020-06-11T21:32:24Z")
Hi @fbeghini, if I set minimum read length to 34, in that case, do you think the profile output will be reliable? Will MetaPhlAn 3.0 be able to classify the reads perfectly? Thanks DC7 ---
**Author:** ![fbeghini](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/fbeghini/32/81_2.png) [@fbeghini](https://forum.biobakery.org/u/fbeghini)\ **Post date:** [June 12, 2020, 8:45am UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/10 "2020-06-12T08:45:45Z")
PRJEB23292 is a project which includes SRAs generated with the SOLiD platform. It is not a good idea converting colorspace reads into basespace, and BowTie2 does not support colorspace input. You should consider analyzing the dataset with a tool able to handle such input. ---
**Author:** ![DEEPCHANDA7](https://avatars.discourse-cdn.com/v4/letter/d/eb9ed0/32.png) [@DEEPCHANDA7](https://forum.biobakery.org/u/DEEPCHANDA7)\ **Post date:** [June 12, 2020, 1:00pm UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/11 "2020-06-12T13:00:30Z")
Sir, in connection with the sequencing platform, I have a query. Can I analyse WGS sequence data generated by 454 GS FLX platform (PRJNA32089) and Illumina Genome analyser II (PRJEB2054)? > NOTE: I am seeing the fastq files from PRJNA32089 are comparatively smaller (mostly less than 300Mb) and single-ended. Thanks and regards, DC7 ---
**Author:** ![fbeghini](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/fbeghini/32/81_2.png) [@fbeghini](https://forum.biobakery.org/u/fbeghini)\ **Post date:** [June 12, 2020, 2:10pm UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/12 "2020-06-12T14:10:54Z")
GAII reads are OK for MetaPhlAn, however, since MetaPhlAn is designed to work with short reads, on longer reads such as 454 reads you can try running the internal Bowtie2 using local alignment (`--bt2_ps very-sensitive-local`) instead of the default end-to-end ---
**Author:** ![DEEPCHANDA7](https://avatars.discourse-cdn.com/v4/letter/d/eb9ed0/32.png) [@DEEPCHANDA7](https://forum.biobakery.org/u/DEEPCHANDA7)\ **Post date:** [June 12, 2020, 7:29pm UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/13 "2020-06-12T19:29:03Z")
Sir, In that case should I set a `--min_alignment_len` value other than the default 100? I’m getting a message like this: `Warning! bt2_ps is set to local mode, and min_alignment_len is None, I automatically set min_alignment_len to 100! If you do not like, rerun the command and set min_alignment_len to a specific value.` ---
**Author:** ![DEEPCHANDA7](https://avatars.discourse-cdn.com/v4/letter/d/eb9ed0/32.png) [@DEEPCHANDA7](https://forum.biobakery.org/u/DEEPCHANDA7)\ **Post date:** [November 23, 2020, 10:11pm UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/14 "2020-11-23T22:11:29Z")
Sir @fbeghini - can you please suggest? ---
**Author:** ![fbeghini](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/fbeghini/32/81_2.png) [@fbeghini](https://forum.biobakery.org/u/fbeghini)\ **Post date:** [November 27, 2020, 10:52am UTC](https://forum.biobakery.org/t/quality-check-before-metaphlan-analysis/539/15 "2020-11-27T10:52:12Z")
I’d suggest to not lower the default value and keep it at 100