BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,842,033	G→T	R67R (CGC→CGA)	hycG ←	hydrogenase activity

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,842,033	0	G	T	99.2%	334.1 / ‑5.6	117	R67R (CGC→CGA)	hycG	hydrogenase activity
Reads supporting (aligned to +/- strand): ref base G (0/1); new base T (51/65); total (51/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 117 total aligned reads displayed.

TGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAGCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCC  >  NC_000913/2841961‑2842105

tGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCAGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGt                                                                       >  1:721140/1‑76 (MQ=255)

  cACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTc                                                                     <  2:2930521/76‑1 (MQ=255)

  cACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTc                                                                     >  1:1932872/1‑76 (MQ=255)

  cACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTc                                                                     >  1:1639436/1‑76 (MQ=255)

   cgggggGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCt                                                                    <  2:3903952/74‑1 (MQ=255)

   aCGGGTGCCCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCt                                                                    <  1:2236158/76‑1 (MQ=255)

   aCGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCt                                                                    >  2:811118/1‑76 (MQ=255)

   aCGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCt                                                                    >  2:561547/1‑76 (MQ=255)

    cGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTg                                                                   >  1:4468909/1‑76 (MQ=255)

     gggTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGc                                                                  >  1:3052916/1‑76 (MQ=255)

      ggTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCa                                                                 <  1:2351904/76‑1 (MQ=255)

      ggTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCa                                                                 <  2:2489774/76‑1 (MQ=255)

       gTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCAt                                                                >  2:85594/1‑76 (MQ=255)

       gTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCAt                                                                <  1:2876489/76‑1 (MQ=255)

       gTGAACGCGCCGGTAAACAGTAAAATATCCGGATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCAt                                                                <  1:4302973/76‑1 (MQ=255)

         gACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCa                                                              <  2:2172323/76‑1 (MQ=255)

           ccGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCaaa                                                            >  2:180050/1‑76 (MQ=255)

            cgcgCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAAc                                                           >  1:1012871/1‑76 (MQ=255)

              cgcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAg                                                         >  1:3728328/1‑76 (MQ=255)

              cgcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAg                                                         <  2:559458/76‑1 (MQ=255)

              cgcCGGTAAACAGAAAAATATCCGCATGACGCGGTGAAGGAACGTCTTTAATGCCGAATCGGTCTGCATCAAACAg                                                         >  2:1728457/1‑76 (MQ=255)

                 cGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcgg                                                      >  1:1739977/1‑76 (MQ=255)

                 cGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcgg                                                      >  1:1642948/1‑76 (MQ=255)

                 cGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcgg                                                      <  1:918958/76‑1 (MQ=255)

                 cGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcgg                                                      >  1:3634112/1‑76 (MQ=255)

                  ggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggc                                                     <  1:3006104/76‑1 (MQ=255)

                  ggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggc                                                     <  2:2328077/76‑1 (MQ=255)

                  ggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggc                                                     <  1:2612660/76‑1 (MQ=255)

                  ggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggc                                                     >  1:4483647/1‑76 (MQ=255)

                    tAAGCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGa                                                   <  1:3526954/76‑1 (MQ=255)

                    tAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGa                                                   >  1:1344018/1‑76 (MQ=255)

                     aaaCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGaa                                                  <  1:2547771/76‑1 (MQ=255)

                      aaCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGaaa                                                 >  2:4202020/1‑76 (MQ=255)

                        cAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGc                                               <  2:579284/76‑1 (MQ=255)

                        cAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGc                                               <  1:1483688/76‑1 (MQ=255)

                          gTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCgg                                             >  2:213996/1‑75 (MQ=255)

                           tAAAATATCTGCATGTCGAGTTGTAGCAACGACTTTAATGCCGAATCGTTCTGCGACAATCTGCGGCGAAAGCgcg                                            <  1:3098875/76‑3 (MQ=255)

                           tAAAATATCCGCATGAGGCGGTGAAGGGACGACTTTAATGCCGAATCGTTCTGCATCAAGCAGCGGCGAAAGCGTg                                            <  2:1288593/76‑1 (MQ=255)

                           tAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTg                                            >  2:1078590/1‑76 (MQ=255)

                           tAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTg                                            <  1:4569337/76‑1 (MQ=255)

                             aaaTATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGc                                          >  1:3023808/1‑76 (MQ=255)

                             aaaTATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGc                                          <  1:3862130/76‑1 (MQ=255)

                              aaTATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCg                                         >  2:3309124/1‑76 (MQ=255)

                                tatCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGACAGCGTGGCgca                                       >  2:4254195/1‑74 (MQ=255)

                                 atCCGCTTGACGCGGTGAGGGAGCGGCTTTAATGCCGAATAGGTCTGCAGCAAGCAGCGGCGAAAGCGTGGCGaaa                                      <  2:4798718/76‑1 (MQ=255)

                                 atCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGaaa                                      >  1:1034016/1‑76 (MQ=255)

                                   ccGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAat                                    <  2:649913/76‑1 (MQ=255)

                                   ccGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAat                                    >  1:3734571/1‑76 (MQ=255)

                                    cGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatt                                   <  2:4209213/76‑1 (MQ=255)

                                    cGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatt                                   <  1:363325/76‑1 (MQ=255)

                                     gCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAattt                                  >  2:2958202/1‑76 (MQ=255)

                                     gCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAattt                                  <  1:1097522/76‑1 (MQ=255)

                                      cATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttc                                 <  1:1607153/76‑1 (MQ=255)

                                       aTGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcg                                >  2:4685189/1‑76 (MQ=255)

                                         gACGGGGGGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgat                              <  2:484229/76‑1 (MQ=255)

                                           cgcgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgattt                            >  2:3950815/1‑76 (MQ=255)

                                            gcgGTGAAGGAACGACTTTAATGCCGAGTCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttc                           <  2:4355432/76‑1 (MQ=255)

                                            gcgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAg                                                <  1:1492361/55‑1 (MQ=255)

                                            gcgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttc                           >  1:4290300/1‑76 (MQ=255)

                                             cgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttcg                          <  1:3585119/76‑1 (MQ=255)

                                              ggTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGc                         >  2:3219129/1‑76 (MQ=255)

                                              ggTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGGGGCGACAATTTCGATTTCGc                         >  2:927846/1‑76 (MQ=255)

                                               gTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCa                        <  1:4363140/76‑1 (MQ=255)

                                               gTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCa                        >  2:2262786/1‑76 (MQ=255)

                                               gTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCa                        <  1:4116149/76‑1 (MQ=255)

                                               gTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCa                        <  1:1055097/76‑1 (MQ=255)

                                               gTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCa                        >  1:878385/1‑76 (MQ=255)

                                                tGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCaa                       >  2:1263344/1‑76 (MQ=255)

                                                tGAAGGAACGACTTTAATGCCGAATCGTGCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATGTCGATTTGGCac                       >  2:2200750/1‑75 (MQ=255)

                                                 gAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAc                      <  1:2624710/76‑1 (MQ=255)

                                                  aaGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAcc                     <  1:4643723/76‑1 (MQ=255)

                                                  aaGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAcc                     <  2:3782526/76‑1 (MQ=255)

                                                   aGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCg                    <  2:1572308/76‑1 (MQ=255)

                                                   aGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCg                    <  2:3957052/76‑1 (MQ=255)

                                                   aGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCg                    <  2:4058060/76‑1 (MQ=255)

                                                    ggAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGt                   <  1:4010079/76‑1 (MQ=255)

                                                    ggAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGAGAACATTTCGATTTCGCAACCGt                   >  2:1850036/1‑76 (MQ=255)

                                                    ggAACGACTTTAAGGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGt                   <  2:1351037/76‑1 (MQ=255)

                                                     gAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGtt                  <  1:1005816/76‑1 (MQ=255)

                                                      aaCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTg                 <  1:1166504/76‑1 (MQ=255)

                                                       aCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGc                <  1:349739/76‑1 (MQ=255)

                                                       aCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGc                <  1:3303843/76‑1 (MQ=255)

                                                       aCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGc                <  1:3024135/76‑1 (MQ=255)

                                                        cggcTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCa               <  2:1473700/73‑1 (MQ=255)

                                                        cGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCa               >  1:2025070/1‑76 (MQ=255)

                                                         gACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAg              <  2:4399464/76‑1 (MQ=255)

                                                          aCTGTAATGCCGAATCGTTCTGCTTCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgc             <  2:4200806/76‑1 (MQ=255)

                                                           cTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgcc            >  2:769546/1‑76 (MQ=255)

                                                            tttAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgccg           <  2:3616964/76‑1 (MQ=255)

                                                             ttAATGCCGAATCGTTCTTCGTCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgccgc          <  1:229308/76‑1 (MQ=255)

                                                               aaTGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAACTTTCGATTTCGCAACCGTCGCTGccgccg        >  2:1442242/1‑76 (MQ=255)

                                                                aTGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCcgccgc       <  1:4686428/76‑1 (MQ=255)

                                                                aTGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCcgccgc       <  2:3762665/76‑1 (MQ=255)

                                                                 tGCCAAATCTTTCTCCACCAACCAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCa      <  1:2490313/76‑1 (MQ=37)

                                                                 tGCCAAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCACCCGTTGCAGCCGCCGCa      <  1:2684974/76‑1 (MQ=38)

                                                                   ccGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGCAACTTTCGCTTTCGCCACCGTTGCAGCCGGCGCAgc    >  2:1887484/1‑75 (MQ=255)

                                                                   ccGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGt    >  2:1862357/1‑76 (MQ=255)

                                                                   ccGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCACAACCGTTGCAGCCGCCGCAGt    <  1:3952855/76‑1 (MQ=255)

                                                                    cGAATCGTTATGCGTCAAACATCGGCGAAGGCGTGGCGAAAATTTCGATTTCTCAACCGTTGCAGCCGCCGCAGTc   <  1:2521374/76‑1 (MQ=37)

                                                                     gAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTcc  <  2:3413089/76‑1 (MQ=255)

TGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAGCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCC  >  NC_000913/2841961‑2842105

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