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.94.4% 141.3 / 0.4 36pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (16/18);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01
*NC_0009132,173,3620C.94.3% 132.9 / 2.1 36pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (16/17);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCG  >  NC_000913/2173221‑2173495
                                                                                                                                            ||                                                                                                                                     
aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                      <  2:375831/143‑2 (MQ=255)
     cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGCTCGCCAAACGGCAGCACCTGGTTACCCGGCACACAAACAGCAATTCAAATGGTGAGTGGGATAAAAATCAGGCCTGCCGATaccaccgccgta                                                                                                                                 >  1:358871/1‑137 (MQ=255)
     cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCGGGCACACAAACAGCAATGAAAATGGTGCGTGGGATAAAAATCAGGCGTGCCGATaccaccgccgtg                                                                                                                                 >  1:190292/1‑137 (MQ=255)
     cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgta                                                                                                                                 >  1:357308/1‑137 (MQ=255)
             aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                       >  2:370534/1‑143 (MQ=255)
             aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGGATGTCCcaa                                                                                                                       >  2:39680/1‑142 (MQ=255)
             aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTCGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                       <  1:671824/143‑1 (MQ=255)
               ggCGACCGCCATCGCCACGAAGATGCCAATGGTGGCAAGATCGCCCAGCGGCGTCACGTGATGACCCTGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                     <  2:320332/143‑1 (MQ=255)
                          tggCCAGGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                          <  2:624450/141‑1 (MQ=255)
                              cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                      >  1:524728/1‑143 (MQ=255)
                                         cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                           <  1:681268/143‑1 (MQ=255)
                                                 gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGACACAGGCCAACCAGGAACTCCTg                                                                                   >  2:147133/1‑143 (MQ=255)
                                                   aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                 >  1:474392/1‑143 (MQ=255)
                                                          ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                          <  1:716772/143‑1 (MQ=255)
                                                          ccAAACGGCAGCACCTGATTACCCGGCAAACAAAGAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAACAAGGCCAATCAGGAAATCATGACCGCCGaa                                                                          >  2:222721/1‑143 (MQ=255)
                                                                  cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGAGATAAAAACAAGACTGTGCGTTGC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGTATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                  <  1:268842/143‑1 (MQ=255)
                                                                    gcaCCTAATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                <  1:616824/143‑1 (MQ=255)
                                                                       ccTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAGCAGGCTTGCCGATACCA‑‑CCGCCGGATGTCCCAGCAGCAACGCGGGAGCAAGGCCAATCAGGAACTCCGGACCGCCGCACTTCGCCTGTaaa                                                             >  1:589033/1‑143 (MQ=255)
                                                                        cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCTCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                            <  1:567930/143‑1 (MQ=255)
                                                                             tACCCGGCACACAAACAGGAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                       <  1:638826/143‑1 (MQ=255)
                                                                               cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                     >  1:242650/1‑143 (MQ=255)
                                                                                ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACCTCGCCTGTACACCACTACGa                                                    >  2:497006/1‑143 (MQ=255)
                                                                                     acacaAACAGCATTTAAAATGGTGAGTGTAATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCTGCAGCCAAGAAGAAATGATGCCATTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTg                                               <  1:568234/143‑1 (MQ=255)
                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTg                                               <  1:263478/143‑1 (MQ=255)
                                                                                       acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTACCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACAACGAgaccgcc                                             >  2:100463/1‑137 (MQ=255)
                                                                                         aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCGTGATTa                                           >  1:19954/1‑143 (MQ=255)
                                                                                         aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                           <  1:518103/143‑1 (MQ=255)
                                                                                               cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATg                                     >  1:695955/1‑143 (MQ=255)
                                                                                                   tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCCTGCTTCGCGATGGGCg                                 >  1:583974/1‑143 (MQ=255)
                                                                                                       aTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                             <  1:447729/143‑1 (MQ=255)
                                                                                                          gTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGACTGCTGAGCGATGGGCGTTAAAcc                          >  1:707086/1‑143 (MQ=255)
                                                                                                                   agaaagaTCAGGCTTGGCGATACCA‑‑CCGCCGTATGTCCCAGCACCAACTCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                 <  2:513012/141‑1 (MQ=255)
                                                                                                                   aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCAAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                 <  1:577153/143‑1 (MQ=255)
                                                                                                                     aaaaaCCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg               <  1:344399/143‑1 (MQ=255)
                                                                                                                          tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGCTTa          >  1:197053/1‑143 (MQ=255)
                                                                                                                           cacgCTCGCCGACACCGCCCGCCTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGCAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa         <  2:274023/140‑1 (MQ=255)
                                                                                                                           cAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGTCCGCCGAACTTCGCCTGTAAACGCCTACGAGTCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa         <  2:221350/143‑1 (MQ=255)
                                                                                                                              gCTTGCCGATACCA‑‑CCGACGTATGTCCCCGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa      <  2:81596/143‑1 (MQ=255)
                                                                                                                                ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc    >  1:368531/1‑143 (MQ=255)
                                                                                                                                  gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg  <  1:39680/143‑1 (MQ=255)
                                                                                                                                            ||                                                                                                                                     
ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCG  >  NC_000913/2173221‑2173495

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: