BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,457	C→T	5.9%	G232S (GGC→AGC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,457	0	C	T	5.9%	1339.2 / 52.6	408	G232S (GGC→AGC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (158/226); new base T (6/18); total (164/244)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

Only 100 of 455 total aligned reads displayed.

TGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATA  >  CP009273/4106390‑4106524

tGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgcc                                                                     <  1:12956045/68‑1 (MQ=255)

  aCCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgc                                                                   <  1:29354601/68‑1 (MQ=255)

   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgctgcc                                                                  >  1:17307067/1‑68 (MQ=255)

   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgcc                                                                  >  1:14570217/1‑68 (MQ=255)

   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgcc                                                                  >  1:16460047/1‑68 (MQ=255)

   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgcc                                                                  >  1:19585530/1‑68 (MQ=255)

   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgcc                                                                  >  1:227392/1‑68 (MQ=255)

   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgcc                                                                  >  1:18161498/1‑68 (MQ=255)

   ccAAACAGCGCGGCCTGCTGCTCACCGGCGCTCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgcg                                                                  >  1:10785320/1‑67 (MQ=255)

    cAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCa                                                                 <  1:24579443/68‑1 (MQ=255)

     aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCaa                                                                >  1:6443309/1‑68 (MQ=255)

     aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCaa                                                                >  1:3198365/1‑68 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAAt                                                               <  1:10088069/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAAt                                                               <  1:6849521/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:4682577/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:25045002/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:20167090/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:8644377/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:21200437/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:29597562/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:10440040/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               >  1:29172013/1‑68 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:8750728/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               >  1:27886440/1‑68 (MQ=255)

       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATg                                                              >  1:9348801/1‑68 (MQ=255)

       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATg                                                              >  1:8362418/1‑68 (MQ=255)

       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATg                                                              >  1:11853314/1‑68 (MQ=255)

       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATg                                                              >  1:12526181/1‑68 (MQ=255)

       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATg                                                              >  1:24242876/1‑68 (MQ=255)

        cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGt                                                             >  1:5022196/1‑68 (MQ=255)

         aGCGCGGCCTGCTGGTCACCGGCGATCAGGGAGATTGGAATACGCGTGCCGCCTTTGACGCCAATGtt                                                            >  1:14463279/1‑68 (MQ=255)

            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTa                                                           <  1:13133667/66‑1 (MQ=255)

            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTa                                                           <  1:20008672/66‑1 (MQ=255)

            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAg                                                          >  1:75070/1‑67 (MQ=255)

            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGt                                                         <  1:14279697/68‑1 (MQ=255)

               gCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTg                                                      >  1:9038425/1‑68 (MQ=255)

                ccTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTg                                                      <  1:27791567/67‑1 (MQ=255)

                ccTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTg                                                      >  1:11578659/1‑67 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAc                                                    <  1:1505952/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGa                                                     >  1:23832940/1‑67 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:3689706/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:6111702/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:22278285/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:28846287/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:20890784/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:19222413/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:15349279/68‑1 (MQ=255)

                  tgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAcc                                                   <  1:23550300/68‑1 (MQ=255)

                   gctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCg                                                  <  1:6718628/68‑1 (MQ=255)

                     tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                                                 <  1:699649/67‑1 (MQ=255)

                     tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                                                 <  1:1901692/67‑1 (MQ=255)

                      ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                                               >  1:21837696/1‑68 (MQ=255)

                        tCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                                               >  1:12816308/1‑66 (MQ=255)

                        tCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                                               >  1:4111527/1‑66 (MQ=255)

                         cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                                            <  1:5647569/68‑1 (MQ=255)

                         cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                                            <  1:13495437/68‑1 (MQ=255)

                         cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                                            <  1:1786331/68‑1 (MQ=255)

                             ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                                        >  1:11264478/1‑68 (MQ=255)

                             ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                                        >  1:17132416/1‑68 (MQ=255)

                             ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                                        >  1:21691072/1‑68 (MQ=255)

                               cGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGga                                       <  1:17675187/67‑1 (MQ=255)

                                gATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaag                                     <  1:29811560/68‑1 (MQ=255)

                                     cGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacg                                 >  1:10022005/1‑67 (MQ=255)

                                      ggAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac                               <  1:7841985/68‑1 (MQ=255)

                                        agaTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgacg                              >  1:8224359/1‑67 (MQ=255)

                                        agaTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                             >  1:26309142/1‑68 (MQ=255)

                                        agaTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAACAACGACGc                             >  1:26892006/1‑68 (MQ=255)

                                          aTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                             >  1:22072358/1‑66 (MQ=255)

                                           ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                          <  1:24825371/68‑1 (MQ=255)

                                           ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                          <  1:18340845/68‑1 (MQ=255)

                                           ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCAGAAGAACGACGCAc                           <  1:9545685/67‑1 (MQ=255)

                                            tGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACtt                         <  1:11076311/68‑1 (MQ=255)

                                            tGGAATACGCGGGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACtt                         <  1:25286644/68‑1 (MQ=255)

                                             ggAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTc                        <  1:285783/68‑1 (MQ=255)

                                             ggAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTc                        <  1:16812669/68‑1 (MQ=255)

                                               aaTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGCCCGTATACTTCGGAAGAACGACGCACTTCTg                      >  1:27647654/1‑68 (MQ=255)

                                                aTACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                      <  1:5893820/67‑1 (MQ=255)

                                                aTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGCACGACGCACTTCtc                      >  1:5471161/1‑66 (MQ=255)

                                                 tACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgc                    <  1:3485189/68‑1 (MQ=255)

                                                 tACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgc                    >  1:22442951/1‑68 (MQ=255)

                                                  aCGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAATAACGACGCACTTCTGgca                   >  1:27456839/1‑68 (MQ=255)

                                                   cgcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGCACGACGCACGTCTGgcc                   >  1:26638844/1‑66 (MQ=255)

                                                    gcgTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc                 >  1:23317294/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                >  1:16152265/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                >  1:25723839/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                >  1:2306087/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                >  1:13843981/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACATCTGgcagca                >  1:10694441/1‑68 (MQ=255)

                                                           gccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgc          <  1:17488277/68‑1 (MQ=255)

                                                            ccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgc          >  1:28415882/1‑67 (MQ=255)

                                                               ttGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg      <  1:3387265/68‑1 (MQ=255)

                                                               ttGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg      <  1:16632006/68‑1 (MQ=255)

                                                               ttGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg      <  1:20640601/68‑1 (MQ=255)

                                                               tcgccgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTACTGGCAGCATCTCGCGCgg      <  1:20527984/66‑1 (MQ=255)

                                                                tGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa     <  1:19454108/68‑1 (MQ=255)

                                                                 gccgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGaa    <  1:23956799/68‑1 (MQ=255)

                                                                   cgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAata  >  1:15823677/1‑68 (MQ=255)

                                                                   cgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAat   <  1:7434619/67‑1 (MQ=255)

                                                                   cgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAat   <  1:19777584/67‑1 (MQ=255)

                                                                   cgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAat   <  1:10557107/67‑1 (MQ=255)

TGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATA  >  CP009273/4106390‑4106524

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