Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,255,616 Δ1 bp intergenic (‑19/+312) yhaO ← / ← yhaP putative transport system permease protein/putative L‑serine dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,255,6150G.97.8% 399.5 / 1.2 93intergenic (‑18/+313)yhaO/yhaPputative transport system permease protein/putative L‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base . (42/48);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTCAGTACTGTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGGCACGACCCGCCGGGGTCGAAGGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTAGGGCTTA  >  NC_000913/3255539‑3255683
                                                                            |                                                                    
gtgtCAGTACTGTCGAATTTGATCGCTTCTCACCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAg                                                                      <  1:1596466/76‑1 (MQ=255)
 tgtCAGTCCTGTCGAATTTGACCGCTTCTCTCCACTCGCGCTCGCTCATTCCAGCACGCCCCGCCGGGGTCGAAgc                                                                     <  1:485840/76‑2 (MQ=255)
  gtCAGTACTGTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAgcg                                                                    >  1:4940898/1‑74 (MQ=255)
  gtCAGTACTGTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAgcg                                                                    >  1:3726105/1‑74 (MQ=255)
    cattACAGTCGAAGTGGAACACCTCCCGCCGCTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAGGCGTc                                                                 <  1:487869/73‑1 (MQ=255)
     aGTACTGTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCt                                                                >  2:2961639/1‑76 (MQ=255)
     aGTACTGTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCt                                                                >  1:3110281/1‑76 (MQ=255)
     aGTACTGTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCt                                                                >  1:4319161/1‑76 (MQ=255)
      gTACTGTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTg                                                               >  1:1247932/1‑76 (MQ=255)
      gTACTGTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTg                                                               <  1:4986068/76‑1 (MQ=255)
       tACTGTCGATTTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGc                                                              <  1:2199397/76‑1 (MQ=255)
        aCTGTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCa                                                             <  1:5592475/76‑1 (MQ=255)
          tctCGAGTGTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAAt                                                           <  2:355273/74‑1 (MQ=255)
          tGTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAAt                                                           >  2:1081212/1‑76 (MQ=255)
          tGTCGAATTTGACCGCTTCTCGCCTGTCGGTCTCACTCATTCCAGCGCGGCCCGCCTGGGTCGAAGCGTCTGCAAt                                                           <  2:353943/76‑1 (MQ=255)
           gTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATg                                                          >  1:3547386/1‑76 (MQ=255)
           gTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGACGCGTCTGCAATg                                                          >  2:269561/1‑76 (MQ=255)
            tCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGa                                                         <  2:589153/76‑1 (MQ=255)
            tCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGa                                                         >  1:2530645/1‑76 (MQ=255)
             cGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGAc                                                        <  1:63218/76‑1 (MQ=255)
              gAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCGCGACCCGCCGGGGTCGAAGCGTCTGCAATGACg                                                       <  2:5760498/76‑1 (MQ=255)
               actttGATCGCTTCTCGCATCTCGCTCTCACTCCTTCCAGCACGATACGGCGGGTTCTAAGCGTCTGCAATGACGc                                                      <  1:3924460/74‑1 (MQ=255)
               aaTTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGc                                                      <  1:3899594/76‑1 (MQ=255)
                aTTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGcc                                                     >  2:631749/1‑76 (MQ=255)
                aTTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGcc                                                     >  2:2756824/1‑76 (MQ=255)
                   tGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCttt                                                  <  1:2861899/76‑1 (MQ=255)
                   tGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCttt                                                  <  1:4772403/76‑1 (MQ=255)
                    gagcGCTTCGCGCCACTCGCTCTCACTCATTCCAGCACGCCCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTa                                                 <  2:5263021/73‑1 (MQ=255)
                    gATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTa                                                 >  2:3398796/1‑76 (MQ=255)
                    gATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTa                                                 >  2:1066736/1‑76 (MQ=255)
                    gATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTa                                                 >  1:1378449/1‑76 (MQ=255)
                    gATCCCTCCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTa                                                 <  2:1084957/76‑1 (MQ=255)
                     aTCGCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTAt                                                >  2:2170896/1‑76 (MQ=255)
                        gCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCg                                             >  2:4883370/1‑76 (MQ=255)
                        gCTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCg                                             >  2:1098627/1‑76 (MQ=255)
                         cTTCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGa                                            >  1:2776116/1‑76 (MQ=255)
                          ttCTCGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGAt                                           >  1:3979067/1‑76 (MQ=255)
                           tctcGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATg                                          <  2:4331413/76‑1 (MQ=255)
                             tcGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGCTGCa                                        >  1:3746763/1‑76 (MQ=255)
                             tcGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCa                                        <  2:2005729/76‑1 (MQ=255)
                             tcGCCACTCGCTCCCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCa                                        >  2:4249379/1‑76 (MQ=255)
                              cGCCACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCaa                                       <  1:1219680/76‑1 (MQ=255)
                                ccACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAAtt                                     >  1:995031/1‑76 (MQ=255)
                                 cACTCGCTCTCACTCATTCCAGCACGCCCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAAttt                                    <  2:3355933/76‑1 (MQ=255)
                                 cACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAAttt                                    <  1:3381408/76‑1 (MQ=255)
                                 cACTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAAttt                                    >  2:1933776/1‑76 (MQ=255)
                                  ccTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTc                                   <  1:1066736/75‑1 (MQ=255)
                                  ccTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTc                                   <  1:657383/75‑1 (MQ=255)
                                  ccTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTc                                   <  1:828991/75‑1 (MQ=255)
                                   cTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTcc                                  >  1:5617110/1‑76 (MQ=255)
                                   cTCGCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTcc                                  >  1:1293475/1‑76 (MQ=255)
                                      gCTCTCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATa                               >  1:1796056/1‑76 (MQ=255)
                                       ctctCCCTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATaa                              <  2:5404519/76‑1 (MQ=255)
                                       ctctCACTCATTCCAGCGCGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATaa                              <  1:1817326/76‑1 (MQ=255)
                                       ctctCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATaa                              >  1:3018608/1‑76 (MQ=255)
                                       ctctCACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATaa                              <  1:5370791/76‑1 (MQ=255)
                                        tctcACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAAt                             >  1:4150499/1‑76 (MQ=255)
                                        tctcACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAAt                             <  1:3334806/76‑1 (MQ=255)
                                        tctcACTCATTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCAATGCAATTTCCATAAt                             >  2:4919626/1‑76 (MQ=255)
                                         ctcacccaTTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAAtt                            >  1:4580703/1‑76 (MQ=255)
                                          tcactcaTTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCCTAAttt                           >  2:469718/1‑76 (MQ=255)
                                          tcactcaTTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAAttt                           >  2:2112486/1‑76 (MQ=255)
                                          tcactcaTTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAAttt                           <  2:2530645/76‑1 (MQ=255)
                                           cactcaTTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTtc                          >  1:3306062/1‑76 (MQ=255)
                                            actcaTTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTtct                         <  2:3892989/76‑1 (MQ=255)
                                             ctcaTTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTtctc                        <  2:180271/76‑1 (MQ=255)
                                             ctcaTTCCAGCACGACCCGCCGGGGTCGAAGCGTCCGCAATGACGCCTTTATTCGATGCAATTTCCATAATTtctc                        <  1:2803/76‑1 (MQ=255)
                                              tcaTTCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCg                       <  1:643058/76‑1 (MQ=255)
                                               caTTCCAGCACGACCCGCGGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGc                      <  1:2683082/76‑1 (MQ=255)
                                                 ttCCAGCACGACGCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTc                    <  2:2919868/76‑1 (MQ=255)
                                                 ttCCAGCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCCTTATTCGATGCAATTTCCATAATTTCTCCCTc                    >  2:5477760/1‑76 (MQ=255)
                                                    ccagCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCaa                  <  1:5373632/73‑1 (MQ=255)
                                                      agCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAata                >  1:2525371/2‑76 (MQ=255)
                                                      agCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAata                <  2:4315114/75‑1 (MQ=255)
                                                       gCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAatat               <  1:5645730/76‑1 (MQ=255)
                                                       gCACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAatat               >  2:1675352/1‑76 (MQ=255)
                                                        cACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATAtt              <  1:628422/76‑1 (MQ=255)
                                                        cACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATAtt              <  2:3619956/76‑1 (MQ=255)
                                                        cACGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATAtt              <  2:3726105/76‑1 (MQ=255)
                                                         aCGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATaa                              >  1:1796047/1‑59 (MQ=255)
                                                         aCGACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATaa                              <  2:1796047/59‑1 (MQ=255)
                                                           gACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTg           <  1:580357/76‑1 (MQ=255)
                                                           gACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTg           >  2:562236/1‑76 (MQ=255)
                                                           gACCCGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTg           <  2:1795220/76‑1 (MQ=255)
                                                             cccGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTa         <  2:3777147/76‑1 (MQ=255)
                                                             cccGCCGGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTa         >  2:2756463/1‑76 (MQ=255)
                                                              ccgccgGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTAg        <  2:527941/76‑1 (MQ=255)
                                                               cgccgGGGTTGCAGCGTCTGCAATGACGCGTTTATTCGATGCAATTTCCATAATTTTTCGCTCAATATTTTGTAgg       <  2:2573937/76‑1 (MQ=255)
                                                                gccgGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTAggg      >  2:523708/1‑76 (MQ=255)
                                                                 ccgGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTAGGGc     <  1:3549037/76‑1 (MQ=255)
                                                                  cgGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTAGGGCt    >  2:2415694/1‑76 (MQ=255)
                                                                  cgGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTAGGGCt    >  2:1181858/1‑76 (MQ=255)
                                                                  cgGGGTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTAGGGCt    <  1:3104072/76‑1 (MQ=255)
                                                                   ggggTGGAAGGGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTAGGGCtt   <  2:144697/76‑1 (MQ=255)
                                                                   ggggTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTAGGGCtt   >  1:1732614/1‑76 (MQ=255)
                                                                   ggggTCGAAGCGTCTGAAATGACGCCTTTATTCGATGCAATTTCCATAAGTTCTCGCTCTATATTTTGTAGGGcgt   >  1:3661789/1‑74 (MQ=255)
                                                                    gggTCGAAGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTAGGGCtta  <  2:4236071/76‑1 (MQ=255)
                                                                            |                                                                    
GTGTCAGTACTGTCGAATTTGATCGCTTCTCGCCACTCGCTCTCACTCATTCCAGGCACGACCCGCCGGGGTCGAAGGCGTCTGCAATGACGCCTTTATTCGATGCAATTTCCATAATTTCTCGCTCAATATTTTGTAGGGCTTA  >  NC_000913/3255539‑3255683

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