Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,220,789 G→A R401C (CGT→TGT)  bglX ← beta‑D‑glucoside glucohydrolase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,220,7890GA96.9% 199.1 / ‑1.9 65R401C (CGT→TGT) bglXbeta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (31/32);  total (32/33)
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.07e-01

GGTCAGTACGGTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACGGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAACCCCA  >  NC_000913/2220658‑2220923
                                                                                                                                   |                                                                                                                                      
ggTCAGTACGGTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTc                                                                                                                                 >  1:940387/1‑139 (MQ=255)
 gTCAGTACGGTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCa                                                                                                                                >  2:939571/1‑139 (MQ=255)
 gTCAGTACGGTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCa                                                                                                                                >  1:60531/1‑139 (MQ=255)
  tCAGTACGGTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAg                                                                                                                               <  1:991159/139‑1 (MQ=255)
      tACGGTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCNATGGTGNNCNANTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAac                                                                                                                           >  1:1348883/1‑139 (MQ=255)
         ggTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACAcc                                                                                                                        <  2:884789/139‑1 (MQ=255)
           tCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCaa                                                                                                                      >  2:801430/1‑139 (MQ=255)
           tCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCaa                                                                                                                      >  2:372070/1‑139 (MQ=255)
                     aTCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTcgcgc                                                                                                            >  1:231081/1‑139 (MQ=255)
                        ggCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAACTTTCGCGCGcc                                                                                                         <  1:939571/139‑1 (MQ=255)
                         gCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTcg                                                                                                               <  1:1252686/132‑1 (MQ=255)
                         gCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTcg                                                                                                               >  2:1252686/1‑132 (MQ=255)
                                cGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcg                                                                                                 <  1:76429/139‑1 (MQ=255)
                                  gCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcgcg                                                                                               >  1:1012406/1‑139 (MQ=255)
                                  gCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcgcg                                                                                               >  1:1266517/1‑139 (MQ=255)
                                                   cATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAgg                                                                              >  2:702151/1‑139 (MQ=255)
                                                     tcacgtcacgTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGATTCTTTACGGTGCAggcg                                                                            <  2:940387/139‑1 (MQ=255)
                                                        cgtcacgtTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCt                                                                         <  2:726820/139‑1 (MQ=255)
                                                          tcacgtTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCttt                                                                       <  2:1266517/139‑1 (MQ=255)
                                                           cacgtTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGc                                                                                                                              >  1:1375165/1‑83 (MQ=255)
                                                           cacgtTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGc                                                                                                                              <  2:1375165/83‑1 (MQ=255)
                                                           cacgtTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTc                                                                      >  2:415713/1‑139 (MQ=255)
                                                               tttACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCa                                                                  <  1:372070/139‑1 (MQ=255)
                                                               tttACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCa                                                                  <  2:1348883/139‑1 (MQ=255)
                                                                   cTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTgg                                                              >  1:299835/1‑139 (MQ=255)
                                                                     gTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGttt                                                                                                                                    <  1:1157789/67‑1 (MQ=255)
                                                                     gTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGttt                                                                                                                                    >  2:1157789/1‑67 (MQ=255)
                                                                      tCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcg                                                                                                 >  1:1185378/1‑101 (MQ=255)
                                                                      tCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcg                                                                                                 <  2:1185378/101‑1 (MQ=255)
                                                                               ggCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGt                                                  <  2:299835/139‑1 (MQ=255)
                                                                                   cAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGa                                              <  1:1653611/139‑1 (MQ=255)
                                                                                        accGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGc                                                                                                                    <  1:1654246/64‑1 (MQ=255)
                                                                                        accGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGc                                                                                                                    >  2:1654246/1‑64 (MQ=255)
                                                                                        accGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCtt                                         <  2:1012406/139‑1 (MQ=255)
                                                                                         ccGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACGGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCttt                                        >  1:1400294/1‑139 (MQ=255)
                                                                                         ccGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCttt                                        <  2:1473448/139‑1 (MQ=255)
                                                                                          cGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcgcg                                                                                               <  1:1244003/83‑1 (MQ=255)
                                                                                          cGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcgcg                                                                                               >  2:1244003/1‑83 (MQ=255)
                                                                                            cAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCgg                                     <  1:801430/139‑1 (MQ=255)
                                                                                              aTGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGcc                                   <  1:415713/139‑1 (MQ=255)
                                                                                               tggtggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCaa                                                                                                                            <  1:86239/49‑1 (MQ=255)
                                                                                               tggtggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCaa                                                                                                                            >  2:86239/1‑49 (MQ=255)
                                                                                               tggtggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCggg                                                   >  1:696578/1‑122 (MQ=255)
                                                                                               tggtggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCggg                                                   <  2:696578/122‑1 (MQ=255)
                                                                                               tggtggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGccc                                  <  2:470341/139‑1 (MQ=255)
                                                                                               tggtggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGccc                                  >  1:470341/1‑139 (MQ=255)
                                                                                                 gtggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCa                                <  2:231081/139‑1 (MQ=255)
                                                                                                  tggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCaa                               >  2:1448378/1‑139 (MQ=255)
                                                                                                  tggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCaa                               <  2:60531/139‑1 (MQ=255)
                                                                                                    gCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAAt                             >  1:1451928/1‑139 (MQ=255)
                                                                                                    gCCGATTTTTTCAGCGGTAACG‑TTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAgg                                                                              <  1:1558000/89‑1 (MQ=255)
                                                                                                    gCCGATTTTTTCAGCGGTAACG‑TTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAgg                                                                              >  2:1558000/1‑89 (MQ=255)
                                                                                                      cGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATgg                           <  1:1448378/139‑1 (MQ=255)
                                                                                                         ttttttCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTg                        >  1:885037/1‑139 (MQ=255)
                                                                                                            tttCAGCGGTAACGTTTCGAGACGGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTAt                     <  2:1400294/139‑1 (MQ=255)
                                                                                                                   ggTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTcgcg                                                                                                             >  2:1525169/1‑44 (MQ=255)
                                                                                                                   ggTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTcgcg                                                                                                             <  1:1525169/44‑1 (MQ=255)
                                                                                                                       aCGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCAt                                                                 >  2:247036/1‑84 (MQ=255)
                                                                                                                       aCGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCAt                                                                 <  1:247036/84‑1 (MQ=255)
                                                                                                                             cGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcg                                                                                                 <  1:735001/46‑1 (MQ=255)
                                                                                                                             cGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcg                                                                                                 >  2:735001/1‑46 (MQ=255)
                                                                                                                             cGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAAccc    >  2:1007376/1‑139 (MQ=255)
                                                                                                                              gagaCAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCttt                                                                       >  1:1347389/1‑71 (MQ=255)
                                                                                                                              gagaCAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCttt                                                                       <  2:1347389/71‑1 (MQ=255)
                                                                                                                               agaCAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAACCCca  >  1:1578874/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                      
GGTCAGTACGGTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACGGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAACCCCA  >  NC_000913/2220658‑2220923

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