BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,438	A→G	87.2%	I238T (ATT→ACT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,438	0	A	G	87.2%	1782.4 / 214.8	636	I238T (ATT→ACT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base A (39/42); new base G (241/312); total (281/355)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 638 total aligned reads displayed.

CTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTG  >  CP009273/4106373‑4106504

  tcttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTAcg                                                                <  1:24553703/68‑1 (MQ=255)

   cttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTAcgc                                                               >  1:27005335/1‑68 (MQ=255)

    tttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTAcgcg                                                              >  1:14064849/1‑68 (MQ=255)

     ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGt                                                             >  1:40461390/1‑68 (MQ=255)

     ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGt                                                             >  1:15506598/1‑68 (MQ=255)

     ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGt                                                             >  1:27977983/1‑68 (MQ=255)

     ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGt                                                             >  1:11050262/1‑68 (MQ=255)

          gCACAACTGACCAAACCGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTgccg                                                         >  1:29143112/1‑67 (MQ=255)

          gCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTgccg                                                         >  1:37589150/1‑67 (MQ=255)

          gCACAACTGACCAAAAAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTgccgc                                                        <  1:687190/68‑1 (MQ=255)

               aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCttt                                                    >  1:12134907/1‑67 (MQ=255)

               aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCttt                                                    >  1:34621941/1‑67 (MQ=255)

               aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCttt                                                    >  1:6391473/1‑67 (MQ=255)

               aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCttt                                                    >  1:34709011/1‑67 (MQ=255)

                 tGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTgcc                                                 <  1:29559809/68‑1 (MQ=255)

                    ccAACCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTgccgc                                               <  1:7707186/67‑1 (MQ=255)

                    ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTgccgcc                                              >  1:23854749/1‑68 (MQ=255)

                    ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTgccgcc                                              >  1:24692633/1‑68 (MQ=255)

                      aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCaa                                            >  1:25196766/1‑68 (MQ=255)

                       aaCAGGGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAAt                                           <  1:33660211/68‑1 (MQ=255)

                       aaCAGCGCGGCCTGCTTGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAAt                                           <  1:32025604/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                        aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATg                                          <  1:14371752/68‑1 (MQ=255)

                        aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATg                                          >  1:41833043/1‑68 (MQ=255)

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

                         cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGt                                         >  1:15794209/1‑68 (MQ=255)

                         cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGt                                         >  1:16250342/1‑68 (MQ=255)

                         cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGt                                         >  1:29962161/1‑68 (MQ=255)

                         cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGt                                         >  1:35410071/1‑68 (MQ=255)

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

                         cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGAATGGAGTACGCGTGCCGCCTTTGCCCCCAATGt                                         >  1:36296103/1‑68 (MQ=255)

                          aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGtt                                        >  1:39248529/1‑68 (MQ=255)

                          aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGtt                                        >  1:16816380/1‑68 (MQ=255)

                             gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTa                                       <  1:8658005/66‑1 (MQ=255)

                             gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAg                                      >  1:8384777/1‑67 (MQ=255)

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

                              cgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAg                                      >  1:3858718/1‑66 (MQ=255)

                                gCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTg                                  >  1:26082586/1‑68 (MQ=255)

                                gCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTg                                  >  1:7396369/1‑68 (MQ=255)

                                  ctgcttGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                <  1:19364855/68‑1 (MQ=255)

                                  ctgctgGTCTCCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                <  1:28399962/68‑1 (MQ=255)

                                  ctgctgGTCACGGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                <  1:10205667/68‑1 (MQ=255)

                                  ctgctgGTCACCGGCGATCCCGGAGTTTGGAGTACGGGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                <  1:17853979/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGGGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                <  1:24236518/68‑1 (MQ=255)

                                      tgGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                             <  1:18813209/67‑1 (MQ=255)

                                      tgGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                             <  1:8975505/67‑1 (MQ=255)

                                       ggTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                           >  1:25541631/1‑68 (MQ=255)

                                          cACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                        <  1:26693683/68‑1 (MQ=255)

                                           aCCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                        >  1:21277539/1‑67 (MQ=255)

                                             cGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCg                     <  1:37721256/68‑1 (MQ=255)

                                             cGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCg                     <  1:35023755/68‑1 (MQ=255)

                                              ggCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                    >  1:36447040/1‑68 (MQ=255)

                                              ggCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                    >  1:43633556/1‑68 (MQ=255)

                                              ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                    >  1:23344512/1‑68 (MQ=255)

                                                 gATCCCGGAGATTTGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaag                 <  1:4019147/68‑1 (MQ=255)

                                                  aTCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaaga                >  1:13003831/1‑68 (MQ=255)

                                                   tCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaaga                >  1:32757896/1‑67 (MQ=255)

                                                   tCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaaga                >  1:42417680/1‑67 (MQ=255)

                                                   tCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaaga                >  1:124422/1‑67 (MQ=255)

                                                    cccGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAac              <  1:152825/68‑1 (MQ=255)

                                                       ggAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac           <  1:42736897/68‑1 (MQ=255)

                                                         agaTTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc         >  1:4797844/1‑68 (MQ=255)

                                                         agaTTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc         >  1:5851647/1‑68 (MQ=255)

                                                         agaTTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc         >  1:38242372/1‑68 (MQ=255)

                                                         agaTTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc         >  1:27699101/1‑68 (MQ=255)

                                                            ttGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAc       <  1:18978985/67‑1 (MQ=255)

                                                            ttGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAc       <  1:27226195/67‑1 (MQ=255)

                                                            ttGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt      <  1:40842272/68‑1 (MQ=255)

                                                            ttGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt      >  1:19765031/1‑68 (MQ=255)

                                                            ttGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt      <  1:9649194/68‑1 (MQ=255)

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

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

                                                             tGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACtt     <  1:11583947/68‑1 (MQ=255)

                                                             tGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACtt     <  1:25889140/68‑1 (MQ=255)

                                                              ggAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTc    <  1:41484709/68‑1 (MQ=255)

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

CTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTG  >  CP009273/4106373‑4106504

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