Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,279,216	G→C	20.2%	G332A (GGT→GCT) ‡	yeeJ →	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,279,216	0	G	C	20.2%	175.3 / 37.7	84	G332A (GGT→GCT) ‡	yeeJ	adhesin
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base C (7/10);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

GAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAAT  >  minE/1279146‑1279286                                                                       |                                                                       gAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCgg                                                                        >  1:1382062/1‑71 (MQ=255)  aaGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGt                                                                       <  1:2814543/71‑1 (MQ=255)  aaGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGt                                                                       <  1:4941040/71‑1 (MQ=255)  aaGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGACGCACCTTGGCgg                                                                        >  1:1552188/1‑70 (MQ=255)   aGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTa                                                                      <  1:3060332/71‑1 (MQ=255)    gCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTaa                                                                     <  1:2677523/71‑1 (MQ=255)      aCGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAAc                                                                   >  1:206194/1‑71 (MQ=255)       cGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACt                                                                  >  1:124550/1‑71 (MQ=255)         cccGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTgg                                                                >  1:851641/1‑71 (MQ=255)              ccAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGt                                                               <  1:1631047/67‑1 (MQ=255)               cAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                          <  1:3226294/71‑1 (MQ=255)               cAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                          <  1:2388528/71‑1 (MQ=255)                aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                          >  1:1558620/1‑70 (MQ=255)                aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                          >  1:2544932/1‑70 (MQ=255)                aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                          >  1:5107887/1‑70 (MQ=255)                aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                          >  1:5648913/1‑70 (MQ=255)                aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                          >  1:5592350/1‑70 (MQ=255)                aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                          >  1:4022573/1‑70 (MQ=255)                 aTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGaa                                                        <  1:1123278/71‑1 (MQ=255)                 aTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGaa                                                        <  1:1320156/71‑1 (MQ=255)                  tGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGaa                                                        >  1:653716/1‑70 (MQ=255)                  tGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAAc                                                       <  1:3253910/71‑1 (MQ=255)                   ggCTGGGATTTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGaa                                                        >  1:670710/1‑69 (MQ=255)                      tGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGta                                                   >  1:1678832/1‑71 (MQ=255)                        ggATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtat                                                  >  1:4901693/1‑70 (MQ=255)                          aTGTACACGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtatta                                                >  1:5159/1‑70 (MQ=255)                           tGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCt                                                             <  1:5445219/56‑1 (MQ=255)                           tGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCt                                                             <  1:4937695/56‑1 (MQ=255)                           tGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATg                                              >  1:1143794/1‑71 (MQ=255)                           tGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATg                                              <  1:3605568/71‑1 (MQ=255)                             tACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGa                                                         <  1:749481/58‑1 (MQ=255)                             tACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGa                                                         <  1:3076240/58‑1 (MQ=255)                             tACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGc                                            <  1:3654211/71‑1 (MQ=255)                               cgcgCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGa                                          >  1:4276887/1‑71 (MQ=255)                                 cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                                        <  1:908132/71‑1 (MQ=255)                                 cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                                        <  1:3825451/71‑1 (MQ=255)                                 cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                                        <  1:173517/71‑1 (MQ=255)                                 cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                                        <  1:1738861/71‑1 (MQ=255)                                  gcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGa                                       <  1:1164305/71‑1 (MQ=255)                                  gcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGa                                       <  1:4741682/71‑1 (MQ=255)                                      aaaGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtatta                                                >  1:3652558/1‑58 (MQ=255)                                       aaGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGaa                                                        >  1:4721063/1‑49 (MQ=255)                                        aGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGtatta                                                >  1:2689715/1‑56 (MQ=255)                                        aGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGc                                 >  1:215755/1‑71 (MQ=255)                                         gCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGcc                                >  1:4672137/1‑71 (MQ=255)                                           tGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCt                              >  1:5298725/1‑71 (MQ=255)                                             tttACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGc                                            <  1:916949/53‑1 (MQ=255)                                             tttACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATg                                        <  1:904699/57‑1 (MQ=255)                                              ttACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGtatta                                                >  1:4626789/2‑50 (MQ=255)                                              ttACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCg                                           <  1:3885847/54‑1 (MQ=255)                                              ttACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCg                                           <  1:4802785/54‑1 (MQ=255)                                               tACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACa                                                      >  1:1987958/1‑43 (MQ=12)                                               tACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtt                           <  1:5211454/70‑1 (MQ=255)                                               tACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                          >  1:789441/1‑71 (MQ=255)                                                aCCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGAt                                         >  1:3713791/1‑55 (MQ=255)                                                aCCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                          >  1:276759/1‑70 (MQ=255)                                                aCCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                          >  1:3480601/1‑70 (MQ=255)                                                  ccGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtttg                         >  1:930567/1‑67 (MQ=255)                                                  ccGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGt                                                    <  1:3013033/42‑1 (MQ=255)                                                  ccGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGt                                                    <  1:2601327/42‑1 (MQ=255)                                                  ccGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGt                                                    <  1:4177020/42‑1 (MQ=255)                                                   cGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGtatta                                                >  1:1030782/1‑45 (MQ=12)                                                   cGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATg                                        <  1:5291644/53‑1 (MQ=255)                                                   cGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATg                                        <  1:1411955/53‑1 (MQ=255)                                                   cGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtttgg                        >  1:582751/1‑66 (MQ=255)                                                   cGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGAt                       >  1:399332/1‑70 (MQ=255)                                                    gCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaa                     <  1:4646674/71‑1 (MQ=255)                                                     ccTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                                        <  1:3169322/51‑1 (MQ=255)                                                     ccTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                                        <  1:3888036/51‑1 (MQ=255)                                                     ccTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATaaa                    <  1:3471370/71‑1 (MQ=255)                                                        ggCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGAc                 >  1:622087/1‑71 (MQ=255)                                                         gCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaacat                 <  1:2447591/70‑4 (MQ=255)                                                         gCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGa                  <  1:5044865/69‑1 (MQ=255)                                                          ccGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa               <  1:212857/71‑1 (MQ=255)                                                            gCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGaa                                      >  1:3186978/1‑46 (MQ=16)                                                            gCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATc             >  1:5469915/1‑71 (MQ=255)                                                            gCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATc             <  1:1670750/71‑1 (MQ=255)                                                              aCCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCgg           >  1:3774442/1‑71 (MQ=255)                                                               gcTTGGCGGTAAAGTTGTCTATGATCAGTATTATGGCGATg                                        <  1:446047/40‑1 (MQ=25)                                                               ccTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCgg           >  1:1661023/1‑70 (MQ=255)                                                                cTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGccc                               >  1:2790648/1‑49 (MQ=255)                                                                cTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATa                      <  1:5462014/58‑1 (MQ=255)                                                                cTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATa                      <  1:4725793/58‑1 (MQ=255)                                                                    gCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAg      <  1:5163512/70‑1 (MQ=255)                                                                    gCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAg      <  1:4415364/70‑1 (MQ=255)                                                                     cGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTa    <  1:4158708/71‑1 (MQ=255)                                                                       gTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGt                            <  1:74015/45‑1 (MQ=255)                                                                       gTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGt                            <  1:3053862/45‑1 (MQ=255)                                                                       gTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAAt  >  1:4577140/1‑71 (MQ=255)                                                                       |                                                                       GAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAAT  >  minE/1279146‑1279286

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