how ro see adapter contamination in Illumina reads

[paa6]

I have illumina read file..which is bacterial DNA sequence...I have used geneious software to assembly it, while assembly I have found that there was vector contamination and it was removed by software since I have given trimming option and I got 1,610 contigs.

but now I am performing the same assembly by using velvet. I have my fastqc report and according to that report sequence duplication level is bad, overrepresented sequences and kmer content showing warning. (I have attached these three files) So, I reached to conclusion that I have adapter contamination on the basis of the sequence I have got in overrepresented sequences. I have seen that GATCGGAAGAGC is adapter contamination because I have seen it in adapter files provided to custmoer given by illumina technology.

Problem is my PI asked me to find that adaptor contamination sequence in my reads, which I was not able to So, he asked me que. that why can't u find it?? I am new to de novo assembly, I dont know what am I supposed to answer and he gave me 1 hrs. to find it. Please help!!!

[relipmoc]

I have illumina read file..which is bacterial DNA sequence...I have used geneious software to assembly it, while assembly I have found that there was vector contamination and it was removed by software since I have given trimming option and I got 1,610 contigs.

but now I am performing the same assembly by using velvet. I have my fastqc report and according to that report sequence duplication level is bad, overrepresented sequences and kmer content showing warning. (I have attached these three files) So, I reached to conclusion that I have adapter contamination on the basis of the sequence I have got in overrepresented sequences. I have seen that GATCGGAAGAGC is adapter contamination because I have seen it in adapter files provided to custmoer given by illumina technology.

Problem is my PI asked me to find that adaptor contamination sequence in my reads, which I was not able to So, he asked me que. that why can't u find it?? I am new to de novo assembly, I dont know what am I supposed to answer and he gave me 1 hrs. to find it. Please help!!!

try

$ grep -c 'GATCGGAAGAGC' reads.fastq
$ grep -c reads.fastq | awk '{print $1/4}'

then you will get an estimation of the contaminant ratio.

for adapter trimming, I suggest using skewer. For your case, you don't need to specify the adapter sequence since it's the same as the default TruSeq3 adapter sequence.

Good luck!

[paa6]

try

$ grep -c 'GATCGGAAGAGC' reads.fastq
$ grep -c reads.fastq | awk '{print $1/4}'

then you will get an estimation of the contaminant ratio.

for adapter trimming, I suggest using skewer. For your case, you don't need to specify the adapter sequence since it's the same as the default TruSeq3 adapter sequence.

Good luck!

THanks for the quick reply!! I have typed $ grep -c 'GATCGGAAGAGC' reads.fastq
and I got 28875..what is this mean??
also I am doing SE assembly while skewer is for PE...

[yueluo]

You can type grep --help for a brief description of OPTIONS for grep.

-c, --count only print a count of matching lines per FILE

The result you got was 28875, suggesting that 28875 reads contained the substring of 'GATCGGAAGAGC' - which is most likely adapter contamination.

[paa6]

you can type grep --help for a brief description of options for grep.



The result you got was 28875, suggesting that 28875 reads contained the substring of 'gatcggaagagc' - which is most likely adapter contamination.

ohh ok thanks!!!