BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,944	A→C	6.7%	D69E (GAT→GAG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,944	0	A	C	6.7%	1716.6 / 82.1	478	D69E (GAT→GAG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base A (191/255); new base C (14/18); total (205/273)
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.64e-01

Only 100 of 559 total aligned reads displayed.

CCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGT  >  CP009273/4106881‑4107010

ccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCg                                                                 >  1:25648666/1‑67 (MQ=255)

    aCAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCtt                                                             >  1:2540415/1‑67 (MQ=255)

    aCAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCtt                                                             >  1:6370004/1‑67 (MQ=255)

     cAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTc                                                           >  1:12789937/1‑68 (MQ=255)

      aaTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTc                                                           >  1:23241480/1‑67 (MQ=255)

      aaTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCg                                                          <  1:22664575/68‑1 (MQ=255)

      aaTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCg                                                          <  1:28434534/68‑1 (MQ=255)

         aGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCa                                                       >  1:25884164/1‑68 (MQ=255)

          gTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCa                                                       >  1:10267132/1‑67 (MQ=255)

           tGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGc                                                     >  1:24025142/1‑68 (MQ=255)

            ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGc                                                     >  1:29007074/1‑67 (MQ=255)

            ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGc                                                     >  1:19609251/1‑67 (MQ=255)

            ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGc                                                     >  1:26649180/1‑67 (MQ=255)

            ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGc                                                     >  1:15708846/1‑67 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATCTCGGCTTTCGCCAGCAc                                                   <  1:4466233/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCAc                                                   <  1:18086389/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                   <  1:22894866/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                   <  1:19137101/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                   <  1:8585847/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                   <  1:1988070/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                   <  1:16531179/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                   <  1:22340718/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                   <  1:8741663/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                   <  1:26443871/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                   <  1:24277959/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                   <  1:24682860/67‑1 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACt                                                  >  1:10994923/1‑68 (MQ=255)

              tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACt                                                  >  1:9312566/1‑68 (MQ=255)

               ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACtt                                                 <  1:13492186/68‑1 (MQ=255)

               ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACtt                                                 <  1:10025968/68‑1 (MQ=255)

               ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACtt                                                 <  1:19457183/68‑1 (MQ=255)

               ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACt                                                  >  1:9627442/1‑67 (MQ=255)

                  aCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTa                                              <  1:17216428/68‑1 (MQ=255)

                    gCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATCTCGGCTTTCGCCAGCACTTCTAcc                                            <  1:14615765/68‑1 (MQ=255)

                       ggTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGc                                         <  1:12668127/68‑1 (MQ=255)

                       ggTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGc                                         <  1:28846873/68‑1 (MQ=255)

                         ttCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGt                                       >  1:16773449/1‑68 (MQ=255)

                          tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTg                                      >  1:27388495/1‑68 (MQ=255)

                          tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTg                                      <  1:11578355/68‑1 (MQ=255)

                              aaTACCGATTGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                  <  1:20736177/68‑1 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                    >  1:15913029/1‑66 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                    >  1:3609642/1‑66 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                    >  1:14341611/1‑66 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                    >  1:6593693/1‑66 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                    >  1:11765472/1‑66 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAg                                   <  1:7276379/67‑1 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                  <  1:1130942/68‑1 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                  <  1:8408329/68‑1 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                  <  1:12160333/68‑1 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                  <  1:10648709/68‑1 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                  <  1:15485665/68‑1 (MQ=255)

                              aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                  <  1:17852890/68‑1 (MQ=255)

                                tACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCtt                                >  1:8027913/1‑68 (MQ=255)

                                 aCCGATAGCTGCAATTTGATCGGAACTGATCTCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                               >  1:3172918/1‑68 (MQ=255)

                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGtttt                               >  1:23696288/1‑68 (MQ=255)

                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                               >  1:5677360/1‑68 (MQ=255)

                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                               >  1:23857767/1‑68 (MQ=255)

                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                               >  1:15048626/1‑68 (MQ=255)

                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                               >  1:21062756/1‑68 (MQ=255)

                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                               >  1:13015089/1‑68 (MQ=255)

                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                               >  1:18857508/1‑68 (MQ=255)

                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                               >  1:25087937/1‑68 (MQ=255)

                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                               >  1:1852989/1‑68 (MQ=255)

                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                               >  1:18271906/1‑68 (MQ=255)

                                  ccGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTg                              >  1:14554353/1‑68 (MQ=255)

                                  ccGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTg                              >  1:3776209/1‑68 (MQ=255)

                                   cGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                             <  1:13426483/68‑1 (MQ=255)

                                   cGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                             >  1:7826577/1‑68 (MQ=255)

                                          tGCAATTTGATCGGAACTGATCTCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGccc                      <  1:24855608/68‑1 (MQ=255)

                                          tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTggt                        <  1:4215397/66‑2 (MQ=255)

                                          tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGAcc                       <  1:7145516/67‑1 (MQ=255)

                                             aaTTTGATCGGAACTGATCTCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                   <  1:26844569/68‑1 (MQ=255)

                                             aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                   <  1:6697618/68‑1 (MQ=255)

                                               tttGATCGCAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAtt                 <  1:1399806/68‑1 (MQ=255)

                                                  gATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGATCTTTGGGTGGCCCAGATTTcc              <  1:16351027/68‑1 (MQ=255)

                                                  gATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc              <  1:1740737/68‑1 (MQ=255)

                                                  gATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc              <  1:22133381/68‑1 (MQ=255)

                                                    tCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAt            >  1:4962963/1‑68 (MQ=255)

                                                        aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg         <  1:21670512/67‑1 (MQ=255)

                                                         aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTggg        >  1:19693819/1‑67 (MQ=255)

                                                         aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt       <  1:17679468/68‑1 (MQ=255)

                                                         aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt       <  1:14853248/68‑1 (MQ=255)

                                                         aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGt       >  1:20128318/1‑68 (MQ=255)

                                                            gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCgt    <  1:16269828/68‑1 (MQ=255)

                                                            gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt    <  1:27651379/68‑1 (MQ=255)

                                                              tatCGGCTTTCGCCAGCCCTTCTACCAGCGTGGAGC‑TTGGGTGACCCAGATTTCCATTGGGTCgtgt  >  1:21923469/1‑67 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTTGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtg   <  1:20018986/67‑1 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGATCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt  >  1:5940847/1‑68 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCgtgt  >  1:21320040/1‑68 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt  <  1:17025244/68‑1 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt  >  1:12332694/1‑68 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt  >  1:8965900/1‑68 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt  >  1:2050352/1‑68 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt  <  1:23086101/68‑1 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt  <  1:27067483/68‑1 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt  <  1:28222637/68‑1 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt  <  1:3031423/68‑1 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt  <  1:4882685/68‑1 (MQ=255)

                                                              tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGCTCgtgt  >  1:13800947/1‑68 (MQ=255)

                                                              tatCGGCGTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGAGACgtgt  >  1:5349945/1‑68 (MQ=255)

CCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGT  >  CP009273/4106881‑4107010

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