Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,896,219 A→G T281T (ACT→ACC)  bglB ← 6‑phospho‑beta‑glucosidase BglB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,896,2190AG96.3% 178.4 / ‑0.5 55T281T (ACT→ACC) bglB6‑phospho‑beta‑glucosidase BglB
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (24/28);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGCAACATATCCT  >  NZ_CP009273/3896084‑3896346
                                                                                                                                       |                                                                                                                               
cGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGt                                                                                                                               <  1:604380/138‑1 (MQ=255)
cGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGt                                                                                                                               <  2:673666/138‑1 (MQ=255)
      ggggATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCAtctc                                                                                                                         <  2:395319/138‑1 (MQ=255)
       gggATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCaa                                                                                                                       <  2:462209/139‑1 (MQ=255)
       gggATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCaa                                                                                                                       <  2:596845/139‑1 (MQ=255)
             atatTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGa                                                                                                                 <  1:600573/139‑1 (MQ=255)
                   aGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATtgt                                                                                                           >  2:471876/1‑139 (MQ=255)
                     tATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTgg                                                                                                         >  2:307383/1‑139 (MQ=255)
                     tATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTgg                                                                                                         >  2:432167/1‑139 (MQ=255)
                     tATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTgg                                                                                                         >  2:441262/1‑139 (MQ=255)
                         ttGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTcgc                                                                                                     <  1:307383/139‑1 (MQ=255)
                          tGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTcgcg                                                                                                    >  1:389342/1‑139 (MQ=255)
                          tGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTcgcg                                                                                                    >  2:510677/1‑139 (MQ=255)
                               tGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGt                                                                                                                               >  2:615373/1‑107 (MQ=255)
                               tGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGt                                                                                                                               <  1:615373/107‑1 (MQ=255)
                                 cgcATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGaaa                                                                                             >  2:208384/1‑139 (MQ=255)
                                   cATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACg                                                                                           <  2:794701/139‑1 (MQ=255)
                                     tttttATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCt                                                                                         <  2:735747/139‑1 (MQ=255)
                                           ttaATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCatat                                                                                   <  1:510677/139‑1 (MQ=255)
                                           ttaATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCatat                                                                                   <  1:525832/139‑1 (MQ=255)
                                                gCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCAt                                                                                                                            <  1:91797/93‑1 (MQ=255)
                                                gCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCAt                                                                                                                            >  2:91797/1‑93 (MQ=255)
                                                gCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCg                                                                              >  2:785252/1‑139 (MQ=255)
                                                 cTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCg                                                                              <  2:33349/138‑1 (MQ=255)
                                                                    aaCCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCt                                                           >  1:42798/1‑138 (MQ=255)
                                                                      ccAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGt                                                                                                                               >  2:205989/1‑68 (MQ=255)
                                                                      ccAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGt                                                                                                                               <  1:205989/68‑1 (MQ=255)
                                                                      ccAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGaa                                                        <  2:409823/139‑1 (MQ=255)
                                                                         gTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACAt                                                     <  1:826697/139‑1 (MQ=255)
                                                                              gTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCa                                                                  <  1:103241/121‑1 (MQ=255)
                                                                              gTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCa                                                                  >  2:103241/1‑121 (MQ=255)
                                                                              gTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCAcc                                                 >  2:865367/1‑138 (MQ=255)
                                                                               taataaCTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCa                                                >  2:665718/1‑138 (MQ=255)
                                                                                 ataaCTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCaaa                                              <  1:208384/138‑1 (MQ=255)
                                                                                  taaCTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGa                                            <  1:240667/139‑1 (MQ=255)
                                                                                    aCTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGt                                                                                                                               <  1:35795/54‑1 (MQ=255)
                                                                                    aCTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGt                                                                                                                               >  2:35795/1‑54 (MQ=255)
                                                                                    aCTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGAt                                                                          <  2:14400/107‑1 (MQ=255)
                                                                                    aCTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGAt                                                                          >  1:14400/1‑107 (MQ=255)
                                                                                    aCTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGaa                                           >  2:61051/1‑138 (MQ=255)
                                                                                         aaGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAg                                                                                 <  1:815843/95‑1 (MQ=255)
                                                                                         aaGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAg                                                                                 >  2:815843/1‑95 (MQ=255)
                                                                                                      gACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTctct                        >  2:709759/1‑139 (MQ=255)
                                                                                                      gACGGTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTctct                        >  2:841066/1‑139 (MQ=255)
                                                                                                          gTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCa                    <  1:665718/139‑1 (MQ=255)
                                                                                                          gTATGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCa                    <  2:16117/139‑1 (MQ=255)
                                                                                                             tGTTTTAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATg                  <  2:839081/138‑1 (MQ=255)
                                                                                                               ttttAAATCTTCTGCGTCACTTTCTGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCGGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGc                <  2:703344/138‑1 (MQ=255)
                                                                                                               ttttAAATCTGCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGc                <  2:591203/138‑1 (MQ=255)
                                                                                                                 ttAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTg             >  2:824297/1‑139 (MQ=255)
                                                                                                                  tAAATCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTg             >  2:593435/1‑138 (MQ=255)
                                                                                                                   aaaTCTTCTGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGc            <  1:61051/138‑1 (MQ=255)
                                                                                                                       cttctGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCat                                                                                     >  1:552928/1‑61 (MQ=255)
                                                                                                                       cttctGCGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCat                                                                                     <  2:552928/61‑1 (MQ=255)
                                                                                                                             cGTCACTTTCGGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGCAACATATCCt  <  2:257326/138‑1 (MQ=255)
                                                                                                                                       |                                                                                                                               
CGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGCAACATATCCT  >  NZ_CP009273/3896084‑3896346

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