Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,279,264 A→G 11.0% D348G (GAT→GGT)  yeeJ → adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,279,2640AG11.0% 243.8 / 16.3 91D348G (GAT→GGT) yeeJadhesin
Reads supporting (aligned to +/- strand):  ref base A (50/31);  new base G (7/3);  total (57/34)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCCTTCCCGCTGA  >  Exported/1279196‑1279332
                                                                    |                                                                    
cGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGAt                                                                     >  1:399332/1‑70 (MQ=255)
 gCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaa                                                                   <  1:4646674/71‑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)
         gCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATc                                                           >  1:5469915/1‑71 (MQ=255)
         gCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATc                                                           <  1:1670750/71‑1 (MQ=255)
           aCCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCgg                                                         >  1:3774442/1‑71 (MQ=255)
            ccTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCgg                                                         >  1:1661023/1‑70 (MQ=255)
             cTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATa                                                                    <  1:4725793/58‑1 (MQ=255)
             cTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATa                                                                    <  1:5462014/58‑1 (MQ=255)
                 gCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAg                                                    <  1:4415364/70‑1 (MQ=255)
                 gCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAg                                                    <  1:5163512/70‑1 (MQ=255)
                  cGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTa                                                  <  1:4158708/71‑1 (MQ=255)
                    gTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAAt                                                >  1:4577140/1‑71 (MQ=255)
                      aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                                     >  1:3612501/1‑64 (MQ=255)
                      aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATcc                                              >  1:3678873/1‑71 (MQ=255)
                      aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATcc                                              >  1:4000099/1‑71 (MQ=255)
                           ggTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATg                                         >  1:5237960/1‑71 (MQ=255)
                            gTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGc                                        <  1:713950/71‑1 (MQ=255)
                             tCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCa                                                       <  1:5048629/55‑1 (MQ=255)
                             tCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCa                                                       <  1:4801790/55‑1 (MQ=255)
                              cTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGc                                        >  1:2481538/1‑69 (MQ=255)
                               tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa                                                             <  1:1009565/47‑1 (MQ=255)
                               tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa                                                             <  1:4473728/47‑1 (MQ=255)
                               tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTa                                                  <  1:1886378/58‑1 (MQ=255)
                                 tGAACAGTATTATGGCTATGAAGTCGGGCTGTTTGGTaaa                                                                  >  1:947842/1‑40 (MQ=16)
                                 tGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGc                                                        >  1:132388/1‑50 (MQ=255)
                                 tGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATaa                                   >  1:1902824/1‑71 (MQ=255)
                                 tGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATa                                    <  1:4411527/70‑1 (MQ=255)
                                 tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                                  >  1:4522852/1‑40 (MQ=16)
                                 tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                                  >  1:644916/1‑40 (MQ=16)
                                 tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                                  >  1:5562060/1‑40 (MQ=16)
                                 tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                                  >  1:247319/1‑40 (MQ=16)
                                 tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGa                                                                >  1:547210/1‑42 (MQ=19)
                                 tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGa                                                                >  1:3504138/1‑42 (MQ=19)
                                  gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACgct                                                            >  1:1865255/1‑43 (MQ=255)
                                  gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                                     >  1:5495788/1‑52 (MQ=255)
                                  gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                                     >  1:4982896/1‑52 (MQ=255)
                                  gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                                     >  1:1861234/1‑52 (MQ=255)
                                   aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa                                                             >  1:1065708/1‑43 (MQ=255)
                                   aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGAt                                                            >  1:1696726/1‑44 (MQ=255)
                                   aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGAt                                                            >  1:390284/1‑44 (MQ=255)
                                   aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAAcc                                 >  1:3649832/1‑71 (MQ=255)
                                    aCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCg                                <  1:2886515/71‑1 (MQ=255)
                                       gTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCa                                           <  1:3520479/57‑1 (MQ=255)
                                         attatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCg                                                          >  1:1046219/1‑40 (MQ=255)
                                         attatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCg                                                          >  1:3319100/1‑40 (MQ=255)
                                          ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCa                                           >  1:1897872/1‑54 (MQ=255)
                                          ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACGGCTgg                            <  1:2634890/69‑1 (MQ=255)
                                          ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTgg                            <  1:3714117/69‑1 (MQ=255)
                                          ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTgg                            <  1:690727/69‑1 (MQ=255)
                                          ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGAc                          <  1:409899/71‑1 (MQ=255)
                                          ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGAc                          >  1:3419483/1‑71 (MQ=255)
                                           tatGGCCATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATaa                                   >  1:4468066/1‑61 (MQ=255)
                                               gCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAAc                     <  1:795407/71‑1 (MQ=255)
                                                cGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAAc                     >  1:3659153/1‑70 (MQ=255)
                                                 gATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACta                   >  1:260368/1‑71 (MQ=255)
                                                  aTGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa                                                 <  1:4151459/40‑1 (MQ=255)
                                                   tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACtata                 <  1:1159844/71‑1 (MQ=255)
                                                    gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCa                                      >  1:5288081/1‑49 (MQ=255)
                                                    gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCa                                      >  1:773287/1‑49 (MQ=255)
                                                    gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCa                                      >  1:1669836/1‑49 (MQ=255)
                                                    gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATAc                >  1:4278948/1‑71 (MQ=255)
                                                       gTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATAcccc             >  1:3522912/1‑71 (MQ=255)
                                                        tGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTc                                            >  1:1573677/1‑39 (MQ=255)
                                                        tGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTg                             <  1:1941665/54‑1 (MQ=255)
                                                        tGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATAccccc            >  1:5304496/1‑71 (MQ=255)
                                                           cccTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACta                   >  1:3542256/1‑61 (MQ=255)
                                                             cTGTTCGATAAAGACTATCGGCAAAGTAATCCTCATg                                         <  1:2206558/37‑1 (MQ=255)
                                                             cTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCCTTccc       <  1:3262022/71‑1 (MQ=255)
                                                             cTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCCTTccc       >  1:5117757/1‑71 (MQ=255)
                                                              tGTTCGATAAAGACGATCGGCAAAGTAATCCTCATg                                         <  1:3521756/36‑1 (MQ=255)
                                                              tGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCCTTCCCg      >  1:3563214/1‑71 (MQ=255)
                                                               gTTCGATAAAGACGATCGGCAAAGTAATCCTCATGcc                                       >  1:3492291/1‑37 (MQ=255)
                                                               gTTCGATAAAGACGATCGGCAAAGTAATCCTCATGcc                                       >  1:5581572/1‑37 (MQ=255)
                                                               gTTCGATAAAGACGATCGGCAAAGTAATCCTCATGcc                                       >  1:2481499/1‑37 (MQ=255)
                                                               gTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCAt                                     <  1:4838388/39‑1 (MQ=255)
                                                               gTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAAc                     <  1:5158740/55‑1 (MQ=255)
                                                                ttCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATa                                    >  1:2539876/1‑39 (MQ=255)
                                                                ttCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATa                                    >  1:1324100/1‑39 (MQ=255)
                                                                ttCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAAc                                  >  1:2315906/1‑41 (MQ=255)
                                                                 tCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCCTTCCCGCt    >  1:3023953/1‑70 (MQ=255)
                                                                  cGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCt                              >  1:1756323/1‑43 (MQ=255)
                                                                  cGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTaa                      >  1:5602586/1‑51 (MQ=255)
                                                                  cGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTa                       <  1:4230057/50‑1 (MQ=255)
                                                                  cGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATAccccc            >  1:1101913/1‑61 (MQ=255)
                                                                  cGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCCTTCCCGCtga  >  1:3359988/1‑71 (MQ=255)
                                                                  cGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCCTTCCCGCtga  >  1:3366479/1‑71 (MQ=255)
                                                                    aTAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACta                   <  1:1911635/52‑1 (MQ=255)
                                                                    aTAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACta                   <  1:5194146/52‑1 (MQ=255)
                                                                    |                                                                    
CGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCCTTCCCGCTGA  >  Exported/1279196‑1279332

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: