Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,216,746 G→C S252S (TCC→TCG yehZ ← putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,216,7460GC99.3% 363.2 / NA 136S252S (TCC→TCGyehZputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (60/75);  minor base T (0/1);  total (60/76)
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.98e-01

Only 100 of 136 total aligned reads displayed.

CAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACGGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCC  >  NC_000913/2216674‑2216818
                                                                        |                                                                        
cAATGTTTTTGCATCGAGGCGGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGAc                                                                       <  1:301388/76‑1 (MQ=255)
  aTGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTc                                                                     <  1:2497721/76‑1 (MQ=255)
  aTGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTc                                                                     <  1:3262380/76‑1 (MQ=255)
  aTGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGGAGGCACTGTGCCCTTTTCGGGTAATCCCCCAAAACCCACTc                                                                     >  1:4714655/1‑76 (MQ=255)
   tGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCa                                                                    <  1:3295410/76‑1 (MQ=255)
    gTTTTTGCATCGAGGCTGGCGAAGACTGGCTGGAGCCACCGTGCCACTTTCCGGTACTCCCCCCACCCCCACTCAc                                                                   >  2:1968546/1‑76 (MQ=255)
       ttggCATCGGGGCGGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGca                                                                <  2:4878741/73‑1 (MQ=255)
          gCAGCGAGGCGGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATCCTCCCTCAACACCGACTCACGctcca                                                             <  1:4737889/76‑1 (MQ=255)
             tCGAGGCTGGCGAAGACTGGCTGGAGCCCCCGTGCCCTTTTTGGGTACTCCCCCAACACCCACTCACCcaccaccg                                                          >  1:1881865/1‑74 (MQ=255)
               gAGGCTGGCGAAGACTGGCTGTAGGCCCTGTGGCCTTTTCGGGTACTCCCCCCACAACCACCCCCCCACCACTggg                                                        >  2:3597140/1‑75 (MQ=255)
               gAGGCTGGCGAAGACTGGCTGGGGGCCCTGTGGCCTTTTTGGGTAATCCCTCAACAACGCCCCACGCACCACCGGt                                                        >  1:4356928/1‑76 (MQ=255)
                tggCTGGCGAAGACTGTCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTg                                                       <  1:1534013/75‑1 (MQ=255)
                 ggCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGc                                                      <  2:3895370/76‑1 (MQ=255)
                  gCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCa                                                     <  2:2015884/76‑1 (MQ=255)
                   cTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAg                                                    <  1:4407709/76‑1 (MQ=255)
                   cTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAg                                                    <  1:4394229/76‑1 (MQ=255)
                    tGGCGAAGACTGGCTGGAGCCACTGTGCCCTTTGCGGGTACTCCCTCCACACCGGCTCCCGCACCACCGGGGCAgg                                                   >  1:3693126/1‑76 (MQ=255)
                     ggCGACGACTGGCTGTAGCCACTGGGCCATTTGCGGATACTCCCCCCACACCGACCCCCCCCGCCCTGGCGCAGgc                                                  >  2:4342034/1‑76 (MQ=255)
                     ggCGAAGACTGGCTGTAGCCACTGTGCCATTTTCGGATACTCCCCCCACACCCACTCACGCACCACCGGGGCAGgc                                                  >  1:1442439/1‑76 (MQ=255)
                     ggCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGgc                                                  <  1:2622953/76‑1 (MQ=255)
                     ggCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGgc                                                  <  1:5689272/76‑1 (MQ=255)
                     ggCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGgc                                                  <  1:2260347/76‑1 (MQ=255)
                     ggCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGgc                                                  <  1:305826/76‑1 (MQ=255)
                      gggAAGACGGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGgcg                                                 <  2:5440389/74‑1 (MQ=255)
                      gCGAAGACTAGCTGTAGCCACTGTGCCATTTGCGGATACTCCCCCAACCCCGACCCACCCACCCCCGGGGCAAgcc                                                 >  2:755763/1‑75 (MQ=255)
                       cGAAGACTGGGTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGgcgc                                                <  2:2954481/76‑1 (MQ=255)
                       cGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACCCCCCCTCCCGCCCCACCCGCGCAGgcgc                                                >  2:2582669/1‑76 (MQ=255)
                       cGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGgcgc                                                <  1:4408580/76‑1 (MQ=255)
                       cGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGgcgc                                                <  1:4138351/76‑1 (MQ=255)
                        gAAGACTGGCTGTAGCCACTGTGCCATTTTCGGATACTCCCCCAACACCCACCCCCGCACCCCCGGTGCCAGcgcg                                               >  2:1588383/1‑76 (MQ=255)
                        gAAGACTGGCTGTAGCCACTGTGCCATTTGCGGGTACTCCCTCAACCCCCCCCCACGCCCCACCGGTGCAGGcgcg                                               >  1:69594/1‑76 (MQ=255)
                        gAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCCACTCACGCACCACCGGTGCAGGcgcg                                               >  1:2737427/1‑76 (MQ=255)
                        gAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCCACTCACGCACCACCCGCGCAGGcgcg                                               >  1:2737842/1‑76 (MQ=255)
                        gAAGACTGGCTGGAGCCACTGTGCCCTTTGCGGATACTCCCTCAACACCGACCCACCCACCACCGGGGCAGGcgcg                                               >  2:1121658/1‑76 (MQ=255)
                         aaNACTGGCGGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGt                                              <  2:2548525/76‑1 (MQ=255)
                         aaGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGt                                              <  1:5058455/76‑1 (MQ=255)
                           gACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAg                                            >  1:3895610/1‑76 (MQ=255)
                              tGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATa                                         >  1:2297019/1‑76 (MQ=255)
                              tGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGCAGATa                                         >  1:1004201/1‑76 (MQ=255)
                              gggCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATa                                         <  1:668015/75‑1 (MQ=255)
                               ggCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAg                                        >  1:3623708/1‑76 (MQ=255)
                                 cTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGt                                      >  1:4052831/1‑76 (MQ=255)
                                 cTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGt                                      <  1:4722420/76‑1 (MQ=255)
                                 cTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGt                                      >  1:2960829/1‑76 (MQ=255)
                                 cTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGt                                      >  2:858360/1‑76 (MQ=255)
                                  tGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGtt                                     >  2:4628860/1‑76 (MQ=255)
                                   gTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTg                                    <  1:4519824/76‑1 (MQ=255)
                                   gTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTg                                    <  2:3512965/76‑1 (MQ=255)
                                    tAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGc                                   <  1:3489400/76‑1 (MQ=255)
                                    tAGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGc                                   <  2:3794976/76‑1 (MQ=255)
                                     aGCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGca                                  >  2:4499972/1‑76 (MQ=255)
                                      gCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGcac                                 >  1:794003/1‑76 (MQ=255)
                                      gCCACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGcac                                 <  1:508225/76‑1 (MQ=255)
                                       ccACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGcaca                                <  2:5386527/76‑1 (MQ=255)
                                       ccACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGcaca                                >  1:3998675/1‑76 (MQ=255)
                                        cgcTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCacac                               <  1:3204358/74‑1 (MQ=255)
                                        cACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGTCGCGTAGATAGGTTGCacac                               <  1:2083832/76‑1 (MQ=255)
                                        cACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGCTAGGTTGCacac                               >  2:1994883/1‑76 (MQ=255)
                                        cACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCacac                               <  2:5602302/76‑1 (MQ=255)
                                        cACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCacac                               >  1:2780009/1‑76 (MQ=255)
                                        cACTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCacac                               <  1:4130353/76‑1 (MQ=255)
                                         aCTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACAcc                              >  2:5184154/1‑76 (MQ=255)
                                         aCTGTGCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACAcc                              >  2:2740610/1‑76 (MQ=255)
                                           tgtgCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCtt                            <  1:372430/76‑1 (MQ=255)
                                           tgtgCCATTTGCGGATACTCCCTCAACACCGAATCACGCCCCACTGGTGAAGGCGCGTAGATAGCTTGCACACCtt                            >  2:2286938/1‑76 (MQ=255)
                                            gtgCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTg                           <  1:4575707/76‑1 (MQ=255)
                                             tgCCATTGGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGc                          <  1:126423/76‑1 (MQ=255)
                                              gCCATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCg                         <  1:580654/76‑1 (MQ=255)
                                               ccATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCgg                        <  1:1470816/76‑1 (MQ=255)
                                               ccATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCgg                        <  1:3819059/76‑1 (MQ=255)
                                               ccATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCgg                        <  1:1949840/76‑1 (MQ=255)
                                               ccATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTCGCgg                        >  2:5332435/1‑76 (MQ=255)
                                                cATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGa                       <  2:1960065/76‑1 (MQ=255)
                                                cATTTGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGa                       >  1:3480015/1‑76 (MQ=255)
                                                    tGCGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGc                   >  2:4840258/1‑76 (MQ=255)
                                                      cGGATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCtt                 <  1:3675782/76‑1 (MQ=255)
                                                       ggATACTCCCCCAACCCCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTa                <  2:5648351/76‑1 (MQ=255)
                                                        gATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTaa               >  2:2326751/1‑76 (MQ=255)
                                                        gATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTaa               >  1:4358893/1‑76 (MQ=255)
                                                        gATACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTaa               <  2:3975646/76‑1 (MQ=255)
                                                          tACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAgg             <  1:5155279/76‑1 (MQ=255)
                                                          tACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAgg             >  2:1507233/1‑76 (MQ=255)
                                                          tACTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAgg             >  1:4906794/1‑76 (MQ=255)
                                                            cTCCCTCAACGCCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGtt           <  2:1599069/76‑1 (MQ=255)
                                                            cTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGtt           <  2:2296929/76‑1 (MQ=255)
                                                            cTCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGtt           <  2:2538440/76‑1 (MQ=255)
                                                             tCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGttt          <  1:5408250/76‑1 (MQ=255)
                                                             tCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGttt          >  2:3699161/1‑76 (MQ=255)
                                                             tCCCTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGttt          >  2:3114581/1‑76 (MQ=255)
                                                              cccTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTTCACACCTTGCGGATCGCTTAAGGTTTg         <  2:4659645/76‑1 (MQ=255)
                                                              cccTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTg         <  2:4920295/76‑1 (MQ=255)
                                                               ccTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGc        <  1:4532543/76‑1 (MQ=255)
                                                               ccTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGc        <  2:3138673/76‑1 (MQ=255)
                                                               ccTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGc        <  1:1336688/76‑1 (MQ=255)
                                                                cTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCa       <  1:1937861/76‑1 (MQ=255)
                                                                cTCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCa       >  2:4357993/1‑76 (MQ=255)
                                                                 tCAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAg      >  1:388350/1‑76 (MQ=255)
                                                                  cAACACCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGc     >  2:1836686/1‑76 (MQ=255)
                                                                    acacCGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGccc   >  2:3991159/1‑76 (MQ=255)
                                                                     cacTGACTCACGCGCCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGcccc  <  1:3556385/76‑1 (MQ=255)
                                                                        |                                                                        
CAATGTTTTTGCATCGAGGCTGGCGAAGACTGGCTGTAGCCACTGTGCCATTTGCGGATACTCCCTCAACACGGACTCACGCACCACTGGTGCAGGCGCGTAGATAGGTTGCACACCTTGCGGATCGCTTAAGGTTTGCAGCCCC  >  NC_000913/2216674‑2216818

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: