Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,279,217	T→A	20.0%	G332G (GGT→GGA) ‡	yeeJ →	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,279,217	0	T	A	20.0%	191.9 / 41.8	85	G332G (GGT→GGA) ‡	yeeJ	adhesin
Reads supporting (aligned to +/- strand):  ref base T (32/36);  new base A (7/10);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

AAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAAT  >  Exported/1279147‑1279286                                                                       |                                                                      aaGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGt                                                                       <  1:4941040/71‑1 (MQ=255) aaGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGt                                                                       <  1:2814543/71‑1 (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:5592350/1‑70 (MQ=255)               aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                          >  1:5648913/1‑70 (MQ=255)               aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                          >  1:5107887/1‑70 (MQ=255)               aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                          >  1:2544932/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:4937695/56‑1 (MQ=255)                          tGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCt                                                             <  1:5445219/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:4741682/71‑1 (MQ=255)                                 gcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGa                                       <  1:1164305/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:4802785/54‑1 (MQ=255)                                             ttACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCg                                           <  1:3885847/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:1411955/53‑1 (MQ=255)                                                  cGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATg                                        <  1:5291644/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)                                                              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)                                                                       |                                                                      AAGCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAAT  >  Exported/1279147‑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: ‑