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%	396.7 / NA	129	R67R (CGC→CGA)	hycG	hydrogenase activity
Reads supporting (aligned to +/- strand): ref base G (0/0); major base T (65/63); minor base C (0/1); total (65/64)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

TGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAGCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGT  >  NC_000913/2841961‑2842103

tGCGCGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGt                                                                     <  2:528919/76‑1 (MQ=255)

tGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGt                                                                     <  1:5476578/76‑1 (MQ=255)

 gCACGGGGGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGtt                                                                    <  2:1791877/76‑1 (MQ=255)

  gaggggTGACCGGGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTc                                                                   <  2:1631588/73‑1 (MQ=255)

   aCGGGGGCCCGCGCCGGGAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCt                                                                  <  2:5526469/76‑1 (MQ=255)

    cGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTg                                                                 <  1:3473020/76‑1 (MQ=255)

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

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

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

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

        tGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATc                                                             >  1:5026999/1‑76 (MQ=255)

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

             gcgcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACa                                                        <  2:2755105/76‑1 (MQ=255)

             gcgcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACa                                                        >  1:4374621/1‑76 (MQ=255)

             gcgcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACa                                                        >  2:3573506/1‑76 (MQ=255)

               gcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCCAACAgc                                                      >  2:5094670/1‑76 (MQ=255)

               gcCGGTAAACAGTAAAATATCCGCATGACGAGGTGAAGGAACGACTTTAATGCCGAATCGTACTCCCTCACACAgc                                                      >  2:765270/1‑76 (MQ=255)

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

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

                 cGGTAAAAACTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcga                                                    <  2:239592/76‑2 (MQ=255)

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

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

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

                   gTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggcg                                                  <  1:3712072/76‑1 (MQ=255)

                   gTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggcg                                                  >  2:5588334/1‑76 (MQ=255)

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

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

                      accAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGaaa                                               <  2:1089881/74‑1 (MQ=255)

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

                      aaCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGaaa                                               <  2:757181/76‑1 (MQ=255)

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

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

                        cAGTCAAATATCCGCATGACGCGGTGAAGGACCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGc                                             <  1:3024860/76‑1 (MQ=255)

                         aGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCg                                            <  1:1349485/76‑1 (MQ=255)

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

                          gTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGt                                           >  1:3775199/1‑76 (MQ=255)

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

                             agatatCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGc                                        <  2:2294858/74‑1 (MQ=255)

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

                               atatCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGa                                      >  1:3382734/1‑76 (MQ=255)

                                tatCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGaa                                     <  1:1706548/76‑1 (MQ=255)

                                tatCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGaa                                     >  1:4248504/1‑76 (MQ=255)

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

                                  tcggCATGACGCGGTGAAGGCACGACTTTAATGCCGACTCGTTCTGCATCAAACAGCGGCGAAAGCGTGTCGaaaa                                   <  1:2353693/73‑1 (MQ=255)

                                   ccGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGCAAat                                  >  1:1543856/1‑76 (MQ=255)

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

                                    cGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatt                                 >  2:2072341/1‑76 (MQ=255)

                                    cGCATGACACGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatt                                 >  1:2726770/1‑76 (MQ=255)

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

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

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

                                      cATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttc                               >  2:4108326/1‑76 (MQ=255)

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

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

                                        tGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcga                             >  1:3232331/1‑76 (MQ=255)

                                        tGACGCGGGGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcga                             >  2:2839850/1‑76 (MQ=255)

                                         gACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgat                            <  1:4589401/76‑1 (MQ=255)

                                         gACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgat                            >  1:5296019/1‑76 (MQ=255)

                                          aCGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgctt                           >  2:4543737/1‑76 (MQ=255)

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

                                            gcgGTGAAGGCACGCCTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCTAAAGCGTGGCGAAAatttcgatttc                         <  1:2026765/76‑1 (MQ=255)

                                            gcgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttc                         <  2:2654981/76‑1 (MQ=255)

                                            gcgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttc                         <  1:5590239/76‑1 (MQ=255)

                                             gggTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttcg                        <  2:5098156/75‑1 (MQ=255)

                                             cgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCGTCAAACAGCGGCGAAAGCGTGGCGAAAatttcgctttcg                        >  2:5144034/1‑76 (MQ=255)

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

                                             cgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttcg                        >  1:4221469/1‑76 (MQ=255)

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

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

                                              ggTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGc                       <  2:2920739/76‑1 (MQ=255)

                                              ggTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGc                       <  2:1092527/76‑1 (MQ=255)

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

                                                  aaGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAcc                   >  1:3391182/1‑76 (MQ=255)

                                                    ggATCGACTTTAAGGCCGAGTCGTTCTGCACCACACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAGCCGt                 <  2:853551/76‑1 (MQ=255)

                                                    ggAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGt                 >  2:754629/1‑76 (MQ=255)

                                                    ggAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGt                 >  1:3946882/1‑76 (MQ=255)

                                                     gAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGtt                >  1:1786397/1‑76 (MQ=255)

                                                     gAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAAGTTCGATTTCGCAACCGtt                <  2:2076114/76‑1 (MQ=255)

                                                      aaCGACTTTAATGCCGAATCGTTCTGCATCAAACCGCGGCGACAGAGTGGCGAAAATTTCGATTTCGCAACCGTTg               >  2:1872310/1‑76 (MQ=255)

                                                      aaCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTg               >  1:1094609/1‑76 (MQ=255)

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

                                                       aCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGc              >  2:2046981/1‑76 (MQ=255)

                                                        cGACTTTCATGCCGAATCGCGCTGCATCAAACACCGGCGCAAGCGTCGCGAAAACTACGACATCGCAAACGCTGCa             >  2:2232486/1‑76 (MQ=255)

                                                        cGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCa             <  2:5296018/76‑1 (MQ=255)

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

                                                         gACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAg            >  2:1416707/1‑76 (MQ=255)

                                                          aCTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgc           <  2:841180/76‑1 (MQ=255)

                                                           cTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgcc          <  1:1815393/76‑1 (MQ=255)

                                                           cTTTAATGCCGAATCGTTCTGCATAAAACAGCGGCGAAAGCGTGGCGAAAATTTCGTTTTCGCAACCGTTTCTgcc          >  2:3030241/1‑76 (MQ=255)

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

                                                            tttAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTAGCGAAAATTTCGATTTCGCAACCGTTGCAgccg         <  1:4108327/76‑1 (MQ=255)

                                                              tAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTCCAgccggc       >  2:570442/1‑74 (MQ=255)

                                                               aaTGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCCACCGTTGCAgccgccc      >  2:1019915/1‑75 (MQ=255)

                                                               aaTGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGccgccg      <  2:4550001/76‑1 (MQ=255)

                                                                atcccGAACCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCcgccgc     <  1:2142973/73‑1 (MQ=255)

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

                                                                   ccGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAgc  >  2:4301863/1‑75 (MQ=255)

                                                                   ccGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGt  <  1:4442369/76‑1 (MQ=255)

                                                                   ccGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGt  <  2:3232331/76‑1 (MQ=255)

                                                                   ccGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGt  <  2:4846104/76‑1 (MQ=255)

TGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAGCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGT  >  NC_000913/2841961‑2842103

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