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.98.5% 612.9 / ‑2.0 138coding (611/1764 nt)dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (71/64);  total (72/66)
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 147 total aligned reads displayed.

TCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGA  >  NC_000913/3870555‑3870701
                                                                           |                                                                       
tCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                         <  2:2378662/76‑1 (MQ=255)
  aTCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActg                                                                       <  2:2424443/76‑3 (MQ=255)
   tCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActgc                                                                      <  2:2640679/76‑4 (MQ=255)
   tCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActgc                                                                      >  2:2694412/1‑73 (MQ=255)
     ggTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCct                                                                   <  2:1463276/76‑1 (MQ=255)
     ggTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGGAATTTTGGTGATCATCATTGCACTGCct                                                                   >  1:3202371/1‑76 (MQ=255)
       tttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCctct                                                                 >  2:429423/1‑76 (MQ=255)
       tttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCctct                                                                 >  2:2239049/1‑76 (MQ=255)
         tCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCg                                                               >  1:2441282/1‑76 (MQ=255)
          cAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                              >  2:409197/1‑76 (MQ=255)
           aaTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGAt                                                             <  1:2329246/76‑1 (MQ=255)
            aTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATa                                                            >  1:1003624/1‑76 (MQ=255)
             ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATAc                                                           <  2:4222663/76‑1 (MQ=255)
             ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATAc                                                           <  2:3456083/76‑1 (MQ=255)
              tttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGacacc                                                          >  2:3654361/1‑75 (MQ=255)
              tttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACg                                                          <  1:1064884/76‑1 (MQ=255)
               tgcAGGAACATCCAGAGGGGAGGCAAACGATACGGGGGAAGTTTGGTGATTTTCATTGCACTGCCTATCGATACGc                                                         >  1:1773009/3‑76 (MQ=255)
                tCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGcc                                                        >  1:3330903/1‑76 (MQ=255)
                tCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGcc                                                        >  2:3062760/1‑76 (MQ=255)
                tCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGcc                                                        >  2:2346719/1‑76 (MQ=255)
                 cAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCt                                                       >  2:4348990/1‑76 (MQ=255)
                  aGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCtt                                                      <  1:4705508/76‑1 (MQ=255)
                   ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTa                                                     <  1:174967/76‑1 (MQ=255)
                    gAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTaa                                                    <  1:3619726/76‑1 (MQ=255)
                     aaCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAAc                                                   <  2:4135584/76‑1 (MQ=255)
                     aaCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAAc                                                   >  2:1890086/1‑76 (MQ=255)
                     aaCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAAc                                                   >  1:984375/1‑76 (MQ=255)
                      aCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACa                                                  >  1:2533266/1‑76 (MQ=255)
                      aCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACa                                                  >  1:4048256/1‑76 (MQ=255)
                      aCATCCAGCGGGGAGGTAAACGAGACGTGGTAATTTTGGTGAGTGTCAATGCACTGCCTCTCGAGACGCCTTAACa                                                  >  1:1763364/1‑76 (MQ=255)
                        aTCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTaaaaaa                                                <  1:409197/76‑1 (MQ=255)
                        aTCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaaa                                                >  1:4251298/1‑76 (MQ=255)
                         tCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAt                                               >  2:1700797/1‑76 (MQ=255)
                         tCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAt                                               <  2:872697/76‑1 (MQ=255)
                          ccAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtg                                              >  1:4227810/1‑76 (MQ=255)
                          ccAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtg                                              >  1:279723/1‑76 (MQ=255)
                           cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgc                                             >  1:923968/1‑76 (MQ=255)
                           cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgc                                             <  1:3200923/76‑1 (MQ=255)
                           cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgc                                             >  1:1880073/1‑76 (MQ=255)
                            aGCGGGGAGGTAAGCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgct                                            <  1:828784/76‑1 (MQ=255)
                            aGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgct                                            <  1:654917/76‑1 (MQ=255)
                             gCGGGGGGGTAACCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAActgctg                                           <  2:1766142/76‑1 (MQ=255)
                             gCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctg                                           >  1:588415/1‑76 (MQ=255)
                             gCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCCCCGACTCTCGATATGCCTTCACAAAtgctg                                           >  1:1124683/1‑76 (MQ=255)
                               ggggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgct                                         >  2:3531126/1‑76 (MQ=255)
                               ggggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgct                                         <  2:3013067/76‑1 (MQ=255)
                                gggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctg                                        <  1:3582488/76‑1 (MQ=255)
                                 gggggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                       <  2:984374/73‑1 (MQ=255)
                                  gAGGTGAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcc                                      >  2:2734614/1‑76 (MQ=255)
                                  gAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcc                                      <  1:2865368/76‑1 (MQ=255)
                                   gggTAAACAATACGGGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcct                                     <  2:1480749/75‑1 (MQ=255)
                                   aGGTAAACGATACGTGGTATTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAACTGctgctgcct                                     <  1:1536238/76‑1 (MQ=255)
                                    ggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctg                                    <  2:3318959/76‑1 (MQ=255)
                                    ggTAAACGATACGTGGTAATTTGTGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctg                                    <  1:2457881/76‑1 (MQ=255)
                                      tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                  <  2:3995072/76‑1 (MQ=255)
                                      gaaaCGATACGTGGTAATTTTGGGGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                  <  2:87231/75‑1 (MQ=255)
                                       aagcGATAAGTGGGTGTTTTGGTGGTTTTGATTGAACTGGCGCTAGATACGACTTAACAAATGctgctgcctgctg                                 <  2:189167/73‑1 (MQ=255)
                                       aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctg                                 >  2:3135566/1‑76 (MQ=255)
                                        aaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctgc                                <  2:2258973/76‑1 (MQ=255)
                                        aaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctgc                                <  1:1430859/76‑1 (MQ=255)
                                         aCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcg                               <  1:4813505/76‑1 (MQ=255)
                                          cGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcgc                              <  1:4207191/76‑1 (MQ=255)
                                          cGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcgc                              <  1:2724254/76‑1 (MQ=255)
                                           gATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGcgcg                             >  2:4721284/1‑76 (MQ=255)
                                            aTACGTTGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCgg                            <  2:1135848/76‑1 (MQ=255)
                                            aTACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCgg                            >  2:1495925/1‑76 (MQ=255)
                                             tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGt                           <  2:2217992/76‑1 (MQ=255)
                                             tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGt                           <  2:4787603/76‑1 (MQ=255)
                                               cGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgc                         >  2:329524/1‑76 (MQ=255)
                                                gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcg                        >  1:536492/1‑76 (MQ=255)
                                                gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcg                        >  1:662025/1‑76 (MQ=255)
                                                  ggTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGct                      <  1:1862232/76‑1 (MQ=255)
                                                   gTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGctc                     <  1:4461876/76‑1 (MQ=255)
                                                    tAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGctct                    >  2:286300/1‑76 (MQ=255)
                                                     aaTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTa                   >  1:2488481/1‑76 (MQ=255)
                                                      aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                  <  1:126526/76‑1 (MQ=255)
                                                      aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                  >  1:4552425/1‑76 (MQ=255)
                                                      aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                  <  1:1464971/76‑1 (MQ=255)
                                                       ttttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACg                 <  1:1350715/76‑1 (MQ=255)
                                                        tttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACgg                >  2:3948138/1‑76 (MQ=255)
                                                        tttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACgg                <  1:4030523/76‑1 (MQ=255)
                                                          tGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCCACCTGCTGCGCGGTGCGCTCTACGGAt              >  2:1323571/1‑76 (MQ=255)
                                                           ggTGATTTTCATTGCACTGCCTCGCGATACGCCTTAACAAATGCTGCTGCCTGCTGCTCGGTGCGCTCTACGGAtt             <  2:3189537/76‑1 (MQ=255)
                                                            gTGGTTTTCATTGCCGTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTg            <  2:4318285/76‑1 (MQ=255)
                                                            gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTg            >  1:2514752/1‑76 (MQ=255)
                                                            gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTg            <  1:2994720/76‑1 (MQ=255)
                                                            gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTg            >  2:1523042/1‑76 (MQ=255)
                                                            gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTg            >  1:3546500/1‑76 (MQ=255)
                                                              gATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGcc          <  2:3820764/76‑1 (MQ=255)
                                                              gATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGcc          >  1:4580368/1‑76 (MQ=255)
                                                                 tttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgg       >  1:3966843/1‑76 (MQ=255)
                                                                 tttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgg       >  1:4263434/1‑76 (MQ=255)
                                                                  ttCATGGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgc      <  2:3607071/76‑1 (MQ=255)
                                                                   tCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcgc                       <  1:3349063/58‑1 (MQ=255)
                                                                   tCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcg     >  2:3971022/1‑76 (MQ=255)
                                                                    cATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc    <  1:2239050/76‑1 (MQ=255)
                                                                    cATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc    <  2:966173/76‑1 (MQ=255)
                                                                     aTTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGcgcg   >  2:2180269/1‑76 (MQ=255)
                                                                      ttGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGa  >  2:3363391/1‑76 (MQ=255)
                                                                      ttGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGa  >  1:4412315/1‑76 (MQ=255)
                                                                           |                                                                       
TCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGA  >  NC_000913/3870555‑3870701

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