Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,218,525 C→A A11A (GCG→GCT yehY ← glycine betaine ABC transporter membrane subunit YehY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,218,5250CA98.1% 180.2 / ‑2.8 53A11A (GCG→GCTyehYglycine betaine ABC transporter membrane subunit YehY
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (29/23);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCAA  >  NC_000913/2218388‑2218656
                                                                                                                                         |                                                                                                                                   
caaccaacGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAgcgc                                                                                                                                   <  1:1930057/140‑1 (MQ=255)
caaccaacGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAgagc                                                                                                                                   <  2:1486134/140‑1 (MQ=255)
caaccaacGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAgagc                                                                                                                                   <  2:1497868/140‑1 (MQ=255)
 aaccaacGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGcc                                                                                                                                  <  2:1989868/140‑1 (MQ=255)
     aacGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGaa                                                                                                                              <  2:2148751/140‑1 (MQ=255)
     aacGCCCACCAGCACCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGaa                                                                                                                              <  2:1277448/140‑1 (MQ=255)
      acGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAAc                                                                                                                             <  2:1690714/140‑1 (MQ=255)
       cGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACa                                                                                                                            >  1:468130/1‑140 (MQ=255)
       cGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACa                                                                                                                            >  1:357441/1‑140 (MQ=255)
           caccaGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCCCCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGAt                                                                                                                        <  2:303332/140‑1 (MQ=255)
           caccaGCATCCAGAACGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACGAAACGATTAGGCGCGTAACTCATAAACGGCAGCGCCGCTGCGATTGCCGTCAGAAACAGCAGCAGAGCCAGAAAAGGAt                                                                                                                        <  2:710150/140‑1 (MQ=255)
             ccaGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTa                                                                                                                      >  1:894956/1‑140 (MQ=255)
             ccaGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTa                                                                                                                      >  1:1929558/1‑140 (MQ=255)
                gCATCCGGATCGTTTGCGGCCACAGCTGCGTGAGATGGCGCCCCTCACCCGAAGCTAAACGATTAGGCGCGTAACTGACAATCGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCTATAACGGGATTAATa                                                                                                                   <  2:375196/140‑1 (MQ=255)
                gCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGATACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATa                                                                                                                   <  1:1091814/140‑1 (MQ=255)
                gCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATa                                                                                                                   <  1:648579/140‑1 (MQ=255)
                        aTCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATa                                                                                                           >  2:1141991/1‑140 (MQ=255)
                           gTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATAc                                                                                                                  <  1:351580/130‑1 (MQ=255)
                           gTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATAc                                                                                                                  >  2:351580/1‑130 (MQ=255)
                           gTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGt                                                                                                        <  2:360793/140‑1 (MQ=255)
                            tttGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTc                                                                                                       >  1:2037617/1‑140 (MQ=255)
                            tttGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTc                                                                                                       >  1:600654/1‑140 (MQ=255)
                             ttGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCa                                                                                                      <  2:1929558/140‑1 (MQ=255)
                               gCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACg                                                                                                    >  1:462057/1‑140 (MQ=255)
                               gCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACg                                                                                                    <  1:503924/140‑1 (MQ=255)
                                  gCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGa                                                                                                 >  2:783200/1‑140 (MQ=255)
                                   ccACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGaa                                                                                                <  1:1233211/140‑1 (MQ=255)
                                    cacaGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGaaa                                                                                               >  2:810999/1‑140 (MQ=255)
                                    cacaGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGaaa                                                                                               >  2:697710/1‑140 (MQ=255)
                                      caGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAAtt                                                                                             >  1:972517/1‑140 (MQ=255)
                                      caGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAAtt                                                                                             >  1:1561479/1‑140 (MQ=255)
                                         cTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACt                                                                                          <  1:810999/140‑1 (MQ=255)
                                              agagATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTc                                                                                     <  2:514227/140‑1 (MQ=255)
                                               gagaTGACGCCCCTCCCCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTcc                                                                                    <  2:2066024/140‑1 (MQ=255)
                                                agaTGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCa                                                                                   <  1:652477/140‑1 (MQ=255)
                                                   tGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCAccc                                                                                >  1:2056019/1‑140 (MQ=255)
                                                   tGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCAccc                                                                                >  1:939456/1‑140 (MQ=255)
                                                       gCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCtttt                                                                            >  2:2060932/1‑140 (MQ=255)
                                                        ccccTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCttttt                                                                           <  2:939456/140‑1 (MQ=255)
                                                          ccTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGt                                                                         <  2:2037617/140‑1 (MQ=255)
                                                           cTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTc                                                                                     >  2:316321/1‑127 (MQ=255)
                                                           cTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTc                                                                                     <  1:316321/127‑1 (MQ=255)
                                                                   aaaCTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAg                                                                >  1:816931/1‑140 (MQ=255)
                                                                       tAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCTCCCTTTTTGTTTCAGGTAGTCgg                                                            <  2:1911590/140‑1 (MQ=255)
                                                                        aaaCGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGc                                                           <  2:1561479/140‑1 (MQ=255)
                                                                         aaCGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCa                                                          >  1:347245/1‑140 (MQ=255)
                                                                          aCGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAg                                                         >  1:2028107/1‑140 (MQ=255)
                                                                                   cgcgTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACttttt                                                >  1:479329/1‑140 (MQ=255)
                                                                                            gATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAAATACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCa                                       >  2:1817771/1‑140 (MQ=255)
                                                                                                 aCGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCt                                  >  1:686258/1‑140 (MQ=255)
                                                                                                    gCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTcc                               >  2:975334/1‑140 (MQ=255)
                                                                                                    gCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTcc                               >  1:1648253/1‑140 (MQ=255)
                                                                                                    gCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTcc                               >  2:2089345/1‑140 (MQ=255)
                                                                                                     cAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCca                              <  2:468130/140‑1 (MQ=255)
                                                                                                     cAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCca                              >  1:1751635/1‑140 (MQ=255)
                                                                                                     cAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCca                              >  1:1752273/1‑140 (MQ=255)
                                                                                                     cAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGAGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCca                              >  1:512214/1‑140 (MQ=255)
                                                                                                                                 agcagcagaGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCaa  <  2:1751635/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                   
CAACCAACGCCCACCAGCATCCAGATCGTTTGCGGCCACAGCTGCCAGAGATGACGCCCCTCACCCGAAACTAAACGATTAGGCGCGTAACTGATAAACGGCAGCGCCGCTGCGATTGCCGTCAGCAACAGCAGCAGCGCCAGAACAGGATTAATACGGAAATAAGTCACGGGAAATTACTTCGTCCACCCTTTTTGTTTCAGGTAGTCGGCAGCCACTTTTTTGGCATCCAGTCCTTCCACAGCAATGCTGGCATTCAGTTGCTGCAA  >  NC_000913/2218388‑2218656

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: