Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,279,220	A→C	20.3%	K333N (AAA→AAC) ‡	yeeJ →	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,279,220	0	A	C	20.3%	135.8 / 45.2	84	K333N (AAA→AAC) ‡	yeeJ	adhesin
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base C (7/10);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 7.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

ACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCC  >  Exported/1279151‑1279288                                                                      |                                                                     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:4022573/1‑70 (MQ=255)           aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                            >  1:5107887/1‑70 (MQ=255)           aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                            >  1:5592350/1‑70 (MQ=255)           aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                            >  1:5648913/1‑70 (MQ=255)           aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                            >  1:2544932/1‑70 (MQ=255)           aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                            >  1:1558620/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: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:1738861/71‑1 (MQ=255)                            cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                                          <  1:3825451/71‑1 (MQ=255)                            cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                                          <  1:173517/71‑1 (MQ=255)                            cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                                          <  1:908132/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:2601327/42‑1 (MQ=255)                                             ccGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGt                                                      <  1:3013033/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:1670750/71‑1 (MQ=255)                                                       gCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATc               >  1:5469915/1‑71 (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)                                                                   tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTg                               <  1:4469042/43‑1 (MQ=255)                                                                    aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                            >  1:4641011/1‑45 (MQ=255)                                                                    aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa         >  1:3612501/1‑64 (MQ=255)                                                                    aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATcc  >  1:3678873/1‑71 (MQ=255)                                                                    aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATcc  >  1:4000099/1‑71 (MQ=255)                                                                     aaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGtt                             >  1:1547429/1‑43 (MQ=255)                                                                      |                                                                     ACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCC  >  Exported/1279151‑1279288

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