BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,459	C→T	14.3%	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	14.3%	1117.1 / 148.3	418	G231D (GGC→GAC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (153/205); new base T (18/42); total (171/247)
Fisher's exact test for biased strand distribution p-value = 6.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 458 total aligned reads displayed.

ACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATAT  >  CP009273/4106392‑4106525

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

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

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

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

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

   aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCaa                                                                 >  1:17957593/1‑68 (MQ=255)

   aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCaa                                                                 >  1:14467666/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

     aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATg                                                               >  1:1369870/1‑68 (MQ=255)

     aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATg                                                               <  1:15752006/68‑1 (MQ=255)

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

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

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

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

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

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

          gcgGCCTGCTGGTCACCAGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTa                                                            <  1:26671299/66‑1 (MQ=255)

             gCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTg                                                       >  1:10480425/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

               ctcctgGTCACCGGCGATCCCGGAGATTGTACTACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGAc                                                     <  1:5436542/65‑1 (MQ=255)

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

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

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

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

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

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

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

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

                      tCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGtat                                                >  1:2773443/1‑66 (MQ=255)

                       caacGGCGATCCCTGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                                             <  1:8767147/65‑1 (MQ=255)

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

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

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

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

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

                        aCCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACt                                             <  1:24351528/67‑1 (MQ=255)

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

                        aCCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCGGACCGTATACt                                             >  1:13219175/1‑67 (MQ=255)

                           ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCgg                                         >  1:16396970/1‑68 (MQ=255)

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

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

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

                                      agaTTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                              >  1:10968122/1‑68 (MQ=255)

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

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

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

                                        aTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                     <  1:19184510/43‑1 (MQ=255)

                                         ttGTAATACGCGTGCCGCATTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                           <  1:4856929/68‑1 (MQ=255)

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

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

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

                                            gAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                           >  1:26209160/1‑65 (MQ=255)

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

                                              aTACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                       >  1:27323402/1‑67 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                 tgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAata   <  1:12519660/67‑1 (MQ=255)

                                                                 cgccACTGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAat    <  1:28467053/67‑1 (MQ=255)

                                                                  gccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAatat  <  1:16564826/68‑1 (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: `‑`