Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,361 +C pseudogene (918/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,3601.C94.0% 179.2 / 3.4 50pseudogene (919/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (22/25);  total (23/27)
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.90e-01

ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCA  >  NC_000913/2173221‑2173491
                                                                                                                                            |                                                                                                                                   
tttcccACGATCCCCGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATacca                                                                                                                                      <  2:787609/136‑1 (MQ=255)
aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATacca                                                                                                                                      >  2:898070/1‑140 (MQ=255)
 tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCAC                                                                                                                                    <  1:898070‑M1/140‑1 (MQ=255)
     cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                                                                                                >  1:887498/1‑137 (MQ=255)
           gCACGGCGTCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgt                                                                                                                           <  1:19681‑M1/140‑10 (MQ=255)
                    ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagc                                                                                                                  <  2:915645‑M1/140‑19 (MQ=255)
                    ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagc                                                                                                                  >  2:678091‑M1/1‑122 (MQ=255)
                    ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagc                                                                                                                  >  1:186243‑M1/1‑122 (MQ=255)
                    ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagc                                                                                                                  <  2:186243‑M1/140‑19 (MQ=255)
                    ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagc                                                                                                                  >  1:946569‑M1/1‑122 (MQ=255)
                      gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaa                                                                                                                >  2:115492‑M1/1‑120 (MQ=255)
                         aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgc                                                                                                             >  1:101279‑M1/1‑117 (MQ=255)
                         aTCGCCAAGATGAACCCGATGGTGGAAAAATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACGGCAATTAAAATGGTGAGCGGGATAAAAACCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCCACGc                                                                                                          <  2:101279/140‑1 (MQ=255)
                          tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgcc                                                                                                            <  1:304388‑M1/140‑25 (MQ=255)
                              cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggat                                                                                                        <  2:152494‑M1/140‑29 (MQ=255)
                                   agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaa                                                                                                     <  2:791986‑M1/138‑32 (MQ=255)
                                   agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaa                                                                                                     >  1:791986‑M1/1‑107 (MQ=255)
                                       gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaa                                                                                               <  1:678091‑M1/140‑38 (MQ=255)
                                       gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaa                                                                                               <  2:887498‑M1/140‑38 (MQ=255)
                                             ggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACgccgtatgtcccagcagcaacgccggatcaaggccaatc                                                                                             >  1:587780‑M1/1‑96 (MQ=255)
                                             ggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACgccgtatgtcccagcagcaacgccggatcaaggccaatc                                                                                             <  2:587780‑M1/135‑40 (MQ=255)
                                              gtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaa                                                                                        <  1:670042‑M1/140‑45 (MQ=255)
                                                   aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcct                                                                                   <  1:115492‑M1/140‑50 (MQ=255)
                                                      aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgac                                                                                >  1:917611‑M1/1‑88 (MQ=255)
                                                        cGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccg                                                                              >  2:496078‑M1/1‑86 (MQ=255)
                                                         gcaaaacggCAGCACCTGATTCCCCGGCACACAAACAGCAATAAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgc                                                                             <  2:146751‑M1/137‑56 (MQ=255)
                                                               cGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaact                                                                       >  2:402552‑M1/1‑79 (MQ=255)
                                                                      aCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctg                                                                <  1:271681‑M1/140‑69 (MQ=255)
                                                                           aTTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaac                                                           <  2:559608‑M1/140‑74 (MQ=255)
                                                                           aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaac                                                           <  2:917611‑M1/140‑74 (MQ=255)
                                                                                acagcacACAAACAGCAATAAAAATGGTGAGTGGAATAAAAATCAGGCTTGCCGATCCCACCgccgtatgtcccagcaacaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgacta                                                      <  2:120514‑M1/137‑79 (MQ=255)
                                                                                  ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacg                                                    >  1:1007450‑M1/1‑60 (MQ=255)
                                                                                    cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagc                                                  >  1:948444‑M1/1‑58 (MQ=255)
                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagc                                                 <  1:109559‑M1/140‑84 (MQ=255)
                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagc                                                 <  1:620142‑M1/140‑84 (MQ=255)
                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagc                                                 >  2:441692‑M1/1‑57 (MQ=255)
                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagc                                                 >  2:277453‑M1/1‑57 (MQ=255)
                                                                                          aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgct                                            >  1:392268‑M1/1‑52 (MQ=255)
                                                                                              gCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagc                                        >  1:594663‑M1/1‑48 (MQ=255)
                                                                                               cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcaacaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttctgt                                       >  1:120514‑M1/1‑47 (MQ=255)
                                                                                                         ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgtta                             >  2:37986‑M1/1‑37 (MQ=255)
                                                                                                          gTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaa                            >  1:926395‑M1/1‑36 (MQ=255)
                                                                                                          gTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccgaatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaa                            >  1:456599‑M1/1‑36 (MQ=255)
                                                                                                           ggAGTGGAAAAAAAATCAGGCTTGCCGATCCCCCCGCCGTATGTCCCAGCAGCAACGCCGAATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                        <  2:456599/139‑1 (MQ=255)
                                                                                                             aGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaaacc                         >  2:371590‑M1/1‑33 (MQ=255)
                                                                                                             aGTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaaacc                         <  1:37986‑M1/140‑108 (MQ=255)
                                                                                                              gTGGGATAAAAATCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaaacca                        <  1:402552‑M1/140‑109 (MQ=255)
                                                                                                                ggaataaAAACCGGGCTTGCCGATCCCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaaaccatc                      <  1:782814‑M1/137‑111 (MQ=255)
                                                                                                                     aaaaaTCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaaaccatccatga                 <  1:371590‑M1/140‑116 (MQ=255)
                                                                                                                         accagGCTTGCCGATCCCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaaaccatccatgatggg             <  2:360806‑M1/138‑120 (MQ=255)
                                                                                                                          tCAGGCTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaaaccatccatgatgggt            <  1:277453‑M1/140‑121 (MQ=255)
                                                                                                                               cTTGCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaaaccatccatgatgggtttaat       <  2:1007450‑M1/140‑126 (MQ=255)
                                                                                                                                  gCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaaaccatccatgatgggtttaatcac    >  1:55003‑M1/1‑12 (MQ=255)
                                                                                                                                  gCCGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaaaccatccatgatgggtttaatcac    >  2:905842‑M1/1‑12 (MQ=255)
                                                                                                                                    ggaTACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctacttagcgatgggcgttaaaccatccatgatgggtttaatcaccc  <  2:822270‑M1/139‑131 (MQ=255)
                                                                                                                                    cGATACCACCgccgtatgtcccagcagcaacgccggatcaaggccaatcaggaactcctgaccgccgaacttcgcctgtaaacgactacgagcctgcttagcgatgggcgttaaaccatccatgatgggtttaatcaccc  <  1:496078‑M1/140‑131 (MQ=255)
                                                                                                                                            |                                                                                                                                   
ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCA  >  NC_000913/2173221‑2173491

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