Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,279,208	C→G	23.2%	H329Q (CAC→CAG)	yeeJ →	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,279,208	0	C	G	23.2%	151.1 / 44.8	82	H329Q (CAC→CAG)	yeeJ	adhesin
Reads supporting (aligned to +/- strand):  ref base C (32/31);  new base G (7/12);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

ACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGG  >  minE/1279138‑1279277                                                                       |                                                                      aCGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCAc                                                                       >  1:1178554/1‑71 (MQ=255)  cGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCAcc                                                                      >  1:4281887/1‑71 (MQ=255)   gATTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCt                                                                     <  1:657543/71‑1 (MQ=255)    aTTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCtt                                                                    <  1:2949821/71‑1 (MQ=255)         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:2544932/1‑70 (MQ=255)                        aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                 >  1:5648913/1‑70 (MQ=255)                        aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                 >  1:5592350/1‑70 (MQ=255)                        aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                 >  1:5107887/1‑70 (MQ=255)                        aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                 >  1:1558620/1‑70 (MQ=255)                        aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                 >  1:4022573/1‑70 (MQ=255)                         aTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTgg                                                                  <  1:2004508/52‑1 (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)                           ggCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCg                                                                >  1:3546268/1‑52 (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:4937695/56‑1 (MQ=255)                                   tGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCt                                                    <  1:5445219/56‑1 (MQ=255)                                   tGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATg                                     <  1:3605568/71‑1 (MQ=255)                                   tGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATg                                     >  1:1143794/1‑71 (MQ=255)                                     tACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGa                                                <  1:3076240/58‑1 (MQ=255)                                     tACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGa                                                <  1:749481/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:173517/71‑1 (MQ=255)                                         cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                               <  1:1738861/71‑1 (MQ=255)                                         cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                               <  1:3825451/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:3480601/1‑70 (MQ=255)                                                        aCCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                 >  1:276759/1‑70 (MQ=255)                                                          ccGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtttg                >  1:930567/1‑67 (MQ=255)                                                          ccGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGt                                           <  1:2601327/42‑1 (MQ=255)                                                          ccGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGt                                           <  1:4177020/42‑1 (MQ=255)                                                          ccGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGt                                           <  1:3013033/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:3888036/51‑1 (MQ=255)                                                             ccTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                               <  1:3169322/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:1670750/71‑1 (MQ=255)                                                                    gCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATc    >  1:5469915/1‑71 (MQ=255)                                                                      aCCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCgg  >  1:3774442/1‑71 (MQ=255)                                                                       ccTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCgg  >  1:1661023/1‑70 (MQ=255)                                                                       |                                                                      ACGATTATGAAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGG  >  minE/1279138‑1279277

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