BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,456	C→T	19.8%	G232D (GGC→GAC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,456	0	C	T	19.8%	937.4 / 229.3	418	G232D (GGC→GAC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (142/193); new base T (29/54); total (171/247)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 446 total aligned reads displayed.

ACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAA  >  CP009273/4106392‑4106522

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

 ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgtc                                                                >  1:8378811/1‑66 (MQ=255)

 ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgc                                                                 <  1:22443271/67‑1 (MQ=255)

 ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgc                                                                 <  1:3342784/67‑1 (MQ=255)

 ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgc                                                                 <  1:8050145/67‑1 (MQ=255)

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

    aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAAt                                                             <  1:26514742/68‑1 (MQ=255)

    aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAAt                                                             <  1:28714665/68‑1 (MQ=255)

    aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAAt                                                             <  1:7570026/68‑1 (MQ=255)

    aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAAt                                                             <  1:10742323/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

      cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGt                                                           >  1:16876578/1‑68 (MQ=255)

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

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

      cAGCGCGGCCTGCTGGTCACCGGAGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGt                                                           >  1:17105070/1‑68 (MQ=255)

       aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGtt                                                          >  1:2125671/1‑68 (MQ=255)

          gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAg                                                        <  1:14907216/67‑1 (MQ=255)

          gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGt                                                       <  1:26859164/68‑1 (MQ=255)

          gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTa                                                         <  1:21003394/66‑1 (MQ=255)

          gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGt                                                       >  1:16449096/1‑68 (MQ=255)

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

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

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

               ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGa                                                   <  1:7166791/67‑1 (MQ=255)

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

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

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

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

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

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

                   tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGt                                               <  1:12036456/67‑1 (MQ=255)

                   tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGt                                               <  1:24829939/67‑1 (MQ=255)

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

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

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

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

                      tCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtata                                            <  1:1292094/67‑1 (MQ=255)

                       cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACt                                          <  1:15085559/68‑1 (MQ=255)

                       cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACt                                          <  1:27516950/68‑1 (MQ=255)

                       cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACt                                          <  1:4353772/68‑1 (MQ=255)

                       cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACt                                          <  1:24423274/68‑1 (MQ=255)

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

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

                        aCCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                                          >  1:23125096/1‑67 (MQ=255)

                           ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCgg                                      >  1:6782407/1‑68 (MQ=255)

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

                                tCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCGAATGTTAGTCTGACCGTATACTTCGgaaga                                  >  1:891637/1‑67 (MQ=255)

                                 cccGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAac                                <  1:27687493/68‑1 (MQ=255)

                                 cccGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAac                                <  1:14499111/68‑1 (MQ=255)

                                  ccGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacg                               <  1:18774648/68‑1 (MQ=255)

                                    ggAGATTGGAATACGCGTGCCTCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac                             <  1:25200064/68‑1 (MQ=255)

                                      agaTTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGtata                                            >  1:6606892/1‑51 (MQ=255)

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

                                         ttGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                        <  1:14590095/68‑1 (MQ=255)

                                         ttGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAc                         <  1:17874662/67‑1 (MQ=255)

                                         ttGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                        <  1:15621295/68‑1 (MQ=255)

                                         ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAc                         <  1:22957665/67‑1 (MQ=255)

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

                                             aaTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                    >  1:15831747/1‑68 (MQ=255)

                                             aaTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                    >  1:29019950/1‑68 (MQ=255)

                                              aTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                    >  1:10311756/1‑67 (MQ=255)

                                              aTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                    >  1:27140573/1‑67 (MQ=255)

                                              aTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                    >  1:22855291/1‑67 (MQ=255)

                                               tACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgc                  <  1:13916574/68‑1 (MQ=255)

                                                aCGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAATTCTGgc                  >  1:10440952/1‑67 (MQ=255)

                                                aCGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                 <  1:21530220/68‑1 (MQ=255)

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

                                                  gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc               >  1:25478633/1‑68 (MQ=255)

                                                  gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc               >  1:5695247/1‑68 (MQ=255)

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

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

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

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

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

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

                                                   cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTAAACTTCGGAAGAACGACGCACTACTGgcagca              >  1:7388404/1‑68 (MQ=255)

                                                   cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGGATACTTCGGAAGAACGCCGCCCTGCTGgcagca              >  1:11019900/1‑68 (MQ=255)

                                                    gTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca              <  1:17961390/67‑1 (MQ=255)

                                                    gTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAt             <  1:27936718/68‑1 (MQ=255)

                                                      gccgccTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAtct           <  1:22158272/68‑1 (MQ=255)

                                                      gccgccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAtc            <  1:19087243/67‑1 (MQ=255)

                                                        cgccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAtctc          >  1:24163917/1‑67 (MQ=255)

                                                         gccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgc        >  1:12254649/1‑68 (MQ=255)

                                                          ccTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgc        >  1:18369934/1‑67 (MQ=255)

                                                             ttGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg    <  1:20471003/68‑1 (MQ=255)

                                                             ttGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg    <  1:292322/68‑1 (MQ=255)

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

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

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

                                                              tGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa   <  1:28161989/68‑1 (MQ=255)

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

ACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAA  >  CP009273/4106392‑4106522

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: `‑`