Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,215,798 C→T 6.9% A227T (GCA→ACA)  yehW ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,215,7980CT 6.9% 181.3 / 2.7 58A227T (GCA→ACA) yehWputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base T (2/2);  total (26/32)
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.60e-01

ATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGAT  >  NC_000913/2215657‑2215936
                                                                                                                                             |                                                                                                                                          
aTGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATc                                                                                                                                     >  1:148776/1‑149 (MQ=255)
          cTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCa                                                                                                                              >  2:235139/1‑146 (MQ=255)
          cTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCa                                                                                                                              <  1:235139/146‑1 (MQ=255)
             aaaTTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGc                                                                                                                          <  6:177844/147‑1 (MQ=255)
                ttGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGccacca                                                                                                                     <  7:61708/149‑1 (MQ=255)
                ttGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGccac                                                                                                                       >  7:2392/1‑147 (MQ=255)
                  ggggTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAg                                                                                                                    <  4:185521/148‑1 (MQ=255)
                   gggTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTa                                                                                                                  <  3:207379/149‑1 (MQ=255)
                     gTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAAc                                                                                                                <  5:43425/149‑1 (MQ=255)
                     gTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAAc                                                                                                                >  1:87574/1‑149 (MQ=255)
                        aTCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGcc                                                                                                             >  3:22368/1‑149 (MQ=255)
                            cccTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCt                                                                                                         <  4:223984/149‑1 (MQ=255)
                                 ttcctgttTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGAt                                                                                                     >  7:167768/1‑148 (MQ=255)
                                    ctgttTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATcac                                                                                                  <  6:155134/148‑1 (MQ=255)
                                     tgttTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGc                                                                                                              >  3:155053/1‑135 (MQ=255)
                                     tgttTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGc                                                                                                              <  4:155053/135‑1 (MQ=255)
                                     tgttTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATcaca                                                                                                 >  1:208535/1‑148 (MQ=255)
                                                ttAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTa                                                                                     <  1:237819/149‑1 (MQ=255)
                                                   ggCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAAt                                                                                  >  3:111441/1‑149 (MQ=255)
                                                    gCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATc                                                                                 <  1:147578/149‑1 (MQ=255)
                                                    gCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATc                                                                                 <  7:111804/149‑1 (MQ=255)
                                                    gCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATc                                                                                 >  1:93492/1‑149 (MQ=255)
                                                          ttGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTaa                                                                           <  6:10837/149‑1 (MQ=255)
                                                                gCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCga                                                                     <  4:22368/149‑1 (MQ=255)
                                                                     aTGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCt                                                                                                         >  2:121327/1‑108 (MQ=255)
                                                                     aTGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCt                                                                                                         <  1:121327/108‑1 (MQ=255)
                                                                        cTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGatgatgatg                                                             <  8:167768/149‑1 (MQ=255)
                                                                          cGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATggg                                                           <  3:148966/149‑1 (MQ=255)
                                                                           gCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATgg                                                            <  2:87574/147‑1 (MQ=255)
                                                                           gCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATgg                                                            <  2:208535/147‑1 (MQ=255)
                                                                           gCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATgg                                                            <  2:148776/147‑1 (MQ=255)
                                                                               gTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTAc                                                      <  3:228834/149‑1 (MQ=255)
                                                                                      cTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGt                                               >  3:230511/1‑149 (MQ=255)
                                                                                       tttATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTg                                              >  4:37630/1‑149 (MQ=255)
                                                                                           ttttGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTgg                                          >  1:24421/1‑149 (MQ=255)
                                                                                           ttttGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTg                                           <  1:16883/148‑1 (MQ=255)
                                                                                           ttttGCGTGCT‑GCTAAGCGCCTGCACCAGCCTTTCAACCAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGc                                         >  5:66049/1‑149 (MQ=255)
                                                                                                      ggCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCa                                                   >  1:98686/1‑129 (MQ=255)
                                                                                                      ggCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCa                                                   <  2:98686/129‑1 (MQ=255)
                                                                                                      ggCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGtt                               >  3:124612/1‑149 (MQ=255)
                                                                                                       gCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGtt                               >  3:34060/1‑148 (MQ=255)
                                                                                                           aGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAg                            <  2:93492/147‑1 (MQ=255)
                                                                                                                   cACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtcgtc                   <  8:2392/148‑1 (MQ=255)
                                                                                                                   cACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtcgtc                   <  3:37630/148‑1 (MQ=255)
                                                                                                                   cACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtcgtc                   <  4:34060/148‑1 (MQ=255)
                                                                                                                    aCCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGc                 >  2:205062/1‑149 (MQ=255)
                                                                                                                      cAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCg               >  5:184451/1‑149 (MQ=255)
                                                                                                                      cAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCg               >  7:63140/1‑149 (MQ=255)
                                                                                                                      cAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCg               <  4:111441/149‑1 (MQ=255)
                                                                                                                           tttCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTAcc           >  2:201962/1‑148 (MQ=255)
                                                                                                                            ttCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCaa         >  1:67992/1‑149 (MQ=255)
                                                                                                                               aaaCAGGCGGTCTGTGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtc                      <  5:242358/133‑1 (MQ=255)
                                                                                                                               aaaCAGGCGGTCTGTGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtc                      >  6:242358/1‑133 (MQ=255)
                                                                                                                               aaaCAGGCGGTCTGTGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtc                      <  7:122586/133‑1 (MQ=255)
                                                                                                                               aaaCAGGCGGTCTGTGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGAtc                      >  8:122586/1‑133 (MQ=255)
                                                                                                                                aaCAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGtt     >  3:81375/1‑149 (MQ=255)
                                                                                                                                 aCAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTg    <  2:24421/149‑1 (MQ=255)
                                                                                                                                   aGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGAt  <  3:185968/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                          
ATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGAT  >  NC_000913/2215657‑2215936

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑