BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,323,988	A→G	42.8%	G122G (GGT→GGC)	pfkA ←	6‑phosphofructokinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,323,988	0	A	G	42.8%	38.6 / 111.0	140	G122G (GGT→GGC)	pfkA	6‑phosphofructokinase I
Reads supporting (aligned to +/- strand): ref base A (15/65); new base G (9/51); total (24/116)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 140 total aligned reads displayed.

CGCAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGAACCGTC  >  minE/2323925‑2324044

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:1830529/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:1526523/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:1627848/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:2680867/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:2620358/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:356830/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:1758871/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:1774343/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:1197555/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:3431107/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:1887719/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:2492382/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:2965222/69‑1 (MQ=255)

cgcAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:225602/69‑1 (MQ=255)

            aCCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGcc                                         <  1:1232839/69‑1 (MQ=255)

            aCCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGcc                                         <  1:1969669/69‑1 (MQ=255)

            aCCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGcc                                         <  1:1442287/69‑1 (MQ=255)

            aCCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGcc                                         <  1:1492874/69‑1 (MQ=255)

            aCCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGcc                                         <  1:1036516/69‑1 (MQ=255)

            aCCGATAGTGTAGTCAGTCCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGcc                                         <  1:2600713/69‑1 (MQ=255)

             ccGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGGCCGAt                                                     <  1:1679263/56‑1 (MQ=255)

                        gTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCa                              <  1:2737151/68‑1 (MQ=255)

                        gTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCa                              <  1:455944/68‑1 (MQ=255)

                        gTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCa                              <  1:1948146/68‑1 (MQ=255)

                        gTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAg                             >  1:863715/1‑69 (MQ=255)

                          cAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGa                            <  1:2065172/68‑1 (MQ=255)

                          cAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGa                            <  1:1391987/68‑1 (MQ=255)

                          cAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGa                            <  1:1378567/68‑1 (MQ=255)

                          cAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGa                            <  1:1149629/68‑1 (MQ=255)

                          cAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGa                            <  1:1048312/68‑1 (MQ=255)

                          cAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGa                            <  1:3104375/68‑1 (MQ=255)

                           aGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGt                                >  1:1271729/1‑63 (MQ=255)

                           aGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGt                                >  1:2022765/1‑63 (MQ=255)

                           aGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACg                          <  1:214949/69‑1 (MQ=255)

                           aGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACg                          <  1:2052233/69‑1 (MQ=255)

                           aGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACg                          <  1:2983608/69‑1 (MQ=255)

                           aGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACg                          <  1:783791/69‑1 (MQ=255)

                           aGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACg                          <  1:2646296/69‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGTAAGcc                                         <  1:2797048/52‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGc                         <  1:3329063/68‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGc                         <  1:2509577/68‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGc                         <  1:2330758/68‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCa                        <  1:2566018/69‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCa                        <  1:174382/69‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCa                        <  1:369072/69‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCa                        <  1:457259/69‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCa                        <  1:492319/69‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCa                        <  1:894999/69‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCa                        <  1:906295/69‑1 (MQ=255)

                             tGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCAATTTCGGTCAGACGc                         <  1:1021370/68‑1 (MQ=255)

                                 tttGATGTCGTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTg                     <  1:2805684/68‑1 (MQ=255)

                                 tttGATGTCGTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTg                     <  1:2958014/68‑1 (MQ=255)

                                 tttGATGTCGTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTg                     <  1:696612/68‑1 (MQ=255)

                                 tttGATGTCGTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTg                     <  1:1021697/68‑1 (MQ=255)

                                 tttGATGTCGTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTg                     <  1:2153687/68‑1 (MQ=255)

                                 tttGATGTCGTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTg                     <  1:2585933/68‑1 (MQ=255)

                                 tttGATGTCGTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTg                     <  1:1481015/68‑1 (MQ=255)

                                 tttGATGTCGTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTg                     <  1:2032772/68‑1 (MQ=255)

                                 tttGATGTCGTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTg                     <  1:2097454/68‑1 (MQ=255)

                                 tttGATGTCGTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTg                     <  1:2272391/68‑1 (MQ=255)

                                 tttGATGTCGTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTg                     <  1:2477968/68‑1 (MQ=255)

                                      tGTCGTTGGCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCAgcc                <  1:15261/68‑1 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCTATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCa               >  1:491103/1‑67 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGaa                                            >  1:3122742/1‑38 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGCCGCATTGCACCCa               >  1:678000/1‑67 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCa               >  1:1086626/1‑67 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCa               >  1:906299/1‑67 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCa               >  1:549846/1‑67 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCa               >  1:3261724/1‑67 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCa               >  1:2536878/1‑67 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATg             <  1:1159536/69‑1 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATg             <  1:805878/69‑1 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATg             <  1:702928/69‑1 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATg             <  1:1629942/69‑1 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATg             <  1:956716/69‑1 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATg             <  1:19338/69‑1 (MQ=255)

                                        tCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATg             <  1:1970027/69‑1 (MQ=255)

                                          gTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGaa                                            >  1:663968/1‑36 (MQ=255)

                                          gTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGt            >  1:1652146/1‑68 (MQ=255)

                                          gTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGt            >  1:2859481/1‑68 (MQ=255)

                                          gTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGt            >  1:2245098/1‑68 (MQ=255)

                                          gTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGt            >  1:2978205/1‑68 (MQ=255)

                                          gTTGTCGATAGTGCCCGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGt            >  1:3306394/1‑68 (MQ=255)

                                          gTTGTCGATAGTGCACGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGt            >  1:318941/1‑68 (MQ=255)

                                                   aGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGAACCGTc  >  1:1432796/1‑69 (MQ=255)

                                                   aGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGAACCGTc  >  1:2177008/1‑69 (MQ=255)

                                                   aGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGAACCGTc  >  1:747737/1‑69 (MQ=255)

                                                       cccGGCAGGCCTATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:1048199/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCTCCCATGTAGGaa       <  1:3345674/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:2395814/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:229515/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:1856244/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:959288/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:1373939/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:745018/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:2946056/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:689388/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:597672/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:561601/60‑1 (MQ=255)

                                                       cccGGCAGGCCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGaa       <  1:1611773/60‑1 (MQ=255)

CGCAGTGAAGAAACCGATAGTGTAGTCAGTGCCTTTGATGTCGTTGTCGATAGTGCCCGGCAGACCGATGCACGGGAAGCCCATTTCGGTCAGACGCATTGCACCCATGTAGGAACCGTC  >  minE/2323925‑2324044

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`