Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,870,631 Δ1 bp coding (610/1764 nt) dgoA ← 2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,870,6300C.99.4% 797.8 / ‑2.7 178coding (611/1764 nt)dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base . (84/91);  minor base G (0/1);  total (84/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 186 total aligned reads displayed.

TCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGA  >  NC_000913/3870555‑3870701
                                                                           |                                                                       
tCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATCTTGGTGATTTTCATTGCAc                                                                         <  1:3495595/76‑1 (MQ=255)
   tCGGTTTCGATTTTCGCGAGCGTCCAGCGGGGCGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActgc                                                                      <  2:5664021/76‑4 (MQ=255)
     ggTTTCAATTTTCAGGAACATCCCTCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCct                                                                   <  1:5374958/76‑1 (MQ=255)
     ggTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCct                                                                   >  2:5016921/1‑76 (MQ=255)
      gTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCctc                                                                  >  2:2485165/1‑76 (MQ=255)
       tttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATAAGTGGTAATTTTGGTGATTTTCATTGCACTGCctct                                                                 >  1:508290/1‑76 (MQ=255)
        ttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCtctc                                                                >  2:433461/1‑76 (MQ=255)
          cAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                              <  1:4103940/76‑1 (MQ=255)
             ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATat                                                           <  1:1352496/76‑2 (MQ=255)
             ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATAc                                                           <  2:4833524/76‑1 (MQ=255)
             ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATAc                                                           <  1:3952270/76‑1 (MQ=255)
              tttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACg                                                          <  1:2960114/76‑1 (MQ=255)
              tttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACg                                                          <  2:1211031/76‑1 (MQ=255)
                tCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGcc                                                        <  2:5322792/76‑1 (MQ=255)
                 cAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCt                                                       <  1:1630929/76‑1 (MQ=255)
                  aGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCtt                                                      <  1:3678469/76‑1 (MQ=255)
                  aGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCtt                                                      >  1:3491738/1‑76 (MQ=255)
                   ggAACATCCAGCGGGGAGGTGAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGTCTTa                                                     >  1:237823/1‑76 (MQ=255)
                   ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTa                                                     >  1:4187948/1‑76 (MQ=255)
                   ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTa                                                     >  1:4126423/1‑76 (MQ=255)
                     aaCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAAc                                                   <  1:2337407/76‑1 (MQ=255)
                      aCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACa                                                  <  2:4760736/76‑1 (MQ=255)
                      aCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACa                                                  >  1:3045690/1‑76 (MQ=255)
                       cATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaa                                                 >  2:2133602/1‑76 (MQ=255)
                       cATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaa                                                 >  1:1718144/1‑76 (MQ=255)
                         tCCAGCGGGGAGGTACACGATACGTGTTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAt                                               <  1:1147555/76‑1 (MQ=255)
                          ccAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtg                                              <  1:411773/76‑1 (MQ=255)
                          ccAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtg                                              <  1:2781763/76‑1 (MQ=255)
                           cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgc                                             <  1:602142/76‑1 (MQ=255)
                           cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgc                                             >  2:2422126/1‑76 (MQ=255)
                           cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgc                                             >  2:890939/1‑76 (MQ=255)
                            aGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgct                                            >  1:213347/1‑76 (MQ=255)
                            aGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgct                                            >  1:4717535/1‑76 (MQ=255)
                              cGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgc                                          >  1:5671404/1‑76 (MQ=255)
                              cGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgc                                          >  1:1716097/1‑76 (MQ=255)
                               ggggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgct                                         >  1:1428300/1‑76 (MQ=255)
                               ggggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgct                                         <  1:2995185/76‑1 (MQ=255)
                                gggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctg                                        >  2:4554767/1‑76 (MQ=255)
                                gggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctg                                        >  2:1904710/1‑76 (MQ=255)
                                gggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAt                                               >  1:5220977/1‑69 (MQ=255)
                                gggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAt                                               <  2:5220976/69‑1 (MQ=255)
                                 ggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                       <  2:3734252/76‑1 (MQ=255)
                                 ggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                       >  1:711032/1‑76 (MQ=255)
                                  gAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcc                                      >  1:2261922/1‑76 (MQ=255)
                                  gAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcc                                      <  2:1706538/76‑1 (MQ=255)
                                  gAGGTAAACGATACGTGGTAATTTTGGTGATATTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcc                                      <  1:5048595/76‑1 (MQ=255)
                                   aGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcct                                     >  1:191652/1‑76 (MQ=255)
                                   aGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcct                                     <  1:5587761/76‑1 (MQ=255)
                                   aGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcct                                     >  2:2659192/1‑76 (MQ=255)
                                     gTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgc                                   <  1:4644889/76‑1 (MQ=255)
                                      tAAACGATGCGTGGTATTTTTTGTGATTTTCATGGCACTGCCTCTCGATACGCCTTACCAAATGctgctgcctgct                                  <  1:1866495/76‑1 (MQ=255)
                                      tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                  <  1:4905128/76‑1 (MQ=255)
                                      tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                  >  1:2059763/1‑76 (MQ=255)
                                      tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                  >  1:4333813/1‑76 (MQ=255)
                                       aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctg                                 >  1:4296840/1‑76 (MQ=255)
                                       aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGccgctgaatgctg                                 >  1:2707976/1‑76 (MQ=255)
                                        aaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctgc                                >  2:362591/1‑76 (MQ=255)
                                        aaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctgc                                <  2:2860162/76‑1 (MQ=255)
                                        aaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctgc                                <  2:646049/76‑1 (MQ=255)
                                        aaCGATACGTGGAAATTTTGGGGATTTTCATTGAACTGCCTCTCGATACGCATTAACAAATGcttctgcctcctgc                                >  2:3791088/1‑76 (MQ=255)
                                         aCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcg                               <  2:4506228/76‑1 (MQ=255)
                                          cGATACGTGGTATTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcgc                              <  2:1876449/76‑1 (MQ=255)
                                          cGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcgc                              >  2:5184762/1‑76 (MQ=255)
                                           gATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGcgcg                             <  1:4980461/76‑1 (MQ=255)
                                            aTACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCgg                            >  2:5040693/1‑76 (MQ=255)
                                            aTACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCgg                            >  1:1520208/1‑76 (MQ=255)
                                              aCGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTg                          >  2:1548007/1‑76 (MQ=255)
                                               cGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAANNCTGCTGCNTGCTGCGCGGTnc                         <  1:2684214/76‑1 (MQ=255)
                                                gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcg                        >  2:40754/1‑76 (MQ=255)
                                                gTGGCATTTTTGGTGATTTGCATTGCCCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcg                        <  1:4396065/76‑1 (MQ=255)
                                                 tGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcgc                       >  1:2928577/1‑76 (MQ=255)
                                                 aggTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcgc                       >  2:2045852/2‑76 (MQ=255)
                                                    tAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGctct                    <  1:1146075/76‑1 (MQ=255)
                                                    tAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGctct                    <  1:4787894/76‑1 (MQ=255)
                                                    tAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGctct                    <  2:2237898/76‑1 (MQ=255)
                                                    tAATTTTGGGGGTTTTCGTTGCGCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGctct                    <  2:5170422/76‑1 (MQ=255)
                                                      aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                  >  2:4890234/1‑76 (MQ=255)
                                                      aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                  <  1:1976529/76‑1 (MQ=255)
                                                       ttttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACg                 <  2:2851744/76‑1 (MQ=255)
                                                        tttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACgg                <  2:1596736/76‑1 (MQ=255)
                                                         ttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGa               >  2:3737933/1‑76 (MQ=255)
                                                         ttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGa               <  1:5367795/76‑1 (MQ=255)
                                                          tGGTGTTTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAt              <  1:260586/76‑1 (MQ=255)
                                                           ggTGCTTTTCATTGCCCTGCCTCACGATACGCCTTAACACCTGCTGCTGCCTGCTCCGCGGTGCGCTCTCCGGAtt             <  1:839680/76‑1 (MQ=255)
                                                           ggTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAtt             <  2:3923434/76‑1 (MQ=255)
                                                            gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTg            >  2:4207296/1‑76 (MQ=255)
                                                             tGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGc           >  1:3399241/1‑76 (MQ=255)
                                                               aTTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGccc         >  2:3564262/1‑76 (MQ=255)
                                                               aTTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGccc         <  2:3491738/76‑1 (MQ=255)
                                                                ttttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCg        >  2:1165678/1‑76 (MQ=255)
                                                                 tttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGCGCGCTCTACGGCTCGccccg       >  2:5705646/1‑74 (MQ=255)
                                                                  ttCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctg                                 <  2:4680491/49‑1 (MQ=38)
                                                                  ttCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctg                                 >  1:4680492/1‑49 (MQ=38)
                                                                  ttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgc      >  2:2681026/1‑76 (MQ=255)
                                                                   tCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcg     <  1:1539199/76‑1 (MQ=255)
                                                                    cATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc    <  1:2276168/76‑1 (MQ=255)
                                                                    cATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc    <  1:3563143/76‑1 (MQ=255)
                                                                     aTTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGcgcg   >  1:3960296/1‑76 (MQ=255)
                                                                     aTTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGcgcg   <  2:2928577/76‑1 (MQ=255)
                                                                      ttGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGa  <  1:2534710/76‑1 (MQ=255)
                                                                           |                                                                       
TCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGA  >  NC_000913/3870555‑3870701

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