Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,279,219 A→G 20.2% K333R (AAA→AGA) ‡ yeeJ → adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,279,2190AG20.2% 166.1 / 36.2 84K333R (AAA→AGA) ‡yeeJadhesin
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base G (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 = 9.08e-01

GCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCC  >  Exported/1279149‑1279288
                                                                      |                                                                     
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:4022573/1‑70 (MQ=255)
            aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                            >  1:5107887/1‑70 (MQ=255)
            aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                            >  1:5592350/1‑70 (MQ=255)
            aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                            >  1:1558620/1‑70 (MQ=255)
            aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                            >  1:5648913/1‑70 (MQ=255)
            aaTGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATg                                                            >  1:2544932/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: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: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: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: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)
                                                                      |                                                                     
GCACGCCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCC  >  Exported/1279149‑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: