BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,815,253	G→T	S6I (AGC→ATC)	ygdH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,815,253	0	G	T	76.0%	691.2 / inf	350	S6I (AGC→ATC)	ygdH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (83/1); new base T (145/121); total (228/122)
Fisher's exact test for biased strand distribution p-value = 2.31e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 356 total aligned reads displayed.

CAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCT  >  REL606/2815179‑2815327

 aggcaggGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATcc                                                                          >  1:6905505/1‑76 (MQ=255)

    caggGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGc                                                                       >  1:400114/1‑76 (MQ=255)

      gggCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCtt                                                                     <  1:376641/76‑1 (MQ=255)

      gggCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCtt                                                                     <  2:7256985/76‑1 (MQ=255)

       ggCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTg                                                                    >  2:3905029/1‑76 (MQ=255)

        gCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTgg                                                                   <  2:4646718/76‑1 (MQ=255)

         cTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGc                                                                  <  2:627944/76‑1 (MQ=255)

         cTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGc                                                                  >  2:3227518/1‑76 (MQ=255)

          tATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCt                                                                 >  1:2711019/1‑76 (MQ=255)

           aTTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTc                                                                >  2:2876198/1‑76 (MQ=255)

           aTTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTc                                                                >  2:5945684/1‑76 (MQ=255)

             tGTAATCAAAGGGAATGATGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCa                                                              >  1:4181392/1‑76 (MQ=255)

             tGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCa                                                              <  1:4979451/76‑1 (MQ=255)

             tGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCa                                                              <  1:6497516/76‑1 (MQ=255)

               tAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATg                                                            >  2:4692201/1‑76 (MQ=255)

               tAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATg                                                            >  1:6703526/1‑76 (MQ=255)

                aaTCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATgg                                                           <  1:7265454/76‑1 (MQ=255)

                aaTCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATgg                                                           <  1:3252453/76‑1 (MQ=255)

                aaTCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATgg                                                           <  1:3097861/76‑1 (MQ=255)

                aaTCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATgg                                                           >  1:7221603/1‑76 (MQ=255)

                  tCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGat                                                         >  1:3583719/1‑76 (MQ=255)

                   cAAAGGGAATGTCGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGata                                                        >  1:6357112/1‑76 (MQ=255)

                   cAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGata                                                        >  1:8187973/1‑76 (MQ=255)

                   cAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGata                                                        >  2:7659702/1‑76 (MQ=255)

                   cAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGata                                                        >  1:3456764/1‑76 (MQ=255)

                   cAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGata                                                        >  1:5899096/1‑76 (MQ=255)

                    aaaGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGatat                                                       >  2:4740752/1‑76 (MQ=255)

                     aaGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATAtg                                                      >  1:1857279/1‑76 (MQ=255)

                     aaGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAtg                                                      >  1:395591/1‑76 (MQ=255)

                     aaGGCAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATAtg                                                      >  2:7519445/1‑76 (MQ=255)

                      aGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATAtgt                                                     <  1:2875859/76‑1 (MQ=255)

                            tGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGc                                               <  2:578333/76‑1 (MQ=255)

                            tGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGc                                               >  1:1423757/1‑76 (MQ=255)

                            tGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGc                                               <  2:41550/76‑1 (MQ=255)

                             gACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCa                                              >  1:5306870/1‑76 (MQ=255)

                             gACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCa                                              <  2:3787204/76‑1 (MQ=255)

                                gATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCt                                           <  2:1178975/76‑1 (MQ=255)

                                 atatTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTg                                          >  1:7115976/1‑76 (MQ=255)

                                 atatTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTg                                          >  2:2515881/1‑76 (MQ=255)

                                 atatTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTg                                          >  1:4475880/1‑76 (MQ=255)

                                 atatTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTg                                          >  1:3366969/1‑76 (MQ=255)

                                 atatTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTg                                          >  1:6744342/1‑76 (MQ=255)

                                   atTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGa                                        <  2:710082/76‑1 (MQ=255)

                                      cGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGt                                     <  2:2693375/76‑1 (MQ=255)

                                      cGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGt                                     >  1:1943264/1‑76 (MQ=255)

                                      cGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGt                                     >  2:8096612/1‑76 (MQ=255)

                                      cGTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGt                                     >  1:1032768/1‑76 (MQ=255)

                                      cGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGt                                     >  1:6082351/1‑76 (MQ=255)

                                      cGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGt                                     >  1:4678796/1‑76 (MQ=255)

                                      cGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGt                                     >  1:8077315/1‑76 (MQ=255)

                                       gTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTg                                    <  2:831491/76‑1 (MQ=255)

                                       gTCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTg                                    <  1:1577742/76‑1 (MQ=255)

                                        tCCCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTgg                                   <  1:5206675/76‑1 (MQ=255)

                                         cccATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGa                                  >  1:5999493/1‑76 (MQ=255)

                                          ccATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGat                                 <  1:4717693/76‑1 (MQ=255)

                                          ccATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGat                                 <  1:4587312/76‑1 (MQ=255)

                                            aTAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGatat                               <  1:7873632/76‑1 (MQ=255)

                                            aTAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGatat                               >  1:8290521/1‑76 (MQ=255)

                                             tAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATg                              >  1:7528082/1‑76 (MQ=255)

                                              aaGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGc                             >  1:5428475/1‑76 (MQ=255)

                                              aaGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGc                             >  2:5641706/1‑76 (MQ=255)

                                                ggAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCtt                           <  2:6223111/76‑1 (MQ=255)

                                                 gAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTa                          >  1:6392077/1‑76 (MQ=255)

                                                  aGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTaa                         <  2:2147165/76‑1 (MQ=255)

                                                  aGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTaa                         <  2:2954457/76‑1 (MQ=255)

                                                   gTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTaaa                        <  1:6212096/76‑1 (MQ=255)

                                                   gTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTaaa                        >  1:2396984/1‑76 (MQ=255)

                                                    tttttCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAAc                       <  2:549450/76‑1 (MQ=255)

                                                    tttttCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAAc                       <  1:5078203/76‑1 (MQ=255)

                                                      tttCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGc                     >  1:1981234/1‑76 (MQ=255)

                                                       ttCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCa                    >  2:3852551/1‑76 (MQ=255)

                                                       ttCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCa                    <  2:1678172/76‑1 (MQ=255)

                                                       ttCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCa                    >  1:1689285/1‑76 (MQ=255)

                                                        tCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCAc                   >  1:1140575/1‑76 (MQ=255)

                                                        tCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCAc                   >  1:6586612/1‑76 (MQ=255)

                                                         cTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCAcc                  >  1:2573780/1‑76 (MQ=255)

                                                         cTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCAcc                  >  2:1543283/1‑76 (MQ=255)

                                                         cTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCAcc                  >  1:5770013/1‑76 (MQ=255)

                                                          ttGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCg                 <  1:859532/76‑1 (MQ=255)

                                                           tGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGc                >  2:6403656/1‑76 (MQ=255)

                                                           tGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAaacgcaacgc                >  1:5640391/1‑76 (MQ=255)

                                                            gATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGcc               >  2:1087591/1‑76 (MQ=255)

                                                            gATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGcc               <  2:7949472/76‑1 (MQ=255)

                                                            gATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGcc               >  1:7661361/1‑76 (MQ=255)

                                                             aTTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCca              >  1:4296556/1‑76 (MQ=255)

                                                             aTTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCca              <  2:2042184/76‑1 (MQ=255)

                                                             aTTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCca              <  2:341913/76‑1 (MQ=255)

                                                             aTTACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCAACGCca              >  2:634605/1‑76 (MQ=255)

                                                               tACACATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCcagc            <  2:425444/76‑1 (MQ=255)

                                                               tACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCcagc            >  1:6274611/1‑76 (MQ=255)

                                                               tACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAAAGCACCGCcagc            >  1:5573712/1‑76 (MQ=255)

                                                                acacATATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCcagca           <  2:6719421/76‑1 (MQ=255)

                                                                 cacaTATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCcagcag          <  1:3671404/76‑1 (MQ=255)

                                                                  acaTATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCagcagc         >  2:4921557/1‑76 (MQ=255)

                                                                   caTATTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCg        <  1:738042/76‑1 (MQ=255)

                                                                   caTATTATCCCGCTCGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCg        <  2:5228211/76‑1 (MQ=255)

                                                                     tattatCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGAc      >  2:1403393/1‑76 (MQ=255)

                                                                      atTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGAcc     <  1:4853956/76‑1 (MQ=255)

                                                                      atTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGAcc     >  1:2220528/1‑76 (MQ=255)

                                                                      atTATCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGAcc     >  1:595870/1‑76 (MQ=255)

CAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCT  >  REL606/2815179‑2815327

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`