Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑1/+1) gatC ← / ← gatC pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.94.1% 219.5 / 3.8 51pseudogene (2/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (24/24);  total (25/26)
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
*NC_0009132,173,3620C.96.0% 219.1 / 0.2 51pseudogene (1/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (24/24);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATT  >  NC_000913/2173220‑2173502
                                                                                                                                              ||                                                                                                                                            
gATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgta                                                                                                                                        >  2:816381/1‑143 (MQ=255)
gATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgta                                                                                                                                        <  2:335393/149‑7 (MQ=255)
 aTTTCCACGATGCACGGCGACCGCCATCGCCCAGCCAGAACGGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                                       >  1:906673/1‑142 (MQ=255)
      cACGATGCACGGCGACCGCCATCGCCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                                                                                                                >  1:780985/1‑149 (MQ=255)
       aCGATGCACGGCGACCGCCATCGCCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACGCAAAGAGAAATTACAACGGGTATTGTGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                               >  1:860849/1‑149 (MQ=255)
        ggATGCAGGGCGAGCGCCATGGCCACCAAGAAGCCTATGGTGGCAAAATCGCCAAAGGGCAGTACCTAATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                              <  2:815868/148‑1 (MQ=255)
        cGATGCACGGCGACCGCCATCGCCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                              <  1:19489/149‑1 (MQ=255)
              aCGGCGACCGCCATCGCCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCaa                                                                                                                        >  2:665169/1‑149 (MQ=255)
              aCGGCGACCGCCATCGCCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCaa                                                                                                                        >  2:664936/1‑149 (MQ=255)
                   gACCGCCATCGCCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                                                                                                   <  2:846301/149‑1 (MQ=255)
                    aCCGCCATCGCCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                  <  2:780985/149‑1 (MQ=255)
                          aTCGCCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATACAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                            >  1:816445/1‑149 (MQ=255)
                               ccACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACACCATTGTGCTGCACCAACATCGCAGTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                        >  1:450538/1‑149 (MQ=255)
                                cACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                       >  2:828625/1‑149 (MQ=255)
                                 aCGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCTGATCAAGGCCAAt                                                                                                      >  2:39705/1‑149 (MQ=255)
                                 aCGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                      >  1:159637/1‑149 (MQ=255)
                                  cGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                     >  2:326827/1‑149 (MQ=255)
                                   gaaggaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                                    <  2:350955/149‑1 (MQ=255)
                                         gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                              <  1:76213/149‑1 (MQ=255)
                                            gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                           <  1:326827/149‑1 (MQ=255)
                                                        aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCCGGAACTCCTGACCGCCGAACt                                                                               >  2:521594/1‑149 (MQ=255)
                                                          cGCCAAACGGCAGCACCTGATTACCCGGTACAAAACCAGCAATTAAAATGGTGGGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                             <  2:640057/149‑1 (MQ=255)
                                                           gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                            >  1:643546/1‑149 (MQ=255)
                                                                   gcagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                                    <  1:580932/149‑1 (MQ=255)
                                                                    cagcaCCTGATTACCCGGCACACAAACAGCAATTTAAATGGTGAGTGGGATAAAACTCAGGCTCGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                   <  1:816381/149‑1 (MQ=255)
                                                                               tACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                   <  1:522872/106‑1 (MQ=255)
                                                                               tACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                   >  2:522872/1‑106 (MQ=255)
                                                                                 cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTg                                                      <  2:703569/149‑1 (MQ=255)
                                                                                  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCTGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGc                                                     <  1:39705/149‑1 (MQ=255)
                                                                                  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGACCCTGc                                                     >  1:260013/1‑149 (MQ=255)
                                                                                  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCAGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGACTGc                                                     >  2:369371/1‑149 (MQ=255)
                                                                                    ggCACACAAACGGCAAATAAAATGGTGAGTGGGATAAAAATCAGGATTGCCGATACCA‑‑CCGCAGTATGTCCCAGCAGCAACGACGGATCAAGGCCAATCACGAACTCCTGACCGCCGAACTTCGCCTTCAGACGCCCACGAGCCTGCtt                                                   >  2:828195/1‑149 (MQ=255)
                                                                                    ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                           >  1:1633/1‑125 (MQ=255)
                                                                                    ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                           <  2:1633/125‑1 (MQ=255)
                                                                                    ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                   >  2:371464/1‑149 (MQ=255)
                                                                                       acacaAACAGCAATTAAAATGGTGAGTAGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCGGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCAGACCGCCGAACTTCGCCTGTAAACGACTACGGGCCTGCTTaac                                                >  1:259482/1‑147 (MQ=255)
                                                                                          caAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                             <  1:730084/149‑1 (MQ=255)
                                                                                          caAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                             <  1:665169/149‑1 (MQ=255)
                                                                                          caAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                             <  1:664936/149‑1 (MQ=255)
                                                                                           aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATg                                            <  1:371464/149‑1 (MQ=255)
                                                                                           aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATg                                            <  2:643546/149‑1 (MQ=255)
                                                                                                 cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                   <  1:720982/120‑1 (MQ=255)
                                                                                                 cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                   >  2:720982/1‑120 (MQ=255)
                                                                                                    ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCCAACTTCGCCTGTAAACGACTACGAGCCTCCTTAGCGATAGGCGTTaca                                   >  1:211694/1‑147 (MQ=255)
                                                                                                     tAAAATGGTGGGTGGGATCGAAATCAGGCTTGCGGAT‑‑GGCCCCGCCGTATGGCCCAGCACCATCGGCTGAACGAGGCTGTTGATGAACTTCGTACCGGCGAACTGGGCGTGTAAACGACGACGAGCCTGCTTAGCGATGGGCGTTAAAc                                  <  2:450538/149‑1 (MQ=255)
                                                                                                        aaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                               <  1:746298/149‑1 (MQ=255)
                                                                                                               aGTGGGCTAAAAATCAGGCTTGCCGATA‑‑GCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                        <  1:890541/149‑1 (MQ=255)
                                                                                                                   ggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt                    <  1:459539/149‑1 (MQ=255)
                                                                                                                    gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt                   <  1:893114/149‑1 (MQ=255)
                                                                                                                    gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCACATCAAGGCCAATCAGGCACTCCTAACCGCCGAACTTCGCCTGTAAACGCCTACGACCCTGCTTAGCGATGGGCGTTAAACCATCCATCATCGGtt                   >  2:337069/1‑149 (MQ=255)
                                                                                                                     aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGatt                  >  1:325809/1‑149 (MQ=255)
                                                                                                                           aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTGTAATCAc            >  1:197695/1‑149 (MQ=255)
                                                                                                                              aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg         >  2:43810/1‑149 (MQ=255)
                                                                                                                                     ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt  <  2:906673/149‑1 (MQ=255)
                                                                                                                                              ||                                                                                                                                            
GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATT  >  NC_000913/2173220‑2173502

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: