Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,661	G→C	6.6%	S531S (TCC→TCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,661	0	G	C	6.6%	272.7 / 6.6	91	S531S (TCC→TCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (46/39);  new base C (2/4);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTT  >  minE/608591‑608731                                                                       |                                                                       ccAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACg                                                                        <  1:827431/71‑1 (MQ=255)   aGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGa                                                                      >  1:4964313/1‑71 (MQ=255)    gACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGa                                                                      >  1:2606141/1‑70 (MQ=255)      ccGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAg                                                                   <  1:3184799/71‑1 (MQ=255)       cGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACa                                                                    >  1:5132053/1‑69 (MQ=255)       cGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGa                                                                  <  1:4881319/71‑1 (MQ=255)           ttCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAg                                                              >  1:2221473/1‑71 (MQ=255)           ttCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAg                                                              >  1:5559560/1‑71 (MQ=255)             cGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                             <  1:829012/70‑1 (MQ=255)             cGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                             <  1:5559538/70‑1 (MQ=255)             cGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                             <  1:525254/70‑1 (MQ=255)              gTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGc                                                             <  1:5654942/69‑1 (MQ=255)               tCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGAt                                                          <  1:4089282/71‑1 (MQ=255)                 aCGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATaa                                                        >  1:4233397/1‑70 (MQ=255)                 aCGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATAc                                                        >  1:1333560/1‑71 (MQ=255)                   gAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACg                                                                        <  1:914416/53‑1 (MQ=255)                   gAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACca                                                      <  1:463864/71‑1 (MQ=255)                   gAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACca                                                      <  1:2480866/71‑1 (MQ=255)                   gAATCGATTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGAc                                                                     <  1:4183711/56‑1 (MQ=255)                     aTCGGTTTAACTTTCGCATATTTTATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACcaca                                                    <  1:3384186/71‑1 (MQ=255)                     aTCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACcaca                                                    <  1:2769197/71‑1 (MQ=255)                      tCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCacac                                                   >  1:5277265/1‑71 (MQ=255)                         gTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGcc                                                <  1:1521978/71‑1 (MQ=255)                         gTTAAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGc                                                 <  1:3387994/70‑1 (MQ=255)                           ttAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCAt                                              >  1:5109093/1‑71 (MQ=255)                            tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                             >  1:5075723/1‑71 (MQ=255)                            tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                             >  1:115620/1‑71 (MQ=255)                            tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                             >  1:3085328/1‑71 (MQ=255)                             aaCTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                                             >  1:2842431/1‑70 (MQ=255)                               cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                                               <  1:4423978/50‑1 (MQ=255)                               cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                                               <  1:4223321/50‑1 (MQ=255)                               cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                                               <  1:758341/50‑1 (MQ=255)                               cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                                               <  1:38620/50‑1 (MQ=255)                               cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTg                                          <  1:117589/71‑1 (MQ=255)                                tttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                                         >  1:2713393/1‑71 (MQ=255)                                 ttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTg                                          >  1:5359604/1‑69 (MQ=255)                                 ttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCg                                        >  1:4107728/1‑71 (MQ=255)                                  tCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                                         <  1:2585078/69‑1 (MQ=255)                                  tCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                                         <  1:2390300/69‑1 (MQ=255)                                  gggCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCga                                       <  1:373464/69‑2 (MQ=255)                                   cGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATgg                                            >  1:1881742/1‑65 (MQ=255)                                   cGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATgg                                            >  1:1824049/1‑65 (MQ=255)                                    gCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGAt                                     <  1:937516/71‑1 (MQ=255)                                     cttatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATa                                    <  1:2431491/69‑1 (MQ=255)                                     cATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATa                                    <  1:920176/71‑1 (MQ=255)                                     cATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGa                                                           >  1:1089330/1‑48 (MQ=255)                                       tatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACa                                                                    <  1:3466546/37‑1 (MQ=255)                                       tatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACa                                                                    <  1:4210037/37‑1 (MQ=255)                                       tatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGAt                                                          <  1:1127564/47‑1 (MQ=255)                                       tatTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGcccg                                                              <  1:1613128/43‑3 (MQ=12)                                        atTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                             <  1:574690/43‑1 (MQ=255)                                        atTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACg                                 >  1:1149481/1‑71 (MQ=255)                                         tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAg                                                                   >  1:1767876/1‑36 (MQ=255)                                         tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAg                                                                   >  1:3417255/1‑36 (MQ=255)                                         tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGa                                                           <  1:218253/44‑1 (MQ=255)                                         tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                                >  1:4624953/1‑71 (MQ=255)                                         tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                                >  1:3464094/1‑71 (MQ=255)                                         tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                                >  1:4120457/1‑71 (MQ=255)                                         tttGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGc                                                             >  1:1975401/1‑42 (MQ=255)                                          ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                                >  1:3549713/1‑70 (MQ=255)                                          ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                                >  1:4482926/1‑70 (MQ=255)                                           tGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTc                               <  1:4495706/70‑1 (MQ=255)                                           tGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACACCATGGTGCGGATAACGTCa                              >  1:1826472/1‑71 (MQ=255)                                              ttGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACgg                           <  1:709857/71‑1 (MQ=255)                                               tGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGt                          >  1:3881570/1‑71 (MQ=255)                                                gCAGACAGGGAGTCAGCAGCGACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                         >  1:1409918/1‑71 (MQ=255)                                                gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATaa                                   >  1:2960014/1‑61 (MQ=255)                                                gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATaa                                   >  1:5666001/1‑61 (MQ=255)                                                gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                         >  1:4431275/1‑71 (MQ=255)                                                gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                         >  1:5365467/1‑71 (MQ=255)                                                gCAGACAGGGAGTCAGCAGCAACGGACAGAACAGCGATACCACACGCCATGGTGCGGATaa                                   >  1:898897/1‑61 (MQ=255)                                                 cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                         <  1:70206/70‑1 (MQ=255)                                                 cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCg                        >  1:3893114/1‑71 (MQ=255)                                                   gACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCgcg                      <  1:994016/71‑3 (MQ=255)                                                   gACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTg                      <  1:2664318/71‑1 (MQ=255)                                                    aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTg                      >  1:4380746/1‑70 (MQ=255)                                                    aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGc                     <  1:259841/71‑1 (MQ=255)                                                     cAGGGAGTCAGCAGCACCGGACAGACCAGCGATACcaca                                                    >  1:4979632/1‑39 (MQ=255)                                                     cAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGc                     <  1:5063168/70‑1 (MQ=255)                                                     cAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCa                    <  1:2870297/71‑1 (MQ=255)                                                      aGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAg                   >  1:4929491/1‑71 (MQ=255)                                                        ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCa                    >  1:4325930/1‑68 (MQ=255)                                                        ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCa                    >  1:767944/1‑68 (MQ=255)                                                            tCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt             <  1:5488704/71‑1 (MQ=255)                                                              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAg                   >  1:4966893/1‑63 (MQ=255)                                                              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt             >  1:2339797/1‑69 (MQ=255)                                                              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt             >  1:3859894/1‑69 (MQ=255)                                                              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCa           <  1:5218487/71‑1 (MQ=255)                                                                   caACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAAc                                  <  1:1283199/43‑1 (MQ=255)                                                                    aaCGGACAGACCAGCGATACCACACGCCATGGTGCGGATAAc                                  <  1:1612146/42‑1 (MQ=255)                                                                     aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagaga       >  1:606469/1‑68 (MQ=255)                                                                     aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagaga       >  1:5597367/1‑68 (MQ=255)                                                                     aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGc    >  1:225641/1‑71 (MQ=255)                                                                       ggACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt  <  1:203083/71‑1 (MQ=255)                                                                       |                                                                       CCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTT  >  minE/608591‑608731

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