Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,818,550 (C)7→8 intergenic (‑50/+55) chbF ← / ← chbR phospho‑chitobiase, general 6‑phospho‑beta‑glucosidase activity/repressor of chb operon for N,N'‑diacetylchitobiose utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,818,5431.C98.6% 243.6 / ‑3.7 69intergenic (‑43/+62)chbF/chbRphospho‑chitobiase, general 6‑phospho‑beta‑glucosidase activity/repressor of chb operon for N,N'‑diacetylchitobiose utilization
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (34/34);  total (34/35)
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.98e-01

GTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTA  >  NC_000913/1818414‑1818668
                                                                                                                                  |                                                                                                                               
gTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccccccg                                                                                                                      <  2:1542956/139‑3 (MQ=255)
                  tCCTTCCAGTAACTCCGGGGTATAGCTGCTCCTGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCa                                                                                                     >  1:1548169/1‑139 (MQ=255)
                  tCCTTCCAGTAACTACGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCa                                                                                                     >  2:493668/1‑139 (MQ=255)
                    cTTCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                                                                   <  2:1804490/139‑1 (MQ=255)
                      tCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAg                                                                                                 >  2:1076746/1‑139 (MQ=255)
                        cAGTAACTCCGGGGTATAGCTGCTCCTGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑CCCACCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  2:1548169/139‑1 (MQ=255)
                        cAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  1:1076746/139‑1 (MQ=255)
                        cAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  2:1123914/139‑1 (MQ=255)
                        cAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  2:1598875/139‑1 (MQ=255)
                        cAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  2:1805913/139‑1 (MQ=255)
                        cAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  2:1895382/139‑1 (MQ=255)
                                 cGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATAtgt                                                                                      >  1:1370405/1‑139 (MQ=255)
                                         aGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGt                                                                               <  2:1182912/139‑1 (MQ=255)
                                         aGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                                                                    <  2:972234/117‑1 (MQ=255)
                                         aGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                                                                    >  1:972234/1‑117 (MQ=255)
                                         aGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccccccggac                                                                                                                   <  2:1471896/101‑6 (MQ=255)
                                         aGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATAACTGGCATATCTG‑‑ccccccccggac                                                                                                                   >  1:1471896/1‑96 (MQ=255)
                                          gctgctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGAACAGATATGTGAATTGTCa                                                                             <  2:1370405/139‑1 (MQ=255)
                                             gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                                                              <  1:330090/119‑1 (MQ=255)
                                             gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                                                              >  2:330090/1‑119 (MQ=255)
                                              ctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACACAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGta                                                                         >  1:548389/1‑139 (MQ=255)
                                                cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                                                         >  2:638431/1‑105 (MQ=255)
                                                cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                                                         <  1:638431/105‑1 (MQ=255)
                                                cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGa                                                                                           <  2:884122/119‑1 (MQ=255)
                                                cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGa                                                                                           >  1:884122/1‑119 (MQ=255)
                                                                acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                                                             <  1:2008223/117‑1 (MQ=255)
                                                                acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                                                             >  2:2008223/1‑117 (MQ=255)
                                                                   acTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGt                                                                               <  1:1931059/113‑1 (MQ=255)
                                                                   acTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGt                                                                               >  2:1931059/1‑113 (MQ=255)
                                                                   acTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                                                                    >  1:1333212/1‑91 (MQ=255)
                                                                   acTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                                                                    <  2:1333212/91‑1 (MQ=255)
                                                                          aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACACAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCaa                                             <  2:548389/139‑1 (MQ=255)
                                                                              cTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAATAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTc                                         >  1:11338/1‑139 (MQ=255)
                                                                              cTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cccccc                                                                                                                         <  2:1705995/58‑1 (MQ=255)
                                                                              cTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cccccc                                                                                                                         >  1:1705995/1‑58 (MQ=255)
                                                                               tGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               >  1:1356619/1‑84 (MQ=255)
                                                                               tGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  2:1356619/84‑1 (MQ=255)
                                                                               tGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCt                                        >  2:1409178/1‑139 (MQ=255)
                                                                               tGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccccc                                                                                                                        >  1:887936/1‑58 (MQ=255)
                                                                               tGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATAGCTGGCATATCTG‑‑ccccccc                                                                                                                        <  2:887936/58‑1 (MQ=255)
                                                                                       aaTTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGAAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAg                                <  1:1409178/139‑1 (MQ=255)
                                                                                             tcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                                                                    <  2:430750/65‑1 (MQ=255)
                                                                                             tcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                                                                    >  1:430750/1‑65 (MQ=255)
                                                                                                 ttCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGAGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCg                      >  2:1339732/1‑139 (MQ=255)
                                                                                                  tCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                                                         <  2:46736/55‑1 (MQ=255)
                                                                                                  tCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                                                         >  1:46736/1‑55 (MQ=255)
                                                                                                  tCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               >  1:818237/1‑65 (MQ=255)
                                                                                                  tCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  2:818237/65‑1 (MQ=255)
                                                                                                    ggTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGAGATAAATAATCCAGCAACAGGa                                                                                               >  2:1263316/2‑63 (MQ=255)
                                                                                                    aGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  1:1263316/63‑1 (MQ=255)
                                                                                                      tACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGAt                                  >  1:2033209/1‑122 (MQ=255)
                                                                                                      tACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGAt                                  <  2:2033209/122‑1 (MQ=255)
                                                                                                           aaTACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAg                                                                                                 >  1:183509/1‑54 (MQ=255)
                                                                                                           aaTACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAg                                                                                                 <  2:183509/54‑1 (MQ=255)
                                                                                                           aaTACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGt            >  1:295036/1‑139 (MQ=255)
                                                                                                              aCTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  1:803845/53‑1 (MQ=255)
                                                                                                              aCTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               >  2:803845/1‑53 (MQ=255)
                                                                                                              aCTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTaa         >  2:1003463/1‑139 (MQ=255)
                                                                                                                tGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGa                                   >  1:1527366/1‑111 (MQ=255)
                                                                                                                tGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGa                                   <  2:1527366/111‑1 (MQ=255)
                                                                                                                tGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCtt                           <  1:1596829/119‑1 (MQ=255)
                                                                                                                tGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCtt                           >  2:1596829/1‑119 (MQ=255)
                                                                                                                tGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAAtt       >  1:259548/1‑139 (MQ=255)
                                                                                                                  atatCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                                                                    <  1:1382188/44‑1 (MQ=255)
                                                                                                                  atatCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                                                                    >  2:1382188/1‑44 (MQ=255)
                                                                                                                  atatCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATAtgt                                                                                      >  1:1886103/1‑58 (MQ=255)
                                                                                                                  atatCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATAtgt                                                                                      <  2:1886103/58‑1 (MQ=255)
                                                                                                                  atatCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                                                             >  1:1129012/1‑67 (MQ=255)
                                                                                                                  atatCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                                                             <  2:1129012/67‑1 (MQ=255)
                                                                                                                     tCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACGGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTa  <  2:11338/139‑1 (MQ=255)
                                                                                                                     tCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               >  2:1450976/1‑46 (MQ=255)
                                                                                                                     tCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  1:1450976/46‑1 (MQ=255)
                                                                                                                       tGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                                                                      <  2:59780/37‑1 (MQ=255)
                                                                                                                       tGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                                                                      >  1:59780/1‑37 (MQ=255)
                                                                                                                       tGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               <  2:749020/44‑1 (MQ=255)
                                                                                                                       tGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                                                               >  1:749020/1‑44 (MQ=255)
                                                                                                                                  |                                                                                                                               
GTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGGCGTAAAGGACGTTAATTTCTTA  >  NC_000913/1818414‑1818668

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: