Select your language

GTでのコンセンサス配列の作り方

  • HKoike
  • HKoike's Avatar Topic Author
  • Visitor
  • Visitor

GTでのコンセンサス配列の作り方

25 Jul 2012 08:12
#
基本的な質問ですいません。マッピングした結果から、コンセンサス配列を作る際、contig.consなどで大文字のATGCと小文字のatgcは区別されているのだと思いますが、これはどのような違いでしょうか?qualityを反映しているのでしょうか。マッピングは、他のプログラムで出力されたBAM fileをインポートしています。どうぞ、よろしくお願いします。

Please Log in to join the conversation.

  • akr
  • akr's Avatar Topic Author
  • Visitor
  • Visitor

Re: Re: GTでのコンセンサス配列の作り方

26 Jul 2012 09:33
#
SAM/BAMファイルからのインポートの場合は、Soft Clipping領域の塩基を小文字(a,c,g,t)で表記しています。
ちなみに、GTのマッピング結果の場合は、Referenceに対してアラインメントできた領域の塩基を大文字(A,C,G,T)で、それ以外の領域の塩基を小文字で表記しています。

Please Log in to join the conversation.

  • HKoike
  • HKoike's Avatar Topic Author
  • Visitor
  • Visitor

Re: Re: GTでのコンセンサス配列の作り方

06 Aug 2012 01:52
#
コンセンサスが書かれた、contig.consでは、a,A,c,C,g,G,t,Gの他に、その他(塩基N)の数が出力されていますが、SAM/BAMファイルからのインポートの場合は、どのような場合にNと数えられるのでしょうか?
SAM/BAMで既にNとなっているのでしょうか? それともGTで変換などがあるのでしょうか? 五月雨式にいろいろと訊いて申し訳ありませんが、どうぞ、よろしくお願いします。

Please Log in to join the conversation.

  • akr
  • akr's Avatar Topic Author
  • Visitor
  • Visitor

Re: Re: GTでのコンセンサス配列の作り方

07 Aug 2012 14:52
#
SAM/BAMインポート時にコンセンサス配列を独自に決定しています。
決定方法は以下の2種類から選択できます。
1.多数決による方法=一番多い塩基をコンセンサス塩基とする。
2.最高Qualityによる方法
そのほかに、
Consensus Ratioというパラメータがあり、これで一番多い塩基の全体に占める必要最低割合を指定します。
もっとも多数の塩基がこの数値以上を占めるとその塩基がコールされますが、それ以下ではNとなります。

Please Log in to join the conversation.

  • saitoh
  • saitoh's Avatar Topic Author
  • Visitor
  • Visitor

Re: Re: GTでのコンセンサス配列の作り方

08 Aug 2012 06:17
#
GTのコンセンサスに関する資料を添付します。

Please Log in to join the conversation.

  • HKoike
  • HKoike's Avatar Topic Author
  • Visitor
  • Visitor

Re: Re: GTでのコンセンサス配列の作り方

09 Aug 2012 09:09
#
お返事、ありがとう御座いました。コンセンサスの配列を決定する以前の、contig.consのcolumn10、その他(塩基 N)の数の「N」も、コンセンサスを作る場合の多数決の結果が反映されているのでしょうか? contig.consのファイルを見た時に、column01 - column10まではコンセンサスの決定の基にするデータ(SAM/BAMから来たデータ)で、GTでのコンセンサス決定の結果が、column11以降に書かれているように思っていました。 それとも、コンセンサスの決定時の結果が、column10に書かれているのでしょうか?
勘違いがあったら申し訳ないのですが、お答え頂いたnは、column12のnについてのように感じたのですが。どうぞ、よろしくお願いします。

Please Log in to join the conversation.

Time to create page: 0.153 seconds
Powered by Kunena Forum
汽水音乐易歪歪HelloGPT易歪歪单机版沙盘纸飞机易翻译比特浏览器mumu模拟器