1 answer
I use cutadapt with a poly-G as adapter, you should allow some errors because the poly-G sometimes combine an occasional A-C-T base.
When I analyze paired-end, the second mate is sometimes a poly-G and they I remove it by testing if the read has more than 80% G's. If that's the case I disregard the entire read (or use it as single-end, depends on what I do with it later).
What I do is run fastqc and then test if poly-G is one of the over-represented sequences (and for what extent).
Then,actually, cutadapt with poly-G as adapter will remove the read but you should give it both mates as input (I think it will remove both of them but I'm not sure)
I checked with FastQC as you suggested, in my data set there is no over-represented sequence and mean quality score is 35. So I hope without any processing the data set I can directly run assembly. What do you think? I used Spades for assembly which also have some error correction steps in ion-hammer.
Log in to answer this question.
-
Bash script error on reading coordinates in a file
written by Ali HEBRA •
-
dbSNP IDs in the online database but not in FTP downloads
written by agd27 •
-
Query Regarding 454 paired end Metagenomic Data
written by subodh •
-
Association analysis with eigengenes
written by iside •
-
How to run Spades For Nextseq data
written by jeccy.J
-
How To Buildup A Genome Assembly Workflow With Simple Scripts???
written by HG
-
How To Increase Contig Length For Gap Filling ?
written by HG
-
How To Run Images??
written by HG
-
Amos Package Instalation Error
written by HG
-
Iterative Mapping And Assembly For Gap Elimination (Images)
written by HG
