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.C93.9% 108.2 / 0.1 33coding (920/1041 nt)galSgalactose‑ and fucose‑inducible galactose regulon transcriptional isorepressor, mgl operon transcriptional repressor, autorepressor
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (12/19);  total (13/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 = 9.96e-01

AGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑G‑CCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGA  >  NC_000913/2240615‑2240880
                                                                                                                                        |                                                                                                                                   
aGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑G‑cccc                                                                                                                                 <  1:381210/139‑1 (MQ=255)
  aTTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑G‑cccccc                                                                                                                               >  1:60359/1‑138 (MQ=255)
     gAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑G‑cccccctgc                                                                                                                            <  2:729907/139‑5 (MQ=255)
          cGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACTAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGc                                                                                                                        >  1:714494/1‑139 (MQ=255)
            gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCa                                                                                                                      >  2:966667/1‑139 (MQ=255)
                     ccTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTgg                                                                                                             >  1:845694/1‑139 (MQ=255)
                                   ttAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCAtt                                                                                               <  2:822055/139‑1 (MQ=255)
                                   ttAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCAtt                                                                                               >  2:118157/1‑139 (MQ=255)
                                   ttAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCAtt                                                                                               >  1:1013750/1‑139 (MQ=255)
                                    tAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTg                                                                                              >  1:1062409/1‑139 (MQ=255)
                                     aGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGa                                                                                             <  2:395603/139‑1 (MQ=255)
                                       tGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAg                                                                                           <  2:270477/139‑1 (MQ=255)
                                       tGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAg                                                                                           <  2:541756/139‑1 (MQ=255)
                                          tCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCaa                                                                                        <  2:845694/139‑1 (MQ=255)
                                          tCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCaa                                                                                        <  2:1062409/139‑1 (MQ=255)
                                           cgccgccgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAAt                                                                                       <  2:245446/139‑1 (MQ=255)
                                           cgccgccgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAAt                                                                                       <  2:817168/139‑1 (MQ=255)
                                               gccgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGa                                                                                   <  2:1013750/139‑1 (MQ=255)
                                               gccgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGa                                                                                   <  1:976696/139‑1 (MQ=255)
                                                      tCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  2:1108482/139‑1 (MQ=255)
                                                      tCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  1:809101/139‑1 (MQ=255)
                                                             tgtTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑G‑cccccctgc                                                                                                                            >  2:296289/1‑79 (MQ=255)
                                                             tgtTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑G‑cccccctgc                                                                                                                            <  1:296289/83‑5 (MQ=255)
                                                                        aaTGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGt                                                          >  2:1108187/1‑139 (MQ=255)
                                                                               aCGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTa                                                   >  2:223624/1‑139 (MQ=255)
                                                                                     tAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGc                                             <  1:223624/139‑1 (MQ=255)
                                                                                                        tggctggcACGAGGATCAATATTGCCTGCTAC‑CCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGa                                                                                             >  2:774214/1‑72 (MQ=255)
                                                                                                        tggctggcACGAGGATCAATATTGCCTGCTAC‑CCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGa                                                                                             <  1:774214/72‑1 (MQ=255)
                                                                                                         ggctggcACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAAc                         <  2:290743/139‑1 (MQ=255)
                                                                                                              gcACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTc                                                         <  2:1014016/102‑1 (MQ=255)
                                                                                                              gcACGAGGATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTc                                                         >  1:1014016/1‑102 (MQ=255)
                                                                                                                    ggATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTgag              <  2:60359/139‑1 (MQ=255)
                                                                                                                    ggATCAATATTGCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTgag              <  1:1108187/139‑1 (MQ=255)
                                                                                                                            atcgCCTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGCGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTaa      <  2:498245/136‑1 (MQ=255)
                                                                                                                                ccTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGa                                                                                             >  1:779777/1‑48 (MQ=255)
                                                                                                                                ccTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGa                                                                                             <  2:779777/48‑1 (MQ=255)
                                                                                                                                ccTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGc                                             <  1:85281/96‑1 (MQ=255)
                                                                                                                                ccTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGc                                             >  2:85281/1‑96 (MQ=255)
                                                                                                                                ccTGCT‑GCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGa  >  1:966066/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                   
AGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCT‑G‑CCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGA  >  NC_000913/2240615‑2240880

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