BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,459	C→T	6.9%	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	6.9%	2447.6 / 118.3	753	G231D (GGC→GAC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (301/397); new base T (22/30); total (325/428)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 835 total aligned reads displayed.

ACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATAT  >  CP009273/4106392‑4106525

caacatccataccctggactgggacgacaaaatgctggaagtgctggatattccgcgcgagatgctgc                                                                    >  1:9018572‑M1/67‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

  acatccataccctggactgggacgacaaaatgctggaagtgctggatattccgcgcgagatgctGCCa                                                                  <  1:12238872‑M1/4‑1 (MQ=255)

    aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAAt                                                                <  1:36630895/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

     aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTgccgcc                                                                   >  1:31163648/1‑64 (MQ=255)

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

      cAGCGCGGCCTGCTGGTCACCGTCGATCCCGGAGATTGGAATACGCGTGCCGCCATTGCCGACaaatt                                                              >  1:24719233/1‑65 (MQ=255)

      cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGt                                                              >  1:19949326/1‑68 (MQ=255)

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

          gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTa                                                            <  1:13198927/66‑1 (MQ=255)

          gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAg                                                           <  1:46684455/67‑1 (MQ=255)

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

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

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

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

               ctgctgGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                     <  1:15203837/68‑1 (MQ=255)

               ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGAc                                                     <  1:30223347/68‑1 (MQ=255)

               ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGAc                                                     <  1:32315609/68‑1 (MQ=255)

               ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGAc                                                     <  1:40303043/68‑1 (MQ=255)

               ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGAc                                                     <  1:4088479/68‑1 (MQ=255)

               ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGAc                                                     <  1:50721227/68‑1 (MQ=255)

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

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

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

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

               ctgctgGTCACCGGCGATCCCGGAGATTGGAATAAGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                     <  1:26820299/68‑1 (MQ=255)

                 gctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCg                                                   >  1:1714113/1‑68 (MQ=255)

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

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

                    ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGtat                                                <  1:5638479/68‑1 (MQ=255)

                    ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGtat                                                >  1:46159594/1‑68 (MQ=255)

                    ggTCACCGGCGATCCCGGAGATTGGAATACGCGTCCCGCCTTTGCCGCCAATGTTAGGCTGACCGtat                                                >  1:9433185/1‑68 (MQ=255)

                      tCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGtat                                                >  1:42719298/1‑66 (MQ=255)

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

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

                       cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATAc                                              <  1:16950770/67‑1 (MQ=255)

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

                       cACCGGAGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                                             <  1:48835675/68‑1 (MQ=255)

                        aCCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACt                                             >  1:21878444/1‑67 (MQ=255)

                                tCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaaga                                     <  1:37578582/67‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                          tGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACtt                          <  1:13910659/68‑1 (MQ=255)

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

                                            gAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                           >  1:6897679/1‑65 (MQ=255)

                                             aaTACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                       >  1:38016768/1‑68 (MQ=255)

                                              aTACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                       >  1:2751445/1‑67 (MQ=255)

                                              aTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                       <  1:41249672/67‑1 (MQ=255)

                                                aCGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                    >  1:20465109/1‑68 (MQ=255)

                                                  gcgTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc                  >  1:42220933/1‑68 (MQ=255)

                                                   cgTGCCGCCTTTGCTGCCATTGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                 >  1:47038445/1‑68 (MQ=255)

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

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

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

                                                   cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGATGACCGCCGCACGTCTGgcagca                 >  1:13504754/1‑68 (MQ=255)

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

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

                                                         gccTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgc           >  1:19393655/1‑68 (MQ=255)

                                                            tttgctGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgcgcg        >  1:18841081‑M1/7‑68 (MQ=255)

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

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

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

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

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

                                                              tGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa      >  1:48208766/1‑68 (MQ=255)

                                                              tGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa      >  1:18332262/1‑68 (MQ=255)

                                                                ctGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGTCAGCATCTCGCGCGGAat    <  1:30612754‑M1/66‑1 (MQ=255)

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

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

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

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

                                                                  GCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAatat  >  1:5801342‑M1/1‑68 (MQ=255)

                                                                  GCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAatat  >  1:31458859‑M1/1‑68 (MQ=255)

                                                                  GCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAatat  >  1:22552712‑M1/1‑68 (MQ=255)

                                                                  GCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAatat  >  1:38742858‑M1/1‑68 (MQ=255)

ACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATAT  >  CP009273/4106392‑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: `‑`