GTでのコンセンサス配列の作り方
- HKoike
- Topic Author
- Visitor
-
24 7月 2012 23:12 #512
by HKoike
GTでのコンセンサス配列の作り方 was created by HKoike
基本的な質問ですいません。マッピングした結果から、コンセンサス配列を作る際、contig.consなどで大文字のATGCと小文字のatgcは区別されているのだと思いますが、これはどのような違いでしょうか?qualityを反映しているのでしょうか。マッピングは、他のプログラムで出力されたBAM fileをインポートしています。どうぞ、よろしくお願いします。
Please ログイン to join the conversation.
- akr
- Topic Author
- Visitor
-
26 7月 2012 00:33 #513
by akr
Replied by akr on topic Re: GTでのコンセンサス配列の作り方
SAM/BAMファイルからのインポートの場合は、Soft Clipping領域の塩基を小文字(a,c,g,t)で表記しています。
ちなみに、GTのマッピング結果の場合は、Referenceに対してアラインメントできた領域の塩基を大文字(A,C,G,T)で、それ以外の領域の塩基を小文字で表記しています。
ちなみに、GTのマッピング結果の場合は、Referenceに対してアラインメントできた領域の塩基を大文字(A,C,G,T)で、それ以外の領域の塩基を小文字で表記しています。
Please ログイン to join the conversation.
- HKoike
- Topic Author
- Visitor
-
05 8月 2012 16:52 #514
by HKoike
Replied by HKoike on topic Re: GTでのコンセンサス配列の作り方
コンセンサスが書かれた、contig.consでは、a,A,c,C,g,G,t,Gの他に、その他(塩基N)の数が出力されていますが、SAM/BAMファイルからのインポートの場合は、どのような場合にNと数えられるのでしょうか?
SAM/BAMで既にNとなっているのでしょうか? それともGTで変換などがあるのでしょうか? 五月雨式にいろいろと訊いて申し訳ありませんが、どうぞ、よろしくお願いします。
SAM/BAMで既にNとなっているのでしょうか? それともGTで変換などがあるのでしょうか? 五月雨式にいろいろと訊いて申し訳ありませんが、どうぞ、よろしくお願いします。
Please ログイン to join the conversation.
- akr
- Topic Author
- Visitor
-
07 8月 2012 05:52 #515
by akr
Replied by akr on topic Re: GTでのコンセンサス配列の作り方
SAM/BAMインポート時にコンセンサス配列を独自に決定しています。
決定方法は以下の2種類から選択できます。
1.多数決による方法=一番多い塩基をコンセンサス塩基とする。
2.最高Qualityによる方法
そのほかに、
Consensus Ratioというパラメータがあり、これで一番多い塩基の全体に占める必要最低割合を指定します。
もっとも多数の塩基がこの数値以上を占めるとその塩基がコールされますが、それ以下ではNとなります。
決定方法は以下の2種類から選択できます。
1.多数決による方法=一番多い塩基をコンセンサス塩基とする。
2.最高Qualityによる方法
そのほかに、
Consensus Ratioというパラメータがあり、これで一番多い塩基の全体に占める必要最低割合を指定します。
もっとも多数の塩基がこの数値以上を占めるとその塩基がコールされますが、それ以下ではNとなります。
Please ログイン to join the conversation.
- saitoh
- Topic Author
- Visitor
-
07 8月 2012 21:17 #516
by saitoh
Replied by saitoh on topic Re: GTでのコンセンサス配列の作り方
GTのコンセンサスに関する資料を添付します。
Please ログイン to join the conversation.
- HKoike
- Topic Author
- Visitor
-
09 8月 2012 00:09 #517
by HKoike
Replied by HKoike on topic Re: GTでのコンセンサス配列の作り方
お返事、ありがとう御座いました。コンセンサスの配列を決定する以前の、contig.consのcolumn10、その他(塩基 N)の数の「N」も、コンセンサスを作る場合の多数決の結果が反映されているのでしょうか? contig.consのファイルを見た時に、column01 - column10まではコンセンサスの決定の基にするデータ(SAM/BAMから来たデータ)で、GTでのコンセンサス決定の結果が、column11以降に書かれているように思っていました。 それとも、コンセンサスの決定時の結果が、column10に書かれているのでしょうか?
勘違いがあったら申し訳ないのですが、お答え頂いたnは、column12のnについてのように感じたのですが。どうぞ、よろしくお願いします。
勘違いがあったら申し訳ないのですが、お答え頂いたnは、column12のnについてのように感じたのですが。どうぞ、よろしくお願いします。
Please ログイン to join the conversation.
Time to create page: 0.090 seconds