BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,634	T→C	5.5%	K173E (AAA→GAA)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,634	0	T	C	5.5%	1948.1 / 57.6	544	K173E (AAA→GAA)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base T (340/173); new base C (19/11); total (359/185)
Fisher's exact test for biased strand distribution p-value = 8.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 628 total aligned reads displayed.

TGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCT  >  CP009273/4106568‑4106701

tGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                                                                     >  1:27731236/1‑67 (MQ=255)

tGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                                                                     >  1:28639818/1‑67 (MQ=255)

tGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                                                                     >  1:1383378/1‑67 (MQ=255)

tGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                                                                     >  1:17462844/1‑67 (MQ=255)

tGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                                                                     >  1:2746476/1‑67 (MQ=255)

tGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                                                                     >  1:13634535/1‑67 (MQ=255)

    gAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGa                                                                >  1:24768892/1‑68 (MQ=255)

      acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATa                                                              >  1:14285832/1‑68 (MQ=255)

      acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATa                                                              >  1:25816855/1‑68 (MQ=255)

      acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATa                                                              >  1:2437402/1‑68 (MQ=255)

      acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATa                                                              >  1:26729502/1‑68 (MQ=255)

      acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACAAGGACACGGCCCTGAGTCATTTTCCAGATa                                                              >  1:3642222/1‑68 (MQ=255)

        aacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATaa                                                             >  1:2609719/1‑67 (MQ=255)

         acaTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAg                                                            >  1:21123022/1‑67 (MQ=255)

         acaTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGc                                                           <  1:25509694/68‑1 (MQ=255)

         acaTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGc                                                           <  1:7984113/68‑1 (MQ=255)

           aTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGcc                                                          >  1:23813758/1‑67 (MQ=255)

           aTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGcc                                                          >  1:15592049/1‑67 (MQ=255)

             ggTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg                                                       <  1:19710406/68‑1 (MQ=255)

             ggTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg                                                       <  1:15320076/68‑1 (MQ=255)

             ggTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg                                                       <  1:14812257/68‑1 (MQ=255)

             ggTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg                                                       <  1:4316373/68‑1 (MQ=255)

             ggTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg                                                       >  1:8435666/1‑68 (MQ=255)

             ggTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg                                                       <  1:1215704/68‑1 (MQ=255)

             ggTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg                                                       <  1:11411060/68‑1 (MQ=255)

                actagagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:7242045/4‑67 (MQ=255)

                actagagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     <  1:2638211/64‑1 (MQ=255)

                actagagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:9060026/4‑67 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:26947764/1‑67 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:23588041/1‑67 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:22322602/1‑67 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:21423955/1‑67 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:28094272/1‑67 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:3117716/1‑67 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     <  1:8694810/67‑1 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:13313820/1‑67 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:14055760/1‑67 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:10474332/1‑67 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:16344916/1‑67 (MQ=255)

                aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa                                                     >  1:6708105/1‑67 (MQ=255)

                   agagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCaa                                                 <  1:5394129/68‑1 (MQ=255)

                   agagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCa                                                  >  1:20184673/1‑67 (MQ=255)

                   agagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCa                                                  >  1:7177831/1‑67 (MQ=255)

                   agagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCa                                                  >  1:16156028/1‑67 (MQ=255)

                   agagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCa                                                  >  1:22585829/1‑67 (MQ=255)

                    gagGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCaa                                                 >  1:16426585/1‑67 (MQ=255)

                     agGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCGACGTATCAAcc                                               >  1:18388492/1‑68 (MQ=255)

                     agGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAcc                                               >  1:9679835/1‑68 (MQ=255)

                     agGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAcc                                               >  1:24954459/1‑68 (MQ=255)

                     agGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAcc                                               >  1:27156561/1‑68 (MQ=255)

                     agGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAcc                                               >  1:17824151/1‑68 (MQ=255)

                     agGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAcc                                               >  1:18254678/1‑68 (MQ=255)

                     agGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAcc                                               >  1:20811464/1‑68 (MQ=255)

                     agGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAcc                                               >  1:3208035/1‑68 (MQ=255)

                      ggCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAcc                                               <  1:4627791/67‑1 (MQ=255)

                      ggCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCg                                              >  1:4235884/1‑68 (MQ=255)

                        cGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGt                                             <  1:10053167/67‑1 (MQ=255)

                          ttGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTAcc                                          >  1:11869735/1‑68 (MQ=255)

                              tgtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAAc                                      >  1:22003348/1‑68 (MQ=255)

                              tgtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAAc                                      >  1:19442919/1‑68 (MQ=255)

                               gtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACa                                     <  1:10882785/68‑1 (MQ=255)

                                   tCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCaa                                 >  1:1911570/1‑68 (MQ=255)

                                   tCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCaa                                 >  1:17099772/1‑68 (MQ=255)

                                   tCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCaa                                 >  1:5189937/1‑68 (MQ=255)

                                   tCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCaa                                 >  1:18108219/1‑68 (MQ=255)

                                   tCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCaa                                 >  1:23087644/1‑68 (MQ=255)

                                   tCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCaa                                 >  1:21504609/1‑68 (MQ=255)

                                    cGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAAt                                >  1:5904232/1‑68 (MQ=255)

                                    cGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAAt                                >  1:20416820/1‑68 (MQ=255)

                                        cacaTGTACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAAGCGTACCAAACAGCAATTcac                            <  1:23385840/68‑1 (MQ=255)

                                        cacaTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcac                            <  1:10303238/68‑1 (MQ=255)

                                        cacaTCTACAAGGCCCTGAGTAATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcac                            <  1:341264/68‑1 (MQ=255)

                                           aTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTcaccac                         >  1:6500960/1‑68 (MQ=255)

                                             ggACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCacga                       <  1:1167993/68‑1 (MQ=255)

                                             ggACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCacga                       <  1:7807463/68‑1 (MQ=255)

                                             ggACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCacg                        >  1:29013078/1‑67 (MQ=255)

                                               acacGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCacgacg                     >  1:112127/1‑68 (MQ=255)

                                               acacGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCacgacg                     >  1:6107422/1‑68 (MQ=255)

                                               acacGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTAACCacgacg                     >  1:24096520/1‑68 (MQ=255)

                                                cacGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGt                    <  1:12848076/68‑1 (MQ=255)

                                                cacGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGt                    <  1:4988546/68‑1 (MQ=255)

                                                cacGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGt                    <  1:1592594/68‑1 (MQ=255)

                                                 acGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGt                    >  1:8315726/1‑67 (MQ=255)

                                                 acGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGt                    >  1:7571053/1‑67 (MQ=255)

                                                 acGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATGCACCACGACGt                    >  1:6376300/1‑67 (MQ=255)

                                                     cccTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCAc                <  1:16603773/67‑1 (MQ=255)

                                                      ccTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCAc                >  1:3768124/1‑66 (MQ=255)

                                                       cTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCt             <  1:9493291/68‑1 (MQ=255)

                                                       cTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCt             <  1:470363/68‑1 (MQ=255)

                                                       cTGAGTCATTTTCCAGATAAGCCACGTACCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCt             <  1:25417156/68‑1 (MQ=255)

                                                           gTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcgcg         >  1:28435300/1‑68 (MQ=255)

                                                           gTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcgcg         >  1:26549331/1‑68 (MQ=255)

                                                           gTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcgcg         >  1:18704348/1‑68 (MQ=255)

                                                           gTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcgcg         >  1:3926640/1‑68 (MQ=255)

                                                           gTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcgcg         >  1:14629538/1‑68 (MQ=255)

                                                                tttCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGagag     >  1:12659951/1‑67 (MQ=255)

                                                                tttCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGc    <  1:29245622/68‑1 (MQ=255)

                                                                 ttCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGcc   <  1:21708030/68‑1 (MQ=255)

                                                                  tCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCt  >  1:9328574/1‑68 (MQ=255)

                                                                  tCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCt  >  1:15188644/1‑68 (MQ=255)

TGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCT  >  CP009273/4106568‑4106701

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