Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,240,754 (C)5→6 coding (915/1041 nt) galS ← galactose‑ and fucose‑inducible galactose regulon transcriptional isorepressor, mgl operon transcriptional repressor, autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,240,7491.C98.1% 180.7 / ‑3.5 52coding (920/1041 nt)galSgalactose‑ and fucose‑inducible galactose regulon transcriptional isorepressor, mgl operon transcriptional repressor, autorepressor
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (23/28);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CACAGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGA  >  NC_000913/2240612‑2240880
                                                                                                                                          |                                                                                                                                    
cacaGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCtgctgc                                                                                                                                     >  2:1125548/1‑139 (MQ=255)
       tGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctg                                                                                                                             <  2:117974/139‑4 (MQ=255)
             cGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGc                                                                                                                        <  2:436756/139‑1 (MQ=255)
              ggTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGcc                                                                                                                       >  2:1347665/1‑139 (MQ=255)
                                        aGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGa                                                                                             >  2:460226/1‑139 (MQ=255)
                                             tCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCaa                                                                                        <  2:619759/139‑1 (MQ=255)
                                                  gccgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGTAGCAATGGGa                                                                                   <  1:810214/139‑1 (MQ=255)
                                                   ccgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGAt                                                                                  >  1:56275/1‑139 (MQ=255)
                                                    cgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATa                                                                                 >  1:23228/1‑139 (MQ=255)
                                                         tCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  1:476268/139‑1 (MQ=255)
                                                         tCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  2:1151628/139‑1 (MQ=255)
                                                         tCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  2:23228/139‑1 (MQ=255)
                                                         tCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  2:56275/139‑1 (MQ=255)
                                                           tGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACg                                                                          <  1:182651/139‑1 (MQ=255)
                                                           tGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACg                                                                          >  1:684404/1‑139 (MQ=255)
                                                                 gtTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGtt                                                                    <  1:1347665/139‑1 (MQ=255)
                                                                  ttGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTa                                                                   >  1:249434/1‑139 (MQ=255)
                                                                   tGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTaa                                                                  >  2:463467/1‑139 (MQ=255)
                                                                     cGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTCCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAAtt                                                                >  1:982871/1‑139 (MQ=255)
                                                                           aaTGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGa                                                                                             >  2:1015294/1‑104 (MQ=255)
                                                                           aaTGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGa                                                                                             <  1:1015294/104‑1 (MQ=255)
                                                                            aTGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTc                                                         <  2:1121094/139‑1 (MQ=255)
                                                                             tGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCg                                                                                                                 <  2:27691/82‑1 (MQ=255)
                                                                             tGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCg                                                                                                                 >  1:27691/1‑82 (MQ=255)
                                                                                   cGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctgc                                                                                                                            >  1:134058/1‑60 (MQ=255)
                                                                                   cGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctgc                                                                                                                            <  2:134058/64‑5 (MQ=255)
                                                                                     cACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAAttt                                                                                                      >  2:688595/1‑85 (MQ=255)
                                                                                     cACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAAttt                                                                                                      <  1:688595/85‑1 (MQ=255)
                                                                                     cACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            >  2:284727/1‑111 (MQ=255)
                                                                                     cACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  1:284727/111‑1 (MQ=255)
                                                                                      aCTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACgg                                               <  1:338784/139‑1 (MQ=255)
                                                                                       cTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACggg                                              <  1:463467/139‑1 (MQ=255)
                                                                                              cGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccc                                                                                                                                 >  1:1153214/1‑48 (MQ=255)
                                                                                              cGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccc                                                                                                                                 <  2:1153214/48‑1 (MQ=255)
                                                                                               ggCATAAAACAGTGGCTGGGACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGc                                             <  1:1012690/132‑1 (MQ=255)
                                                                                               ggCATAAAACAGTGGCTGGGACGAGGATCAATATTGCCTGCTGCACCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGc                                             >  2:1012690/1‑132 (MQ=255)
                                                                                               ggCATAAAACAGTGGCTGGCACGAGGATCAATATTGCctgctg                                                                                                                                       >  2:804722/1‑43 (MQ=255)
                                                                                               ggCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTc                                                                                                     <  1:477243/76‑1 (MQ=255)
                                                                                               ggCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTc                                                                                                     >  2:477243/1‑76 (MQ=255)
                                                                                               ggCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTc                                                         <  2:1287817/120‑1 (MQ=255)
                                                                                               ggCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTc                                                         >  1:1287817/1‑120 (MQ=255)
                                                                                               ggCATAAAACAGTGGCTGGCACAAGGATCAATATTGCctgctg                                                                                                                                       <  1:804722/43‑1 (MQ=255)
                                                                                                gCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCg                                                                                                                 >  1:300540/1‑63 (MQ=255)
                                                                                                gCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCg                                                                                                                 <  2:300540/63‑1 (MQ=255)
                                                                                                 cATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAat                                    <  1:379821/139‑1 (MQ=255)
                                                                                                           tggctggcACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGaaa                          >  2:451845/1‑139 (MQ=255)
                                                                                                              ctggcACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCg                       >  2:587862/1‑139 (MQ=255)
                                                                                                               tggcACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCga                      <  2:249434/139‑1 (MQ=255)
                                                                                                                  cACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCgatga                   >  2:43036/1‑139 (MQ=255)
                                                                                                                       ggATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTgag              <  2:484810/139‑1 (MQ=255)
                                                                                                                          tCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTgagaga           <  1:460226/139‑1 (MQ=255)
                                                                                                                             atatTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGAt                                          >  1:168523/1‑105 (MQ=255)
                                                                                                                             atatTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGAt                                          <  2:168523/105‑1 (MQ=255)
                                                                                                                                 tGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTaa      >  1:1158015/1‑137 (MQ=255)
                                                                                                                                 tGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTaa      <  2:1158015/137‑1 (MQ=255)
                                                                                                                                  gCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAg                                                                                           >  1:1221737/1‑51 (MQ=255)
                                                                                                                                  gCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAg                                                                                           <  2:1221737/51‑1 (MQ=255)
                                                                                                                                  gCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  2:138537/66‑1 (MQ=255)
                                                                                                                                  gCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            >  1:138537/1‑66 (MQ=255)
                                                                                                                                   ccTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGa  >  1:874019/1‑139 (MQ=255)
                                                                                                                                    ctgctgcCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            >  2:635084/1‑64 (MQ=255)
                                                                                                                                    ctgctgcCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  1:635084/64‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                    
CACAGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGA  >  NC_000913/2240612‑2240880

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