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.96.0% 210.0 / 0.9 50pseudogene (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 (0/2);  new base . (23/25);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
*NC_0009132,173,3620C.96.0% 207.9 / 0.1 50pseudogene (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 . (23/25);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA  >  NC_000913/2173221‑2173488
                                                                                                                                                ||                                                                                                                              
aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                               >  1:1077531/1‑142 (MQ=255)
      aCGATGCCCGGCGCCCGCCATCCCCACGAAGACCCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                         <  1:81663/143‑8 (MQ=255)
            cACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGACAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                 <  2:500186/143‑1 (MQ=255)
             aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                >  1:197764/1‑143 (MQ=255)
             aCGGCGACCGCCAGCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCGATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                <  1:875501/143‑1 (MQ=255)
                    ccGCCATCGCAACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                         <  1:781196/143‑1 (MQ=255)
                       ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAAACAAACCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                      >  1:628005/1‑143 (MQ=255)
                         aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCGGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                    >  1:474484/1‑143 (MQ=255)
                          tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                   <  2:165198/143‑1 (MQ=255)
                           cGCCAAGAAGAAGCCGCGGGTGGCAAGAGGGCCAAACCGCAGCACCGGATTACCCGGCAAACAAACAGCGATTAAAAT‑‑GGTGAGTGGGATAAGAATCAGGCTTGGCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                  <  1:217345/143‑1 (MQ=255)
                              cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                               >  1:1036709/1‑143 (MQ=255)
                              cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                               >  1:513359/1‑143 (MQ=255)
                               aCGAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                              <  1:232257/143‑1 (MQ=255)
                                cgacgaAGCCGATGGTGGCAAGATCGCCAAACGGTACCACCTTTTTACCCGGCACACAAACAG‑‑ATCAACTTTAGGTGACTTCGCTAAAAAACATTTCGTTGAGTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                             <  1:803918/143‑1 (MQ=255)
                                     ccgCCGATGGTGGCAAGCTCGCCAAACGGCAGCACGCGACTACTCGGCCCACAAACAGCACTCACAAA‑‑GCCGAGTGTTGCTCTGGACATCCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCTGATCAAGGCCAAt                                                                                        <  2:231348/141‑1 (MQ=255)
                                      agCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                       >  1:437731/1‑143 (MQ=255)
                                      agCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                       <  1:483384/143‑1 (MQ=255)
                                       gcagATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTAGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                      >  2:1022375/4‑143 (MQ=255)
                                         cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGTCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                    >  1:227897/1‑143 (MQ=255)
                                          gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTACCAGCAGCAACGCCGGACCAAGGCCGATCAGGa                                                                                   >  2:462716/1‑143 (MQ=255)
                                           aTGGTGTCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                  <  1:678412/143‑1 (MQ=255)
                                              caggCAAGAACCCGAAACGGCAGCACCCGAGTACCCGGCACACAAACAGCAAGTAAAAT‑‑GGTGAGTAGGACAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                               <  1:121391/141‑1 (MQ=255)
                                               tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc                                                                              >  1:108372/1‑143 (MQ=255)
                                                 gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAACTCCTg                                                                            >  2:126685/1‑143 (MQ=255)
                                                      aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                       >  1:1031436/1‑143 (MQ=255)
                                                      aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGCGAGTGGGATAACAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACACCGGACCAAGGCCAATCAGGACCCCCTGAcccg                                                                       >  2:299370/1‑141 (MQ=255)
                                                       tCGCCAAACGGCAGCACCTGATTACCCAGCACACAAACAGCAACTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                      <  1:404088/143‑1 (MQ=255)
                                                         gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGa                                                                    >  1:521497/1‑143 (MQ=255)
                                                            aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGCACt                                                                 >  2:218561/1‑143 (MQ=255)
                                                                ggCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                             <  2:252335/143‑1 (MQ=255)
                                                                ggCAGCACCTCATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                             <  1:483269/143‑1 (MQ=255)
                                                                 gcagcaCCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGGGGGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                            <  2:619903/143‑1 (MQ=255)
                                                                     caCCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGCCTTGCCGAAACCA‑‑CCGCCGGATGTCCCAGCAGCAACGCCGGATCAAGGCCAGTCAGGCACTCCGGACGGCCGCGCTTCGccccca                                                        >  2:1098757/1‑139 (MQ=255)
                                                                        cTGATGTCCCGGCACATCAACAGCAATGGAAGT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACTGCCGAACTTCGCCTGTAAAc                                                     <  2:742593/143‑1 (MQ=255)
                                                                                ccGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                             <  1:778335/143‑1 (MQ=255)
                                                                                         aaaCAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                    <  1:1054222/143‑1 (MQ=255)
                                                                                           aCAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCGGTAAACGACTACGCGCGTGCTTAGc                                  >  1:932381/1‑143 (MQ=255)
                                                                                           aCAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGAGCCTGCTTAGc                                  >  2:651961/1‑143 (MQ=255)
                                                                                             aGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTCCGCCTGTACACGACTACGCCCCTGCTTAGCGa                                >  2:664980/1‑143 (MQ=255)
                                                                                             aGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTAAGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACCTCACCTGTAAACGACGCCGAGCCTCCTACGCGa                                >  2:1005202/1‑143 (MQ=255)
                                                                                               cAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGCCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTACGTCGCGATg                              >  2:208566/1‑143 (MQ=255)
                                                                                                   tAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCg                          <  2:203143/143‑1 (MQ=255)
                                                                                                     aaaT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATTCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACCTCGCCTGTAAACGACTACGAGCCTGCTCAGCGATGGGCGtt                        >  2:488033/1‑143 (MQ=255)
                                                                                                       aT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCATGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                      <  1:777019/143‑1 (MQ=255)
                                                                                                       aT‑‑GGTGAGGGGGATAAAAATCGGGCTTGCCGCGAC‑‑CCCCGCCGTGCGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                      <  1:851089/143‑1 (MQ=255)
                                                                                                        t‑‑GGTCAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                     <  1:629830/143‑1 (MQ=255)
                                                                                                               aGTGTGACAACAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCCATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                <  1:98842/143‑1 (MQ=255)
                                                                                                                 ggggAGAAACATGACGCTGGCCGATGCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTGCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcca              <  1:384638/142‑1 (MQ=255)
                                                                                                                     gATAAAAATCAGGCGTGCCGACACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga           <  2:692423/143‑1 (MQ=255)
                                                                                                                        aaacaTCAGCCTTGCCGATACCA‑‑CCGCCGTATGTCCCACCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg        <  1:23305/143‑1 (MQ=255)
                                                                                                                             tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTCGACGACTACGAGCCTGGTTAGCGATGGGCGGTCAGCCATCCATGCTGGGTTTa   >  1:435763/1‑143 (MQ=255)
                                                                                                                              cAGGCTTGCCGCTACCA‑‑CCTCCGTATGTCCCAGCAGCAACGCCTGATCAAGGCCAATCAGGAACTCCTGCCCGTCGTAGTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa  <  1:153409/143‑1 (MQ=255)
                                                                                                                                                ||                                                                                                                              
ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA  >  NC_000913/2173221‑2173488

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