BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,744	T→C	7.4%	Y136C (TAC→TGC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,744	0	T	C	7.4%	2245.9 / 104.6	660	Y136C (TAC→TGC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base T (276/333); new base C (22/27); total (299/361)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 679 total aligned reads displayed.

ACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATG  >  CP009273/4106677‑4106808

acgacgTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGt                                                                  >  1:27370357/1‑68 (MQ=255)

acgacgTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGt                                                                  <  1:5842065/68‑1 (MQ=255)

acgacgTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGt                                                                  >  1:20600797/1‑68 (MQ=255)

acgacgTGCACGCTCGCGAGAGCCTTCCACATCGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGt                                                                  >  1:11016521/1‑68 (MQ=255)

    cgTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                                                              >  1:31346887/1‑68 (MQ=255)

     gTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACggg                                                             <  1:15312457/68‑1 (MQ=255)

       gCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGt                                                            <  1:33856699/67‑1 (MQ=255)

       gCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTc                                                           <  1:21863141/68‑1 (MQ=255)

       gCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTc                                                           <  1:10823142/68‑1 (MQ=255)

         aCGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCaa                                                         <  1:1550012/68‑1 (MQ=255)

          cGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCaa                                                         >  1:29719334/1‑67 (MQ=255)

          cGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGCACGGGTCaa                                                         >  1:3791723/1‑67 (MQ=255)

             tCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCa                                                      <  1:2382617/67‑1 (MQ=255)

             tCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCa                                                      <  1:17598469/67‑1 (MQ=255)

              cgcgAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGCACGGGTCAATCAcc                                                    >  1:22407931/1‑68 (MQ=255)

                 gagaGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGa                                                 <  1:6975200/68‑1 (MQ=255)

                 gagaGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGa                                                 <  1:1685886/68‑1 (MQ=255)

                 gagaGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGa                                                 <  1:17996949/68‑1 (MQ=255)

                  agagCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGa                                                 >  1:3911821/1‑67 (MQ=255)

                   gagCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGAcc                                               >  1:32451927/1‑68 (MQ=255)

                     gCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGAcc                                               <  1:30036105/66‑1 (MQ=255)

                     gCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCgg                                             >  1:31902387/1‑68 (MQ=255)

                     gCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCgg                                             >  1:16438897/1‑68 (MQ=255)

                     gCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCg                                              <  1:9705345/67‑1 (MQ=255)

                     gCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCg                                              <  1:4578790/67‑1 (MQ=255)

                      ccTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGt                                            >  1:11872905/1‑68 (MQ=255)

                       cTTCCACATGGTCGAGTATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTa                                           <  1:1454737/68‑1 (MQ=255)

                       cTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGt                                            <  1:18626345/67‑1 (MQ=255)

                       cTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGt                                            <  1:12540787/67‑1 (MQ=255)

                       cTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGCACGGGTCAATCACCAGACCGGt                                            <  1:2438399/67‑1 (MQ=255)

                           cacaTGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgc                                       <  1:6236783/68‑1 (MQ=255)

                            acaTGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTAtt                                         >  1:24127874/1‑65 (MQ=255)

                              aTGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctgc                                    >  1:26909525/1‑68 (MQ=255)

                              aTGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctgc                                    >  1:17403588/1‑68 (MQ=255)

                              aTGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctgc                                    >  1:27974672/1‑68 (MQ=255)

                              aTGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctgc                                    >  1:28993678/1‑68 (MQ=255)

                              aTGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctgc                                    >  1:12196184/1‑68 (MQ=255)

                              aTGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctgc                                    >  1:32716141/1‑68 (MQ=255)

                              aTGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctgc                                    >  1:16014978/1‑68 (MQ=255)

                              aTGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgctgc                                    >  1:17169116/1‑68 (MQ=255)

                               tGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCg                                   >  1:13582294/1‑68 (MQ=255)

                               tGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCg                                   >  1:23235324/1‑68 (MQ=255)

                               tGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCg                                   >  1:1964110/1‑68 (MQ=255)

                                ggTCGAGGATCCACTTCACTTTGGTGCCTGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCgg                                  <  1:15154082/68‑1 (MQ=255)

                                ggTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCgg                                  <  1:29513089/68‑1 (MQ=255)

                                ggTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGCACGGGGCAATCACCAGACCGGTATTGCTGCgg                                  <  1:21613931/68‑1 (MQ=255)

                                 gTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGa                                 >  1:19271965/1‑68 (MQ=255)

                                  tCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGa                                 <  1:7217891/67‑1 (MQ=255)

                                  tCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGa                                 <  1:7960029/67‑1 (MQ=255)

                                  tCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGa                                 <  1:25635867/67‑1 (MQ=255)

                                   cGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGata                               <  1:25225982/68‑1 (MQ=255)

                                      ggATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGCCCGGTATTGCTGCGGATATaa                            >  1:9458157/1‑68 (MQ=255)

                                      ggATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATaa                            >  1:31663729/1‑68 (MQ=255)

                                       gATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATaa                            >  1:19990898/1‑67 (MQ=255)

                                        aTCCACTTCACTTTGGTGCCAGAAAAGTCCGGGTCAATCACCAGACCGGTATTGCTGCGGATATaa                            >  1:497905/1‑66 (MQ=255)

                                        aTCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATaa                            >  1:27729644/1‑66 (MQ=255)

                                        aTCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATaa                            >  1:33640905/1‑66 (MQ=255)

                                        aTCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATaa                            <  1:34283719/66‑1 (MQ=255)

                                        aTCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATaa                            >  1:22742212/1‑66 (MQ=255)

                                        aTCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtc                          >  1:8318492/1‑68 (MQ=255)

                                        aTCCACTTCACTTTGGTGCCAGAAAAGCACGGGTCAATCACCAGACCGGTATTGCTGCGGATATaa                            >  1:28588184/1‑66 (MQ=255)

                                        aTCCACTTCACTTTGGTGACAGAAAAGTACGGGACAATCACCAGACCGGTATTTATGCGGATATaa                            >  1:4315687/1‑66 (MQ=255)

                                        aTCCACTACACTTTGGTCCCACAAAACTACGGGTCAATCACCAGACCGGTATTGCACCGGATATaa                            >  1:25430261/1‑66 (MQ=255)

                                        aTCAACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATaa                            >  1:14307842/1‑66 (MQ=255)

                                         tCCACTTCACTTTGGTGGCAGAAAAGTACGGGTCAATCACCAGTCCGGTATTGCTGCGGATATAAtct                         <  1:6786268/68‑1 (MQ=255)

                                         tCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtct                         <  1:18215537/68‑1 (MQ=255)

                                         tCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtct                         <  1:6481197/68‑1 (MQ=255)

                                         tCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtct                         <  1:32274016/68‑1 (MQ=255)

                                         tCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtct                         <  1:2046984/68‑1 (MQ=255)

                                         tCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtct                         <  1:26161783/68‑1 (MQ=255)

                                         tCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtct                         <  1:31203841/68‑1 (MQ=255)

                                         tCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtct                         <  1:28301072/68‑1 (MQ=255)

                                         tCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtct                         <  1:30665555/68‑1 (MQ=255)

                                         tCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtct                         <  1:2088492/68‑1 (MQ=255)

                                         tCCACTTCACTTTGGTGCCAGAAAAGCACGGGTCAATCACCAGGCCGGTATTGCTGCGGATATAAtct                         <  1:13106837/68‑1 (MQ=255)

                                         tCCACTTCACTTTGGTGCCAGAAAAGCACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtct                         <  1:523112/68‑1 (MQ=255)

                                          ccACTTCACTTTGGTGCCAGAAAAGCACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAAtctt                        <  1:22983464/68‑1 (MQ=255)

                                             cttcacttTGGTGCCAGAATAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTa                     >  1:29337349/1‑68 (MQ=255)

                                               tcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTaa                    >  1:10200767/1‑67 (MQ=255)

                                                 acttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAAcc                 >  1:754421/1‑68 (MQ=255)

                                                 acttTGGTGCCAGAAAAGCACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAAcc                 >  1:9238018/1‑68 (MQ=255)

                                                    ttGGTGCCAGAAAAGTACGGGTCAATCTCCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTc              >  1:18994181/1‑68 (MQ=255)

                                                    ttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTc              >  1:34437980/1‑68 (MQ=255)

                                                    ttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTc              >  1:3455765/1‑68 (MQ=255)

                                                    ttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTc              >  1:14007653/1‑68 (MQ=255)

                                                        tGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACg           <  1:26312938/67‑1 (MQ=255)

                                                        tGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACg           <  1:33106049/67‑1 (MQ=255)

                                                        tGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACg           <  1:11526920/67‑1 (MQ=255)

                                                        tGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACg           <  1:1288173/67‑1 (MQ=255)

                                                        tGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACg           <  1:24749240/67‑1 (MQ=255)

                                                        tGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTCCTGCGGATATAATCTTCTAAACCGTCACGt          >  1:20211886/1‑68 (MQ=255)

                                                        tGCCAGAAAAGTACGGGTCAATCACCAGACCGGTAATGCTGCGGATATAATCTTCTAAACCGTCACg           <  1:34075658/67‑1 (MQ=255)

                                                         gCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGtt         <  1:5636049/68‑1 (MQ=255)

                                                         gCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGtt         <  1:14852424/68‑1 (MQ=255)

                                                         gCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGt          <  1:26654288/67‑1 (MQ=255)

                                                         gCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGt          <  1:1122787/67‑1 (MQ=255)

                                                         gCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGt          <  1:31034949/67‑1 (MQ=255)

                                                              aaaaGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTaaa    >  1:23506823/1‑68 (MQ=255)

                                                              aaaaGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTaaa    >  1:22202361/1‑68 (MQ=255)

                                                                 aGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAAt   >  1:10544721/1‑66 (MQ=255)

ACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATG  >  CP009273/4106677‑4106808

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