This is a test version of Biostars. For the public version, visit https://www.biostars.org.
deal with the 16S rDNA sequence from ion torrent

Now I got some Fastq file from the ion torrent machine. And the person who give me this Fastq file tell me that the fastq doesn't have barcode , they are the v3 and v6 sequence from the fecal microorganism 16S pcr. I have to distinguish between V3 and V6, and use this file to try to use the PICRUSt(in QIIME).

I am totally a freshman in this. I have just tried the QIIME 454 tutorial, and Know the Mothur was used to deal with 16S sequence. Now, I am try to asking what should I do in this. Tell me whats the "begin"

genome 16s

2873@SA1.129: 1.0
7772@SA1.132: 1.0
1602@SA1.131: 1.0
0268@HSANC.3: 1.0
54162@KIDC.225: 1.0
4089@HSAA2.103: 1.0
4061@KIDC.217: 1.0
8596@KIDC.226: 1.0
7016@HSAA2.140: 1.0
0640@HSANC.4: 1.0
39528@HSAA2.142: 1.0
078@HSANC.2: 1.0
6358@HSANC.1: 1.0
7655@KIDC.218: 1.0
24878@HSAA2.141: 1.0
31890@SA1.135: 1.0
0110@KIDC.216: 1.0

HSANC.1: 7076.0
HSANC.4: 9315.0
HSANC.3: 10638.0
HSANC.2: 30266.0
KIDC.226: 40108.0
KIDC.225: 47014.0
KIDC.218: 64059.0
KIDC.217: 68594.0
SA1.135: 71600.0
SA1.132: 72871.0
HSAA2.143: 74088.0
SA1.131: 77770.0
KIDC.216: 86356.0
HSAA2.142: 87653.0
HSAA2.140: 124876.0
HSAA2.103: 131888.0
HSAA2.141: 139526.0
SA1.129: 154160.0

The result from the ion torrent machine are some fastq file. In the first line they have same title @OJWK, so I use the bash commands sed and awk) to change their first line to

@A143_#Number#__CTCCTACGGGAGGCAGCAG_Adult_2_ATTACCGCGGCTGCTGG_20160131_Adult_2_chengjun_barcode=2
______^^^^^^^^

and

@A143_1__CTCCTACGGGAGGCAGCAG_Adult_2_ATTACCGCGGCTGCTGG_20160131_Adult_2_chengjun_barcode=2
______^

and use bash command cat make a merge.fastq, and use Qiime script(Convert) let the fastq file become two file== fna and qual. It means I split library by myself, but I dont filter the fastq, means they would have many low level sequence...

my question:

  1. it work? "split library by myself(use bash command)
  2. how to filter the low level sequence? or if I don't filter the low level sequence what would happen?
  3. anybody analysis ion torrent data? I want some tutorial or example~~~

0 answers

No answers yet.

Log in to answer this question.