Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	709,001	A→G	9.4%	Y102H (TAT→CAT)	yceI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	709,001	0	A	G	9.4%	283.6 / 16.8	96	Y102H (TAT→CAT)	yceI	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/56);  new base G (7/2);  total (38/58)
Fisher's exact test for biased strand distribution p-value = 2.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AGATCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGACGCTG  >  minE/708931‑709071                                                                       |                                                                       aGATCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGata                                                                        >  1:4210891/1‑71 (MQ=255)   aTCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAtt                                                                      <  1:856362/71‑1 (MQ=255)   aTCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAtt                                                                      <  1:1424643/71‑1 (MQ=255)    tCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAttt                                                                     >  1:4747887/1‑71 (MQ=255)    tCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAttt                                                                     >  1:1386783/1‑71 (MQ=255)             atatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATGTTTTGCGGTAtt                                                            <  1:812136/71‑1 (MQ=255)             atatCCAGTTCGGCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            <  1:3623531/71‑1 (MQ=255)              tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAttt                                                                     <  1:1325690/61‑1 (MQ=255)              tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTg                                                                   <  1:2690065/63‑1 (MQ=255)              tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTg                                                                   <  1:3753033/63‑1 (MQ=255)              tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTg                                                                   <  1:3214466/63‑1 (MQ=255)              tatCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTg                                                           <  1:1668441/71‑1 (MQ=255)                tCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCg                                                                 >  1:1753245/1‑63 (MQ=255)                tCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGt                                                               <  1:1820984/65‑1 (MQ=255)                tCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGt                                                               <  1:3858419/65‑1 (MQ=255)                tCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGa                                                          <  1:1024008/70‑1 (MQ=255)                tCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAg                                                         <  1:3721345/71‑1 (MQ=255)                 ccAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAg                                                         >  1:3140962/1‑70 (MQ=255)                  cAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAgg                                                        <  1:2629709/70‑1 (MQ=255)                   aGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaa                                                      <  1:2385814/71‑1 (MQ=255)                   aGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaa                                                      <  1:4196178/71‑1 (MQ=255)                   aGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaa                                                      <  1:4195083/71‑1 (MQ=255)                   aGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaa                                                      <  1:3119284/71‑1 (MQ=255)                    gTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            <  1:2741610/64‑1 (MQ=255)                    gTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            <  1:323022/64‑1 (MQ=255)                    gTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            <  1:2058575/64‑1 (MQ=255)                    gTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGaaa                                                     <  1:71023/71‑1 (MQ=255)                       cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATAtttt                                                                    >  1:2710125/1‑53 (MQ=255)                       cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAg                                                         >  1:2919258/1‑64 (MQ=255)                       cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAg                                                         >  1:3136263/1‑64 (MQ=255)                       cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAg                                                         >  1:190927/1‑64 (MQ=255)                       cGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCt                                                  >  1:269371/1‑71 (MQ=255)                        gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 >  1:2969289/1‑71 (MQ=255)                        gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:537186/71‑1 (MQ=255)                        gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:2665645/71‑1 (MQ=255)                        gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:4188940/71‑1 (MQ=255)                        gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:1695755/71‑1 (MQ=255)                        gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 >  1:2274470/1‑71 (MQ=255)                        gTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:2154732/71‑1 (MQ=255)                         tCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGc                                                <  1:594416/71‑1 (MQ=255)                          cACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCa                                               >  1:3960237/1‑71 (MQ=255)                           aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGa                                                          <  1:2257355/59‑1 (MQ=255)                           aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGa                                                          <  1:1304137/59‑1 (MQ=255)                           aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGa                                                          <  1:2070267/59‑1 (MQ=255)                           aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCAc                                              >  1:1701364/1‑71 (MQ=255)                           aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCAc                                              >  1:4747249/1‑71 (MQ=255)                           aCCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCAc                                              >  1:3566381/1‑71 (MQ=255)                             cGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGt                                                               <  1:602545/52‑1 (MQ=255)                               tctttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          <  1:3381919/71‑1 (MQ=255)                               tctttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          <  1:3599000/71‑1 (MQ=255)                                ctttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          <  1:3544634/70‑1 (MQ=255)                                 tttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:2768733/62‑1 (MQ=255)                                 tttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:3852228/62‑1 (MQ=255)                                 tttcttCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTg                                                 <  1:3863720/62‑1 (MQ=255)                                   tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            >  1:650069/1‑49 (MQ=255)                                   tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            >  1:4561198/1‑49 (MQ=255)                                   tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTAtt                                                            >  1:766149/1‑49 (MQ=255)                                   tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTg                                            >  1:1344090/1‑65 (MQ=255)                                   tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          >  1:877286/1‑67 (MQ=255)                                   tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          >  1:2842289/1‑67 (MQ=255)                                   tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCg                                          >  1:4643470/1‑67 (MQ=255)                                   tcttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGaa                                        <  1:1540179/69‑1 (MQ=255)                                     ttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAAt                                                    >  1:643488/1‑55 (MQ=255)                                     ttcACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAAt                                                    >  1:3804013/1‑55 (MQ=255)                                       cACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATg                                    <  1:3591480/69‑1 (MQ=255)                                        aCGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGttt                                 <  1:963352/71‑1 (MQ=255)                                        aCGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGttt                                 <  1:2468774/71‑1 (MQ=255)                                        aCGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGttt                                 <  1:11483/71‑1 (MQ=255)                                        aCGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGttt                                 <  1:1793961/71‑1 (MQ=255)                                           cTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTAt                               <  1:873697/70‑1 (MQ=255)                                                  aGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCg                       <  1:3627036/71‑1 (MQ=255)                                                  aGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCg                       <  1:1389041/71‑1 (MQ=255)                                                   tgTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCg                       <  1:881780/69‑1 (MQ=255)                                                   ggTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCgg                      >  1:3021797/1‑71 (MQ=255)                                                      gAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTAt                               >  1:3441264/1‑59 (MQ=255)                                                      gAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTAt                               >  1:55913/1‑59 (MQ=255)                                                      gAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTAt                               >  1:2713551/1‑59 (MQ=255)                                                      gAATGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTAt                               >  1:2464745/1‑59 (MQ=255)                                                      gAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCgg                      >  1:1692551/1‑68 (MQ=255)                                                      gAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCgg                      >  1:4597172/1‑68 (MQ=255)                                                      gAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCgg                      >  1:2630899/1‑68 (MQ=255)                                                       aaTGTTGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTg                  <  1:2010665/71‑1 (MQ=255)                                                       aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   <  1:1122831/70‑1 (MQ=255)                                                       aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   <  1:3878876/70‑1 (MQ=255)                                                       aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   <  1:2640879/70‑1 (MQ=255)                                                       aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   <  1:4458898/70‑1 (MQ=255)                                                       aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   >  1:1448243/1‑70 (MQ=255)                                                       aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTg                  <  1:999372/71‑1 (MQ=255)                                                       aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTg                  <  1:3731636/71‑1 (MQ=255)                                                       aaTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTg                  <  1:1366492/71‑1 (MQ=255)                                                        aTGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGa                 >  1:867164/1‑71 (MQ=255)                                                         tGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGt                                   <  1:2044292/52‑1 (MQ=255)                                                           ttGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGa                 >  1:2420128/1‑68 (MQ=255)                                                           ttGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGAt                >  1:2293915/1‑69 (MQ=255)                                                           ttGCCTGTGGATGTTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGAt                >  1:951578/1‑69 (MQ=255)                                                           ttGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGt                   <  1:2280829/66‑1 (MQ=255)                                                           ttGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGAt                >  1:769369/1‑69 (MQ=255)                                                                     atatTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGACGc    >  1:2360724/1‑71 (MQ=255)                                                                       atTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGACGCtg  >  1:3213359/1‑71 (MQ=255)                                                                       atTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGACGCtg  >  1:3914724/1‑71 (MQ=255)                                                                       |                                                                       AGATCGCCGGTAATATCCAGTTCGTCACCGTCTTTCTTCACGCTGGTGGAGGTGAATGTTGCCTGTGGATATTTTGCGGTATTGAGGAAATCTGCACTGCGAAGATGTTTATCGCGTTCGGCGTGATTAGTATCGACGCTG  >  minE/708931‑709071

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