Hi!
I’m trying to make contigs from my paired end MiSeq reads following your MiSeq SOP. I examined the output of make.contigs in stability.trim.contigs.fasta file and I found out that the contigs still contain sequences for my primers, barcodes and adaptors. I should get rid of these extras, right? Next I tried joining the contigs and using an .oligos file with my primers and bar code sequences (as suggested in the SOP), but this resulted in an empty join.trim.contigs.fasta file. Obviously I’m doing something wrong. I’m not sure if I compiled my .oligos file correctly. Any suggestions?
Thanks!
Here is an example of a contig that contains my barcodes and primers.
M01168_43_000000000-A3F87_1_1101_19338_7496
GTGCACCTCCCCCCCCCTCACATTATTTGCCTTTTTTTGTTTTAATGATACGGCGCCCACCGAGATCTACACCATGACGTACGTCGTGCACCAGGCGCGAATCATTCACAATGGGCGCAAGCCTGATGATGCCACGCCGCGTGCGGGATGAAGGCCCTCTGGGTTGTAAACCGCTGTCAATGAAGGGGAAACGTATGGTGTTAACAGCGCCATATTTGACTGCTTCGTAAAGGAAGCCCCGGCAAAACATGTGCCAGCCGCCGCGGTAAAACCAGCCACTGACTGACTGGCTCCATATCGCATGCCGTCTTCTGCTTGGAAAAAAAAAAAAACAGCCAGGCACATAGACACAT
And here is a sample of my .oligos file.
primer GYGCASCAGKCGMGAAW GCBGGTDTTACCGCGGCGGCTGRCA
BARCODE CATGACGTACGTC GAGTCA UH-1
BARCODE CATGACGTACGTC CTATGA UH-2
Could you post the forward and reverse reads? What is the trash code mothur is giving you in the scrap file for this sequence?
Here are some sequences that are hopefully a better illustration of this problem. The first 2 sequences are individual fwd and rev reads and the 3rd one is a contig that I got by using make.contigs(file=stability.file). The primer sequences (below) can be found in the individual reads + in the contig. The way the contig is constructed doesn’t seem right to me (the amplicons I’m trying to analyze should be about 200 bp long). What can I do to get rid of the primer sequences from the reads?
@M01168:43:000000000-A3F87:1:1101:14789:1828 1:N:0:13
CTCACAGCCCTGGTCATGGCCCGAACCTTGTGGCGATTGCGGTTTCGGGGATTGCGGCACGTGTCAGCCGCCGCGGTAAAACCAGCCACTGACTGACAGGCGCCATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAATAAAACAAAGCAATAAGCGACTCTAAACAGACAGGATGACATATCAACCATTAAAAAAGTGACAAGAAAAAGGACAAAAGTTGCTAAAAAAAAAGAACGACAAAGCTGTGG
@M01168:43:000000000-A3F87:1:1101:14789:1828 2:N:0:13
CGTGCCGCAATACCCGAAACCGCAAGCGCCACAAGGTTCGGGCCATGACCAGGGCTGTGATATTCGCGCCTGGGGCGCGACGTACGTCATGGTGTAGATCTCGGTGGTCGCCGTATCATTAAAAATAAAAAAACTTCACAACTACCGACTGCATTATTAAAATAGCTCCAATGAACCGATAGATAAACGAAAAGAATAGAACATGACGAGATAGGGGAGACTGGGAGACGAAAGAAGGGCGGAAGAGGAAA
M01168_43_000000000-A3F87_1_1101_14789_1828
TTTCCTCTTCCGCCCTTCTTTCGTCTCCCAGTCTCCCCTATCTCGTCATGTTCTATTCTTTTCGTTTATCTATCGGTTCATTGGAGCTATTTTAATAATGCAGTCGGTAGTTGTGAAGTTTTTTTATTTTTAATGATACGGCGACCACCGAGATCTACACCATGACGTACGTCGCGCCCCAGGCGCGAATATCACAGCCCTGGTCATGGCCCGAACCTTGTGGCGNTTGCGGTTTCGGGNATTGCGGCACGTGTCAGCCGCCGCGGTAAAACCAGCCACTGACTGACAGGCGCCATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAATAAAACAAAGCAATAAGCGACTCTAAACAGACAGGATGACATATCAACCATTAAAAAAGTGACAAGAAAAAGGACAAAAGTTGCTAAAAAAAAAGAACGACAAAGCTGTGG
primer GYGCASCAGKCGMGAAW GCBGGTDTTACCGCGGCGGCTGRCA
BARCODE CATGACGTACGTC GGAGCC EH-1
Sorry, I forgot to mention earlier that the scrap file error code in the above cases was bf.