Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,870,765 +G coding (476/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,7641.G99.2% 342.1 / ‑3.8 125coding (477/1764 nt)dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (58/66);  total (58/67)
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.94e-01

Only 100 of 131 total aligned reads displayed.

TTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTG  >  NC_000913/3870689‑3870834
                                                                             |                                                                      
ttGCCCGGCGCGATAGCAATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCt                                                                         <  1:731311/76‑1 (MQ=255)
 tGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTg                                                                       <  1:3591457/76‑1 (MQ=255)
 tGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTg                                                                       >  2:4546810/1‑76 (MQ=255)
  gCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTgg                                                                      >  1:255656/1‑75 (MQ=255)
   cccGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggc                                                                     >  1:470313/1‑74 (MQ=255)
    ccGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggcg                                                                    <  1:2304156/76‑4 (MQ=255)
       gcgcGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCa                                                                  >  1:2401292/1‑76 (MQ=255)
        gccgATAGAGGCCCCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAc                                                                 <  2:243631/74‑1 (MQ=255)
        gccgATAGAGATCGCGCCATGAGCCCGCCTCTGTGCAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAc                                                                 <  1:831113/74‑1 (MQ=255)
        cgcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAc                                                                 <  1:5630502/76‑1 (MQ=255)
         gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcg                                                                <  1:2428346/76‑1 (MQ=255)
         gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcg                                                                >  1:1071826/1‑76 (MQ=255)
          tgATAGAGATCGGTTCCTAATCCCGCCCCTGCACAACCTCCGTCTATCCACGGCGCCAGGTTTTCTGGCGTCAcgc                                                               <  1:5050975/75‑1 (MQ=255)
          cgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGCTTTCTGGCGTCAcgc                                                               >  2:2937933/1‑76 (MQ=255)
           gATAGGGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgcc                                                              <  2:5082564/76‑1 (MQ=255)
           gATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgcc                                                              <  1:491380/76‑1 (MQ=255)
             tAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgccgc                                                            <  1:253251/76‑1 (MQ=255)
              agagATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACgccgcc                                                           <  1:2268710/76‑1 (MQ=255)
              agagATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACgccgcc                                                           <  1:1621073/76‑1 (MQ=255)
               gagaTCGCTGCCTACGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCa                                                          <  2:4434132/76‑1 (MQ=255)
               gagaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCa                                                          >  2:2806886/1‑76 (MQ=255)
               gagaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCa                                                          >  2:4853596/1‑76 (MQ=255)
                agaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCaa                                                         >  1:1691182/1‑76 (MQ=255)
                agaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTCCTGGCGTCACGCCGCCaa                                                         >  1:5192540/1‑76 (MQ=255)
                 gaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAAc                                                        >  2:635103/1‑76 (MQ=255)
                  aTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACg                                                       <  1:4870148/76‑1 (MQ=255)
                  aTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACg                                                       >  1:2158220/1‑76 (MQ=255)
                    cGCTGCCTACGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                     <  2:1146927/76‑1 (MQ=255)
                    cGCTCCCTACGCCCCCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                     <  1:1263148/76‑1 (MQ=255)
                       tGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCCCCAACGGCaaa                                                  >  2:766977/1‑76 (MQ=255)
                        gCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAg                                                 >  1:5698598/1‑76 (MQ=255)
                         ccTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGa                                                <  2:4589678/76‑1 (MQ=255)
                          cTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                               <  1:1015185/76‑1 (MQ=255)
                           tacgccggcccCTGCACAGCCTGCGTCTATCCACTGCGCCAGGTTTTCGGGCGTCACGCCGCCAACGGCAAAGACt                                              <  2:508290/73‑1 (MQ=255)
                           tAAGCCCGCCCCTGCACACCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACt                                              <  2:867254/76‑1 (MQ=255)
                           tAAGCCCGCCCCTGCACACCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACt                                              <  1:1258814/76‑1 (MQ=255)
                             accccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                            <  1:5237484/74‑1 (MQ=255)
                             accccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                            <  1:1018905/74‑1 (MQ=255)
                             aGCCCGCCCCTGCACACCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                            <  1:3564262/76‑1 (MQ=255)
                             aGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                            <  2:5511589/76‑1 (MQ=255)
                              gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCg                                           >  1:5363068/1‑76 (MQ=255)
                                ccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                            <  2:2357273/73‑1 (MQ=255)
                                ccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                            >  1:2357273/1‑73 (MQ=255)
                                ccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGAt                                         >  2:4566318/1‑76 (MQ=255)
                                ccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGAt                                         <  1:2343182/76‑1 (MQ=255)
                                ccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGAt                                         <  2:1933154/76‑1 (MQ=255)
                                 cgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCACAGACTCCGATg                                        >  2:5260132/1‑76 (MQ=255)
                                  gcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                       >  2:4726132/1‑76 (MQ=255)
                                  gcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                       >  1:3332671/1‑76 (MQ=255)
                                  gcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                       >  1:3272640/1‑76 (MQ=255)
                                  gcccCTCCACATCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                       <  1:5016922/76‑1 (MQ=255)
                                   ccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCCACGGCAACGACTGCGATgcc                                      >  2:227879/1‑74 (MQ=255)
                                    cccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                     <  2:1733024/76‑1 (MQ=255)
                                       tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGag                                  >  1:810253/1‑75 (MQ=255)
                                        gcccaACCGGCGTCTACCCACTGCGCCAGGTTTTCTGGCGGCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                 <  2:5385071/73‑1 (MQ=255)
                                        gCACACGCTGCGTCTATCCTCTGCGCAGGATTTTCGGGCTCCACGCCGCCAACGGCAAAGACTTCGATGTCCGATg                                 <  1:5219467/76‑1 (MQ=255)
                                        gCACAACCGGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                 <  2:3394060/76‑1 (MQ=255)
                                         cacaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATgg                                <  1:3705313/76‑1 (MQ=255)
                                          acaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGc                               >  2:2211461/1‑76 (MQ=255)
                                            atccTGCGTCTGGCCCCTGGGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCaa                             <  2:671032/74‑1 (MQ=255)
                                             aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAAt                            <  1:2333255/76‑1 (MQ=255)
                                              ccTGCGTCTGTCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGGCCGATGGCAATa                           <  1:4800355/76‑1 (MQ=255)
                                              ccTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCACAGACTGCGATGTCCGATGGCAATa                           >  1:1179477/1‑76 (MQ=255)
                                                 gCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCg                        <  1:3707530/76‑1 (MQ=255)
                                                 gCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCg                        >  1:4294265/1‑76 (MQ=255)
                                                 gCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCg                        <  1:4566319/76‑1 (MQ=255)
                                                   gTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                      >  1:2222771/1‑76 (MQ=255)
                                                   gTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                      <  2:2401292/76‑1 (MQ=255)
                                                   gTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                      <  1:3737933/76‑1 (MQ=255)
                                                     cTATACACTGCGCCAGGTTTTCTGGCGTAACGCCGCCAACGGCAAAGACTGCGATGGCCGATGGCAATACCGCttt                    >  1:2491287/1‑76 (MQ=255)
                                                      tagccACTGCGCGAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCtttt                   <  2:742742/73‑1 (MQ=255)
                                                       agccACGGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTACGATCGCAAGACCGACTTTa                  >  2:4268270/3‑76 (MQ=255)
                                                       aTCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAACGACTGCGATGTCCGATAGCCATCCCGCTTTTa                  >  2:604124/1‑76 (MQ=255)
                                                       aTCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                  <  1:4097674/76‑1 (MQ=255)
                                                        tCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                 >  1:3007504/1‑76 (MQ=255)
                                                        tCCACGGCGCCGGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                 <  2:4333812/76‑1 (MQ=255)
                                                         ccACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAAc                >  1:4286121/1‑76 (MQ=255)
                                                         ccACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCCTTTAAc                >  2:1238261/1‑76 (MQ=255)
                                                          cACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACg               >  2:4694729/1‑76 (MQ=255)
                                                          cACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACg               >  1:3826759/1‑76 (MQ=255)
                                                          cACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACg               <  2:4186274/76‑1 (MQ=255)
                                                          cACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACg               >  2:3453629/1‑76 (MQ=255)
                                                          cACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACg               <  1:5040694/76‑1 (MQ=255)
                                                           aCTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGc              >  1:3042330/1‑76 (MQ=255)
                                                           aCTGCGCCAGGGTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGCTGACAATACCGCTTTTAACGc              >  1:4319549/1‑76 (MQ=255)
                                                             tGCGCCAGGTTTTCTGGCGTTACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTTACGCtt            >  1:3118419/1‑76 (MQ=255)
                                                              gcgcCAGGTTTTCTGGTGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt           <  1:4554768/76‑1 (MQ=255)
                                                              gcgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt           >  1:68581/1‑76 (MQ=255)
                                                              gcgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt           >  2:297282/1‑76 (MQ=255)
                                                              gcgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt           <  1:4539934/76‑1 (MQ=255)
                                                               cgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTg          <  2:1718144/76‑1 (MQ=255)
                                                               cgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACCCTTTg          >  1:2372458/1‑76 (MQ=255)
                                                                  cAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATg       <  1:407257/76‑1 (MQ=255)
                                                                   gggTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGt      <  2:2793747/75‑1 (MQ=255)
                                                                   aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGt      >  2:5530307/1‑76 (MQ=255)
                                                                    ggTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTa     <  1:1583287/76‑1 (MQ=255)
                                                                    ggTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTa     >  1:1250291/1‑76 (MQ=255)
                                                                     gTTTTCTGGCGTCACGCCGCCAACGGCAAAAACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAt    <  2:4717534/76‑1 (MQ=255)
                                                                      ttttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt   <  2:2261922/76‑1 (MQ=255)
                                                                       tttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTg  >  1:96966/1‑76 (MQ=255)
                                                                             |                                                                      
TTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTG  >  NC_000913/3870689‑3870834

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