BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,766,423	G→A	R231R (CGG→CGA)	b2632 →	putative GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,766,423	0	G	A	100.0%	350.6 / NA	119	R231R (CGG→CGA)	b2632	putative GTP‑binding protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (60/59); total (60/59)

Only 100 of 119 total aligned reads displayed.

AGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGT  >  NC_000913/2766351‑2766491

gaagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGGAGTAGCGCTACTGCAACATCCCTTTCGAACa                                                                   <  2:5683627/74‑1 (MQ=25)

agagCGGATGATTAAGTGTCTGCCGCGTGAGGCCCGCAGCCCGGTAGTAGCGCTACTGCCACCTCCCTTTCGAACa                                                                   >  2:209357/1‑76 (MQ=21)

agagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGCAGCGCTACTGCAACATCCCTTTCGAACa                                                                   >  2:4547804/1‑76 (MQ=37)

 gagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaa                                                                  >  2:1670509/1‑75 (MQ=37)

  agCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaac                                                                 <  1:4077165/76‑3 (MQ=37)

  agCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaac                                                                 <  2:2586285/76‑3 (MQ=37)

   gCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaacg                                                                <  2:5390128/76‑4 (MQ=37)

   gCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaacg                                                                <  2:5059141/76‑4 (MQ=25)

   gCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaacg                                                                <  1:4638339/76‑4 (MQ=37)

   gCGGATGATTAAGTGACAGCCGCTGGAGGCCAGCAGACCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaacg                                                                <  2:5104454/76‑4 (MQ=255)

       atgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAg                                                            >  1:3165541/1‑76 (MQ=255)

       atgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAg                                                            >  1:1532754/1‑76 (MQ=255)

       atgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAg                                                            >  2:2039758/1‑76 (MQ=255)

        tgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGc                                                           >  2:2388854/1‑76 (MQ=255)

        tgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGc                                                           <  2:1262962/76‑1 (MQ=255)

        tgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGc                                                           >  1:146607/1‑76 (MQ=255)

         gatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGcc                                                          >  1:4009649/1‑76 (MQ=255)

           ttAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCg                                                        >  2:1933245/1‑76 (MQ=255)

           ttAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCg                                                        >  1:3746567/1‑76 (MQ=255)

              aGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGa                                                     <  1:1211456/76‑1 (MQ=255)

              aGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGa                                                     <  2:2538494/76‑1 (MQ=255)

              aGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGa                                                     <  2:4808919/76‑1 (MQ=255)

               gtgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAg                                                    <  1:5553148/76‑1 (MQ=255)

               gtgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAg                                                    >  2:3001934/1‑76 (MQ=255)

                tgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGc                                                   <  1:1316123/76‑1 (MQ=255)

                tgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGc                                                   <  2:4591180/76‑1 (MQ=255)

                 gtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCa                                                  >  2:3991811/1‑76 (MQ=255)

                   cTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAgg                                                >  2:2973290/1‑76 (MQ=255)

                   cTGCCGCGTCACGCCAGCAGCCCGGCACTAGCGCAACTCCAACATCCCTTTCGAACAACGGAAGCCTGGcagcagc                                                >  2:630378/1‑75 (MQ=255)

                    tGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGc                                               <  2:1354789/76‑1 (MQ=255)

                    tGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGc                                               >  1:818832/1‑76 (MQ=255)

                     gCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCa                                              <  2:836/76‑1 (MQ=255)

                      ccGCGTGAGGCCATCAGCCCGGTAGTAGCTCTACTGCAACATCCCTTTCGAACAACGGTAGCCCTGGAGCAGGCAc                                             <  1:785791/76‑1 (MQ=255)

                      ccGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACCACGGTAGCCCGGGAgcacgcac                                             >  2:1629627/1‑76 (MQ=255)

                      ccGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCAc                                             >  2:1750000/1‑76 (MQ=255)

                      ccGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCAc                                             >  1:2367311/1‑76 (MQ=255)

                      ccGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCAc                                             <  2:5660997/76‑1 (MQ=255)

                      ccGCGGGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCAc                                             <  2:276349/76‑1 (MQ=255)

                       cgcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                            <  1:4573746/76‑1 (MQ=255)

                        gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                           >  1:4599904/1‑76 (MQ=255)

                        gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                           <  1:210534/76‑1 (MQ=255)

                        gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                           >  2:2459856/1‑76 (MQ=255)

                        gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                           >  1:430958/1‑76 (MQ=255)

                          gTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGa                                         >  1:3719836/1‑76 (MQ=255)

                          gTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCCGGAGCAGGCACGTGa                                         >  2:1042085/1‑76 (MQ=255)

                          gTGAGGCCAGCAGCCCGGTAGGAGCGCTACTGCAACATCCCTGTCGAACAACGGTAGCCCAAGAGCAGGCACGTGa                                         >  2:3036058/1‑76 (MQ=255)

                             aggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGggagga                                      >  1:25500/4‑76 (MQ=255)

                              ggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACGACGGTAGCCCGGGAGCAGGCACGTGACGAt                                     <  1:3566203/74‑1 (MQ=255)

                              ggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAt                                     >  2:2207369/3‑76 (MQ=255)

                               gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAtt                                    <  1:1207221/75‑1 (MQ=255)

                               gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAtt                                    <  2:2369433/75‑1 (MQ=255)

                               gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAtt                                    <  1:2207369/75‑1 (MQ=255)

                               gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAtt                                    <  2:3520868/75‑1 (MQ=255)

                               gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAtt                                    <  2:964118/75‑1 (MQ=255)

                               gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAtt                                    >  1:949462/2‑76 (MQ=255)

                               gacagcagCCCGGGAGGAGCGCTACGGCAACATCCCGTTCGAACAACGGGAGCCCGGGAGCAGGCACGTGACGAtt                                    <  2:3539888/74‑1 (MQ=255)

                                ccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAttt                                   >  2:4237604/1‑76 (MQ=255)

                                 cagcagCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGCCGATTTc                                  >  2:4322462/1‑76 (MQ=255)

                                 cagcagCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTc                                  <  1:4411552/76‑1 (MQ=255)

                                 cagcagCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTc                                  >  1:530170/1‑76 (MQ=255)

                                  agcagcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCg                                 >  1:722981/1‑76 (MQ=255)

                                  agcagcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCg                                 >  1:3263094/1‑76 (MQ=255)

                                   gcagcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGCTTTCgg                                >  1:5654019/1‑76 (MQ=255)

                                   gcagcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCgg                                <  1:3702168/76‑1 (MQ=255)

                                    cagcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGt                               <  1:5749668/76‑1 (MQ=255)

                                     agcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCCGGCACGTGACGATGTCGGTg                              >  2:2164792/1‑76 (MQ=255)

                                     agcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTg                              <  1:4225689/76‑1 (MQ=255)

                                         cGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGCTTTCGGTGAAAc                          >  2:524877/1‑76 (MQ=255)

                                          ggTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAcc                         <  2:3648871/76‑1 (MQ=255)

                                           gtagtaGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCg                        <  1:5329281/76‑1 (MQ=255)

                                            tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGt                       <  1:1986809/76‑1 (MQ=255)

                                              gtagCGCTACTGCAACATCCCTTTCGAACAGCGGGAGCCCGGGAGCAGGCACGTGACGAGTTCGGTGAAAACGTAg                     >  1:4766858/1‑76 (MQ=255)

                                               tagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAgg                    <  1:753267/76‑2 (MQ=255)

                                               tagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAgg                    >  2:1384827/1‑75 (MQ=255)

                                               tagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGCAgg                    >  2:488381/1‑75 (MQ=255)

                                                 gcgcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggtg                  <  2:265493/76‑4 (MQ=255)

                                                 gcgcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGGAggtg                  >  1:357846/1‑73 (MQ=255)

                                                     tACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCAt              >  2:1266174/1‑76 (MQ=255)

                                                     tACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCAt              >  1:4794514/1‑76 (MQ=255)

                                                       cTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc            >  1:1294711/1‑76 (MQ=255)

                                                        tGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACt           <  1:1750000/76‑1 (MQ=255)

                                                        tGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACt           <  2:4859376/76‑1 (MQ=255)

                                                        tGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACt           >  1:1342234/1‑76 (MQ=255)

                                                         gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc          >  1:3841074/1‑76 (MQ=255)

                                                         gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc          >  1:363675/1‑76 (MQ=255)

                                                         gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc          <  2:2675459/76‑1 (MQ=255)

                                                         gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGGGAAACCGGAGGTGCCATCCTc          >  2:438540/1‑76 (MQ=255)

                                                          caccATCCCTTTCGATCACCGGTAGCCCGGGAGCAGGCACGTGCCGATTTCGGTGAATCCGTAGGTGCCATACTCg         <  1:2133798/73‑1 (MQ=255)

                                                          cAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCg         >  1:3063847/1‑76 (MQ=255)

                                                          cAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCg         >  1:522012/1‑76 (MQ=255)

                                                           aaCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGa        >  1:1215049/1‑76 (MQ=255)

                                                            aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGAt       <  2:5437396/76‑1 (MQ=255)

                                                            aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGCAGGCGCCATACTCGAt       >  2:761070/1‑76 (MQ=255)

                                                             cATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa      <  2:3063847/76‑1 (MQ=255)

                                                             cATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa      <  1:3001934/76‑1 (MQ=255)

                                                               tCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg    >  2:4660875/1‑76 (MQ=255)

                                                               tCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg    >  2:3723999/1‑76 (MQ=255)

                                                               tCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg    <  1:640122/76‑1 (MQ=255)

                                                               tCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg    <  2:785076/76‑1 (MQ=255)

                                                                cccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACgg   <  1:1544014/76‑1 (MQ=255)

AGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGT  >  NC_000913/2766351‑2766491

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 7 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`