Has any one created a MetaSIM error profile file (.mconf) for Illumina reads >100bp?
I am considering using MetaSIM to create a metagenomic synthetic reads dataset. I want to use an Illumina error profile. I see that MetaSIM includes Illumina error profile files for 36bp and 80bp read lengths. In their manual they say you can create your own error profile, the .mocnf file. I need an error profile for reads that are longer than 100bp.
Does anyone have experience with creating an error profile file for MetaSIM that covers longer reads?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Can ART reads generator create NextSeq 550 synthetic reads?
written by Tawny
-
K-mer selection for correcting long reads with short reads
written by zack.saud •
-
How to extract sequences from fasta that have a certain length ?
written by v.berriosfarias •
-
Expanding gene family using hmmer
written by el97004 •
-
Efficient read mapping and alignment to reference set
written by endre.sebestyen •
-
BBMap randomreads.sh exact path to ref error: Could not find or load main class
written by Tawny
-
Running HMMER via Python
written by Ghoti •
-
art_illumina use custom error profile to simulate reads
written by jjmmii •
-
How to calculate organism relative abundance from Illumina metagenomic data?
written by nkuyfq
-
Why Longer Reads Must Be Trimmed Or Divided Into 36Bp ?
written by Liuyunlong •
Hello Tawny!
I have no the same problem as you had several months ago. I would like to create a MetaSim empirical error model file (.mconf) for Ilumina reads 100 bp. Did you solve this problem? What did you do?
Thank you so much
Beatriz
I ended up using the error profile that is provided with MetaSIM and just expanding it. So I added new rows with similar values until I reached the read length I needed for my project. Then I plotted my revised profile file so I could look at it and make sure it looked like what I would expect for an Illumina run of the same read length. My plot had a few unusual dips so I went back and tweaked the values I added to the error model until I got a plot that reflected what I have seen in real data.
Thank you so much Tawny. I see what you mean. I was thinking in doing something similar but I was not sure if it was right. In the FAQ section of MetaSim there is a light explanation on this aspect, but not enough to understand it.
Thank you so much
Bea