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.7% 88.1 / ‑4.9 30R401C (CGT→TGT) bglXbeta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (15/14);  total (15/15)
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

TCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACGGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAACCCC  >  NC_000913/2220669‑2220922
                                                                                                                        |                                                                                                                                     
tCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGTAACACCaa                                                                                                                     >  2:300978/1‑139 (MQ=255)
              gCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCAATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCa                                                                                                       >  2:512868/1‑139 (MQ=255)
                    ccGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAc                                                                                                 >  1:498303/1‑139 (MQ=255)
                                             cgtcacgtTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACAcc                                                                                                                       >  1:381733/1‑92 (MQ=255)
                                             cgtcacgtTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACAcc                                                                                                                       <  2:381733/92‑1 (MQ=255)
                                                  cgtTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCATTTCACGCGCTTCTTTACGGTGCAGGCGGCttcggc                                                                   >  1:506816/1‑135 (MQ=255)
                                                   ttttccggTCCGCCAGTGGCCCACCCACCGCATTGGTGGCCGTTTTTTTCAGCGTTAACTTTTCGGGACGGTTTTTCCGCACCACCAGCCTTTCCCGCACCACTTCACGCGCTTCTTTACGGGGCAGGCGGCTTCGGCa                                                                 <  2:506816/132‑1 (MQ=255)
                                                    tttACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGTAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCa                                                                 <  1:300978/139‑1 (MQ=255)
                                                    tttACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCa                                                                 <  2:498303/139‑1 (MQ=255)
                                                        cTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcgcg                                                                                              >  1:364534/1‑106 (MQ=255)
                                                        cTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcgcg                                                                                              <  2:747155/106‑1 (MQ=255)
                                                        cTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcgcg                                                                                              <  2:364534/106‑1 (MQ=255)
                                                        cTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcgcg                                                                                              >  1:747155/1‑106 (MQ=255)
                                                          gTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTa                                                           <  1:223338/139‑1 (MQ=255)
                                                               ccAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCAc                                                      <  1:496771/139‑1 (MQ=255)
                                                                cAGTGGCCCAACCATCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCTACTTCACGCGCTTCTTTACGGTGCAGCCGGCTTTCGGCATTGGTATCCAcc                                                     <  1:512868/139‑1 (MQ=255)
                                                                           ccaccGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTAc                                                                                     <  1:735073/96‑1 (MQ=255)
                                                                           ccaccGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTAc                                                                                     >  2:735073/1‑96 (MQ=255)
                                                                               cGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTc                                      >  2:634280/1‑139 (MQ=255)
                                                                                    tggtggCCGATTTTTTCAGCGGTAACGTTTAGAGACAGCTTTTCAGCCACCCCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCTGCATTGGTATCAACCGGGTCAGACTCTTTAGGccc                                 >  2:179477/1‑139 (MQ=255)
                                                                                        ggCCGATTTTTTCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCaaa                             >  2:762933/1‑139 (MQ=255)
                                                                                              ttttttCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcg                                                                                                >  2:545617/1‑66 (MQ=255)
                                                                                              ttttttCAGCGGTAACGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCAcg                                                                                                <  1:545617/66‑1 (MQ=255)
                                                                                                            aCGTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTaaa         >  1:215945/1‑139 (MQ=255)
                                                                                                              gTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGt                                                 <  2:390075/97‑1 (MQ=255)
                                                                                                              gTTTCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGt                                                 >  1:390075/1‑97 (MQ=255)
                                                                                                                 tCGAGACAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAAcc    <  1:762933/139‑1 (MQ=255)
                                                                                                                   gagaCAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGTCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAAcccc  <  2:514806/139‑1 (MQ=255)
                                                                                                                   gagaCAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATg                   >  1:827283/1‑122 (MQ=255)
                                                                                                                   gagaCAGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATg                   <  2:827283/122‑1 (MQ=255)
                                                                                                                        |                                                                                                                                     
TCACGGATTGATCGGCAACACCGGCTGCGGACCAGCTGCCCATCACGTCACGTTTACTGTCCGCCAGTGGCCCAACCACCGCAATGGTGGCCGATTTTTTCAGCGGTAACGTTTCGAGACGGTTTTTCAGCAACACCAAGCTTTCGCGCGCCACTTCACGCGCTTCTTTACGGTGCAGGCGGCTTTCGGCATTGGTATCCACCGGGTCAGACTCTTTCGGCCCCAAATGGCTGTATGGGTCGTTAAACAACCCC  >  NC_000913/2220669‑2220922

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