Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp pseudogene (915‑916/1358 nt) gatC ← galactitol‑specific PTS enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.92.2% 203.3 / 7.1 51pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base . (17/30);  total (18/33)
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.94.0% 201.0 / 4.7 51pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (17/30);  total (18/33)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173222‑2173497
                                                                                                                                           ||                                                                                                                                       
tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATCAAAATGGTGAGCGGGATAAAAATCAGGCTTGCCGATcccaccgc                                                                                                                                       >  2:758788/1‑141 (MQ=255)
 caccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGACGGTGGCAAGAGCGCCAAACGGCAGCACCTGAGTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                                                                                                      <  1:1068711/141‑4 (MQ=255)
       cgtGCCCGGCGACCGCCATTGCCACGAGGAAGCCCGTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCCGCACACAAACAGCAATTAAAATGGTGAGTGGGACAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                              <  1:706525/141‑1 (MQ=255)
        ctGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGTCAAGGTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTc                                                                                                                             <  2:9333/142‑1 (MQ=255)
           cACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                          >  1:549125/1‑143 (MQ=255)
            aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACCCAAACAGCAGTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                         >  1:832210/1‑143 (MQ=255)
            aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGGCCcag                                                                                                                         >  2:153885/1‑143 (MQ=255)
              ggCGCCCGCCATCGCCACGAAGCAGCCGATGGTGTCAAGATCGCCAAACGGCAGCACCTCATTACCCGGCACACAAACAGCAATCAAAATTGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                       <  1:696077/143‑1 (MQ=255)
              ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAGAGTGGCGACCGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTCTGTCCcagca                                                                                                                       >  1:829066/1‑143 (MQ=255)
              ggCGACCGCCATCGCCACGAAGAAGCCGATGGTCGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACATCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                       <  1:436263/143‑1 (MQ=255)
                cggccGCCAGCGCCACGAAGGAGCCGATGGTGGCAAGAGCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                     <  2:816721/140‑1 (MQ=255)
                   ccGCCAGCGCCACGAAGGAGCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                                  <  1:846701/143‑1 (MQ=255)
                    cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                                                                                                 >  2:832458/1‑143 (MQ=255)
                      ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                               <  2:734109/143‑1 (MQ=255)
                              aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                       <  2:836948/143‑1 (MQ=255)
                               cGAAAAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                      <  1:25431/143‑1 (MQ=255)
                                  agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCAAACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCCTGCCGAAA‑‑ACACCGACGTAGGTCCCAGCCGCAAGGCCGGATCACGGcaa                                                                                                   >  2:744418/1‑141 (MQ=255)
                                    aagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                 <  2:1061676/143‑1 (MQ=255)
                                      gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGACACCA‑‑CCGCCGTATGGCACAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                               >  2:10011/1‑143 (MQ=255)
                                        cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                             <  1:583611/143‑1 (MQ=255)
                                        cGACGGTGGCAAGTTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATCAgg                                                                                             <  1:812492/143‑1 (MQ=255)
                                         gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTACGTCCCAGCAGCAACGCCGGAGCGACGCCACTCAGGa                                                                                            >  2:494743/1‑143 (MQ=255)
                                          atagtgcCAAGATCGCCAAGCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                           <  2:1079627/140‑1 (MQ=255)
                                                  aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                   <  1:494200/143‑1 (MQ=255)
                                                             gcGGCAGCACCTGATTACCCGGCACACAAACAGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                        <  1:875909/142‑1 (MQ=255)
                                                             aCGGCAGCACCTGAATACCCGGCACACAAACAGCAAAGAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                        <  1:323990/143‑1 (MQ=255)
                                                               ggCAGCACCCGATCACCCGGCACACAAACCGCAATTAACATGGTGAGTGCGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      <  1:844488/143‑1 (MQ=255)
                                                                 cagcaCCGGATGACCAGGCCCACAAACAGCAAGTAAACGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCC‑‑CAGTATGTCCCAGAAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                    <  1:497728/143‑1 (MQ=255)
                                                                         gATTACCCGGCACACAAACAGCCATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                            <  2:658192/143‑1 (MQ=255)
                                                                         gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                            <  1:317666/143‑1 (MQ=255)
                                                                               ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGAAAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                      >  1:737790/1‑143 (MQ=255)
                                                                                  gCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTCCggccc                                                   >  2:646576/1‑139 (MQ=255)
                                                                                      acaAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGGAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCAGTTCGCCGGTACACGCCTACGAGCCGGCt                                               >  1:331058/1‑143 (MQ=255)
                                                                                        aaaCAGCAATTAAACTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                             <  1:106822/143‑1 (MQ=255)
                                                                                        aaaCAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                             <  1:976775/143‑1 (MQ=255)
                                                                                         aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                            <  1:1002470/143‑1 (MQ=255)
                                                                                            cgCAACTAAAAGGGTGGGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCCGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                         <  2:439998/142‑1 (MQ=255)
                                                                                            aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGGGTCTTCGTAGCGa                                         >  1:703870/1‑143 (MQ=255)
                                                                                              cgactAACATGGTGAGTGGGATAAAAATCAGGCGTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATg                                       <  2:109430/139‑1 (MQ=255)
                                                                                              aaaTTAAAGTGGCGAGTGGGATAAAAATCAGGGGGGCACAGACCA‑‑CCGGCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTGCGCCTGTAAACGACTACGAGCCTGCTTAGCGATg                                       <  2:499553/142‑1 (MQ=255)
                                                                                                    aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAGCCGGCTACGAGCCTGCTTGGCGCTGTGCGtt                                 >  1:508002/1‑143 (MQ=255)
                                                                                                        ggTGAGTGGGATAAACATCCGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCGGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                             <  2:991688/143‑1 (MQ=255)
                                                                                                             gTGGCATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGCCCCAGCAGCAACGCCGGATCAGGGGTCCCCATCAACTACAGACACCCGTGCTTCTCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcc                        <  1:834854/143‑1 (MQ=255)
                                                                                                                 gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                          <  2:215746/41‑1 (MQ=255)
                                                                                                                 gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                          >  1:215746/1‑41 (MQ=255)
                                                                                                                  aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCctgat                   >  1:122257/1‑143 (MQ=255)
                                                                                                                  aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCACCGAACTTCGCCTGTAACCGACTACGCGCCTGCTTAGCGCTGGGCGTTAAACCCTCCatgat                   >  2:11270/1‑143 (MQ=255)
                                                                                                                   tAAAAACCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                  <  1:316219/143‑1 (MQ=255)
                                                                                                                       aaccggccTTGCCGACACCA‑‑CCGCCGTATGGCCCAGCAGCAACGCCGGAGCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGGAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt              <  2:401836/136‑1 (MQ=255)
                                                                                                                         gcAGGCTTGCCGAGACCA‑‑CCGCCGTCTGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTa            <  2:335433/142‑1 (MQ=255)
                                                                                                                             gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa        <  1:832458/143‑1 (MQ=255)
                                                                                                                                   cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg  >  1:608556/1‑143 (MQ=255)
                                                                                                                                   cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACATCGACTATAAACGAATACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCCTg  >  1:1060386/1‑143 (MQ=255)
                                                                                                                                           ||                                                                                                                                       
TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173222‑2173497

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: