Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 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
*Exported1,279,2080CG23.2% 151.1 / 44.8 82H329Q (CAC→CAGyeeJadhesin
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  >  Exported/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  >  Exported/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: