Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,738	T→G	22.0%	A172A (GCA→GCC) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,738	0	T	G	22.0%	122.4 / 50.1	82	A172A (GCA→GCC) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (38/26);  new base G (10/8);  total (48/34)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACA  >  minE/609668‑609808                                                                       |                                                                       cTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACAGCCACGGCCatat                                                                        >  1:3121457/1‑71 (MQ=255)    tcatcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCa                                                                     >  1:2098277/1‑71 (MQ=255)     catcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCAt                                                                    >  1:3428461/1‑71 (MQ=255)       tcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCt                                                                  <  1:2305776/71‑1 (MQ=255)         aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTgg                                                                >  1:1033102/1‑71 (MQ=255)         aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTgg                                                                >  1:3538684/1‑71 (MQ=255)          ggTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTgg                                                                >  1:1976427/1‑70 (MQ=255)           gTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCa                                                              <  1:2268395/71‑1 (MQ=255)             aGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGa                                                            >  1:878512/1‑71 (MQ=255)              gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATc                                                                   >  1:2225540/1‑63 (MQ=255)              gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATc                                                                   >  1:3429532/1‑63 (MQ=255)              gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGa                                                            <  1:1835672/70‑1 (MQ=255)              gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                                                           >  1:1429495/1‑71 (MQ=255)              gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                                                           >  1:2999190/1‑71 (MQ=255)               ctgATGCCATACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTg                                                                 >  1:941580/3‑64 (MQ=255)               ctgATGCCATACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTg                                                                 >  1:208337/3‑64 (MQ=255)                 gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                                        >  1:772263/1‑71 (MQ=255)                 gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                                        <  1:651234/71‑1 (MQ=255)                 gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                                        >  1:1975028/1‑71 (MQ=255)                 gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCg                                                         <  1:3947610/70‑1 (MQ=255)                   tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                                                           >  1:1150230/1‑66 (MQ=255)                    aCCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                                        >  1:3559219/1‑68 (MQ=255)                    aCCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGt                                                       >  1:2029271/1‑69 (MQ=255)                     ccGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGAATCTGGCAGACCGGTCa                                                     <  1:2138738/70‑1 (MQ=255)                      cGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGa                                                   <  1:614699/71‑1 (MQ=255)                       gTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                                  <  1:2199452/71‑1 (MQ=255)                          cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                                  <  1:1353960/68‑1 (MQ=255)                           aGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGcc                                              <  1:3794751/71‑1 (MQ=255)                              gctaCGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGa                                           <  1:2344681/68‑1 (MQ=255)                              gcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCt                                                                  <  1:1224188/48‑1 (MQ=255)                                 aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACAcc                                              <  1:3513508/65‑1 (MQ=255)                                 aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAttt                                        <  1:3367157/71‑1 (MQ=255)                                  cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa                                       <  1:3746778/71‑1 (MQ=255)                                  cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa                                       >  1:1437964/1‑71 (MQ=255)                                      aCGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACgg                                    >  1:3381941/1‑70 (MQ=255)                                       cGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGc                                   >  1:2347202/1‑70 (MQ=255)                                        ggTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCaa                                 >  1:277157/1‑71 (MQ=255)                                          tAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACg                               <  1:3363619/71‑1 (MQ=255)                                            gTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCg                               >  1:1254477/1‑69 (MQ=255)                                             tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAtt                                         <  1:2711657/58‑1 (MQ=255)                                             tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                                      <  1:3535560/61‑1 (MQ=255)                                             tCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAg                            <  1:2118094/71‑1 (MQ=255)                                               aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCa                                                     >  1:2631320/1‑44 (MQ=255)                                               aCCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAgc                                                   <  1:3941157/46‑2 (MQ=255)                                                ccGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                         <  1:1036060/71‑1 (MQ=255)                                                ccGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGt                                                       <  1:2146723/41‑1 (MQ=255)                                                 cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                            >  1:2603756/1‑67 (MQ=255)                                                 cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                            >  1:762822/1‑67 (MQ=255)                                                 cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg                        >  1:1591625/1‑71 (MQ=255)                                                 cGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGAtt                                         <  1:1316290/54‑1 (MQ=255)                                                  gatgatACGGCCACGGCCATATGGATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                       >  1:3618208/1‑71 (MQ=255)                                                  gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                       >  1:2490979/1‑71 (MQ=255)                                                  gataatGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAtat                       <  1:3795941/71‑3 (MQ=255)                                                   atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTc                                                      >  1:2591318/1‑39 (MQ=255)                                                   atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTc                                                      >  1:1744646/1‑39 (MQ=255)                                                   atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                              >  1:563614/1‑63 (MQ=255)                                                   atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                              >  1:1713932/1‑63 (MQ=255)                                                   atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                      >  1:8738/1‑71 (MQ=255)                                                   atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAAGAtag                       >  1:1777735/1‑68 (MQ=255)                                                      atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCgg                   <  1:2110042/71‑1 (MQ=255)                                                      atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                    <  1:264934/70‑1 (MQ=255)                                                       tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                    >  1:1819363/1‑69 (MQ=255)                                                       tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                    >  1:403082/1‑69 (MQ=255)                                                          ggccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAg                 <  1:3154099/69‑1 (MQ=255)                                                           gccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCa                             <  1:681140/56‑1 (MQ=255)                                                           cccacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAgc                           >  1:2368655/2‑57 (MQ=255)                                                            ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAgtgt              >  1:2180042/1‑70 (MQ=255)                                                            ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAgt                >  1:2812195/1‑68 (MQ=255)                                                             cacggccaTATGCATGTGGCAGACCGGTCAGAACACCAGAtt                                         <  1:1626566/42‑1 (MQ=255)                                                             cacggccaTATGCATCTGTCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa            <  1:5535/71‑1 (MQ=255)                                                             cacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGAtt                                         <  1:3645305/42‑1 (MQ=255)                                                             cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaa            <  1:3006292/71‑1 (MQ=255)                                                              acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa           >  1:1205796/1‑71 (MQ=255)                                                               cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa             >  1:2006983/1‑68 (MQ=255)                                                               cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa             <  1:1759126/68‑1 (MQ=255)                                                               cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa           >  1:2389073/1‑70 (MQ=255)                                                               cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa           >  1:2683014/1‑70 (MQ=255)                                                               cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa           >  1:2393813/1‑70 (MQ=255)                                                               cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa           >  1:3015164/1‑70 (MQ=255)                                                                tgccaGATGCATATGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg         >  1:2894211/2‑71 (MQ=255)                                                                ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg         >  1:96175/1‑71 (MQ=255)                                                                ggccaTATGCATCTGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg         <  1:1154638/71‑1 (MQ=255)                                                                  ccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAAcat        >  1:3254109/1‑68 (MQ=255)                                                                     tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACgg                                    <  1:735503/39‑1 (MQ=255)                                                                     tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                            <  1:1244351/47‑1 (MQ=255)                                                                      atGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                  >  1:2681569/1‑56 (MQ=255)                                                                      atGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                  >  1:3667396/1‑56 (MQ=255)                                                                       tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa           >  1:89121/1‑62 (MQ=255)                                                                       tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca  <  1:3863764/71‑1 (MQ=255)                                                                       tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca  <  1:3713503/71‑1 (MQ=255)                                                                       tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca  <  1:1557305/71‑1 (MQ=255)                                                                       |                                                                       CTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACA  >  minE/609668‑609808

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