Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,279,224 G→A 19.7% V335I (GTC→ATC) ‡ yeeJ → adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,279,2240GA19.7% 180.6 / 34.8 86V335I (GTC→ATC) ‡yeeJadhesin
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base A (7/10);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

CCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGC  >  Exported/1279154‑1279294
                                                                      |                                                                      
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: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:3825451/71‑1 (MQ=255)
                        cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                                                <  1:173517/71‑1 (MQ=255)
                        cgcAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATg                                                <  1:1738861/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:3053862/45‑1 (MQ=255)
                                                              gTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGt                                    <  1:74015/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)
                                                                     ggTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATg   >  1:5237960/1‑71 (MQ=255)
                                                                      gTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGc  <  1:713950/71‑1 (MQ=255)
                                                                      |                                                                      
CCCGGCCAATGGCTGGGATGTACGCGCAGAAAGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGC  >  Exported/1279154‑1279294

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: