BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,459	0	C	T	5.1%	1627.0 / 46.4	484	G231D (GGC→GAC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (186/272); new base T (12/13); total (198/286)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 536 total aligned reads displayed.

CCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATAT  >  CP009273/4106393‑4106525

ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgtcgcc                                                                   >  1:15210145/1‑68 (MQ=255)

ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgtcgcc                                                                   <  1:27576172/68‑1 (MQ=255)

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

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

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

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

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

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

   cccAGCGCGGCCTGCTTGTCACCGGCGATCCCGGAGATTGGAATAAGCGTGCCGCCTTTGCCGCCAAt                                                                <  1:22976485/66‑1 (MQ=255)

   aaCAGCGCGGGGTGCTGGTCAGCGGCGATCCGGGAGATTGGAATACGCGTGCCGCCTTTgccgccttt                                                                <  1:31558698/68‑4 (MQ=255)

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

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

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

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

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

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

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

    aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATg                                                               >  1:9351385/1‑68 (MQ=255)

    aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATg                                                               >  1:25963914/1‑68 (MQ=255)

     cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGt                                                              >  1:8540968/1‑68 (MQ=255)

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

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

      aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGtt                                                             >  1:19310026/1‑68 (MQ=255)

      aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCGTTGCCGTCAATGtt                                                             >  1:2179536/1‑68 (MQ=255)

       gcgcGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTa                                                            >  1:24800057/1‑68 (MQ=255)

         gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTa                                                            <  1:17508597/66‑1 (MQ=255)

         gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTa                                                            <  1:7300352/66‑1 (MQ=255)

         gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAg                                                           >  1:30067669/1‑67 (MQ=255)

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

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

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

            gCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTg                                                       >  1:17079085/1‑68 (MQ=255)

              ctgttgGTCACGGGCGATCCCGGAGATTGGAATACGCGTGCCTCCTTTGCCGCGAATGTTAGTCTGAc                                                     <  1:4025961/68‑1 (MQ=255)

              ctgcttGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                     <  1:11323503/68‑1 (MQ=255)

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

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

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

              ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGAc                                                     <  1:21893533/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

               tgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAcc                                                    <  1:6817915/68‑1 (MQ=255)

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

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

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

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

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

                     tCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGtat                                                >  1:25801688/1‑66 (MQ=255)

                     tCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGtat                                                >  1:27212040/1‑66 (MQ=255)

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

                      cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATAc                                              <  1:28791634/67‑1 (MQ=255)

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

                        ccGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTc                                           <  1:13414161/68‑1 (MQ=255)

                          ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACAGTATACTTCgg                                         >  1:26474708/1‑68 (MQ=255)

                            cGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGga                                        <  1:21113533/67‑1 (MQ=255)

                                     agaTTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                              >  1:28423695/1‑68 (MQ=255)

                                     agaTTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                              >  1:6247010/1‑68 (MQ=255)

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

                                        ttGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAc                            <  1:30717881/67‑1 (MQ=255)

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

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

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

                                         tGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACtt                          <  1:6554262/68‑1 (MQ=255)

                                         tGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACtt                          <  1:4875505/68‑1 (MQ=255)

                                          ggAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTc                         <  1:14717719/68‑1 (MQ=255)

                                          ggAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTc                         <  1:5482934/68‑1 (MQ=255)

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

                                           gAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACAACGCACt                           >  1:25217042/1‑65 (MQ=255)

                                             aTACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                       <  1:14346941/67‑1 (MQ=255)

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

                                                 gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGCCCGTATACTTCGGAAGAACGACGCACTTCTGgcag                   >  1:11553164/1‑67 (MQ=255)

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

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

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

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

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

                                                  cgTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                 >  1:26335466/1‑68 (MQ=255)

                                                       cgccTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcg            <  1:13854113/68‑1 (MQ=255)

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

                                                            ttGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg       <  1:11932071/68‑1 (MQ=255)

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

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

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

                                                              gtcgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGaa     <  1:26764639/66‑1 (MQ=255)

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

                                                               tcgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAat    >  1:30359021/2‑68 (MQ=255)

                                                                ctccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAat    <  1:16680776/65‑1 (MQ=255)

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

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

                                                                 gccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAatat  <  1:20238655/68‑1 (MQ=255)

                                                                 gccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAatat  <  1:19444564/68‑1 (MQ=255)

                                                                 gccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAatat  <  1:17953055/68‑1 (MQ=255)

CCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATAT  >  CP009273/4106393‑4106525

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