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,7890GA94.7% 108.3 / ‑2.3 38R401C (CGT→TGT) bglXbeta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (17/19);  total (18/20)
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

GTCAGTACGGTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACGGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAACCCC  >  NC_000913/2220659‑2220922
                                                                                                                                  |                                                                                                                                     
gTCAGTACGGTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCa                                                                                                                               >  2:134686/1‑139 (MQ=255)
   aGTACGGTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCa                                                                                                                            <  1:309563/139‑1 (MQ=255)
                      cGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCgcgcgc                                                                                                         <  2:529954/139‑1 (MQ=255)
                               cGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcg                                                                                                <  2:833965/139‑1 (MQ=255)
                                 gCTGCGGACCAGCTGCCCATCACTTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcgcg                                                                                              <  1:203416/139‑1 (MQ=255)
                                 gCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcgcg                                                                                              >  2:513899/1‑139 (MQ=255)
                                  cTGCGGACCAGCTGCCCATCGCGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACgcgc                                                                                             <  2:489033/139‑1 (MQ=255)
                                             cTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGt                                                                                  >  2:655871/1‑139 (MQ=255)
                                                              tttACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCgg                                                                   <  1:108043/137‑1 (MQ=255)
                                                              tttACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCgg                                                                   >  2:108043/1‑137 (MQ=255)
                                                                     tCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTAt                                                          >  1:152730/1‑139 (MQ=255)
                                                                       cGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATcc                                                        >  2:214756/1‑139 (MQ=255)
                                                                        gCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCa                                                       <  1:25150/139‑1 (MQ=255)
                                                                          cAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCTGGCGGCTTTCGGCATTGGTATCCAcc                                                     >  2:641408/1‑139 (MQ=255)
                                                                             tGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTCCGGTGCAGGCGGCTTTCGGCATTGGTATCCACCggg                                                  >  2:132616/1‑139 (MQ=255)
                                                                                    accaccGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGAct                                           >  2:395914/1‑139 (MQ=255)
                                                                                        ccGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCttt                                       >  1:183163/1‑139 (MQ=255)
                                                                                         cGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACGGTTTTTCAGCAACACCAAGCTGTCGCGCGCCACTTCACGCGCTTCTTTACGGTACAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTc                                      >  2:604820/1‑139 (MQ=255)
                                                                                         cGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTcg                                                                                                              <  1:428156/67‑1 (MQ=255)
                                                                                         cGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTcg                                                                                                              >  2:428156/1‑67 (MQ=255)
                                                                                         cGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTc                                      >  1:192395/1‑139 (MQ=255)
                                                                                           cAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCgg                                    <  1:699104/139‑1 (MQ=255)
                                                                                           cAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCgg                                    >  1:98114/1‑139 (MQ=255)
                                                                                                gtggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcg                                                                                                <  2:930507/74‑1 (MQ=255)
                                                                                                gtggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcg                                                                                                >  1:930507/1‑74 (MQ=255)
                                                                                                gtggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCa                               <  2:183163/139‑1 (MQ=255)
                                                                                                  ggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCaaa                             <  2:182484/139‑1 (MQ=255)
                                                                                                   gCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCa                                                                 >  2:427160/1‑102 (MQ=255)
                                                                                                   gCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCa                                                                 <  1:427160/102‑1 (MQ=255)
                                                                                                   gCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAAt                            >  1:444053/1‑139 (MQ=255)
                                                                                                           tttCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATg                   <  1:134686/139‑1 (MQ=255)
                                                                                                               aGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGt                <  1:214756/139‑1 (MQ=255)
                                                                                                                         tttCGAGACGGATTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCATTAaacaa      <  1:604820/139‑1 (MQ=255)
                                                                                                                           tCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAggcg                                                                           <  1:614228/68‑1 (MQ=255)
                                                                                                                           tCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAggcg                                                                           >  2:614228/1‑68 (MQ=255)
                                                                                                                           tCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAAcc    <  2:192395/139‑1 (MQ=255)
                                                                                                                           tCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAAcc    <  1:655871/139‑1 (MQ=255)
                                                                                                                             gagaCAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAAcccc  <  2:152730/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                     
GTCAGTACGGTCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACGGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAACCCC  >  NC_000913/2220659‑2220922

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