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.97.5% 169.0 / ‑2.2 40pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (13/26);  total (13/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.01e-01
*NC_0009132,173,3620C.97.5% 167.8 / ‑3.3 40pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (13/26);  total (13/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 = 1.00e+00

ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173221‑2173496
                                                                                                                                            ||                                                                                                                                      
aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                       >  1:228364/1‑142 (MQ=255)
  ttCCGCGCTGCACGGCGACCGCCATCGGCACGAAGAAGCCGGTGGTGACAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                                                                                                     <  2:1169545/143‑4 (MQ=255)
                   actgCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAACATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCaa                                                                                                                  <  1:611842/140‑1 (MQ=255)
                    ccGTCATAGCCACGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                                 <  1:630931/143‑1 (MQ=255)
                      gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                               <  1:1173112/143‑1 (MQ=255)
                       ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                              >  1:651799/1‑143 (MQ=255)
                                cGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                     <  2:655312/143‑1 (MQ=255)
                                cGAAGAAGCCGATGGAGGCAAGAACGCCAAACGGCAGCACCTGATTACCCGGCACACAGACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                     <  1:842354/143‑1 (MQ=255)
                                    gaagCCGGTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                 <  2:30290/143‑1 (MQ=255)
                                     aagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                <  1:143496/143‑1 (MQ=255)
                                           atgtcggCAACATCGCCAAACGGCAGCACTACAGTACCCGGCACAAACACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCTCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                          <  2:126665/138‑1 (MQ=255)
                                           aTGGTGGCAAGATCGCCAAACCGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                          <  1:1058707/143‑1 (MQ=255)
                                             gctggAAAGATCTCCAAATGGCAGGACCTGAGTTCCCAGCAAGCAATCAGCAATTAAACTGGTGAGTGGGATTAAAAGCAGGCTGGCCGATACCA‑‑CCGCCTTATGTCCCTGCAGCAACGCCTGATCAAGGCCAATCAGGAACt                                                                                        <  2:821211/141‑1 (MQ=255)
                                               tggCCAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc                                                                                      <  1:59441/143‑1 (MQ=255)
                                               tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCAACGCCGGATCAAGGCCAATCAGGGACTcc                                                                                      >  2:1153193/1‑143 (MQ=255)
                                                     gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑TCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGACAATCAAGAACTCCTgcacg                                                                                >  2:916324/1‑139 (MQ=255)
                                                      aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAACATGGTGAGTGGGAGAAAAATCAGGGTTGCCGATACCA‑‑ACGCCGTATGACCCAGCCGAAACGCCGGATCACTGGCCACCAGAAACTCCTCAccga                                                                               >  2:357033/1‑142 (MQ=255)
                                                           cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATGAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCAc                                                                          >  2:944877/1‑143 (MQ=255)
                                                              aCGGCAGCACCCGATATCCCAGCACGCAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCTGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                       <  1:620017/143‑1 (MQ=255)
                                                                  cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                   <  2:592765/143‑1 (MQ=255)
                                                                    gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGTAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                 <  1:646695/143‑1 (MQ=255)
                                                                       ccTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                              <  1:37743/143‑1 (MQ=255)
                                                                         tGATTACCCGGCACACAAACAGCAATTAAAATGGTAAGAGGAATCAAAATCAGGCCTGCCGATACCA‑‑CCGCCGCATGTCCCACCAGCAACGCCGGCTCAAGGCCACTCACGAACACCCGACCCCCGAACGTCGCCCGTAAACg                                                            >  2:1209392/1‑143 (MQ=255)
                                                                           atgaCCCGGCACACAAGCAGCAATTAAGATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc                                                          <  2:228364/140‑1 (MQ=255)
                                                                                  ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                   <  1:925445/143‑1 (MQ=255)
                                                                                  ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACGACGAGc                                                   <  1:364229/143‑1 (MQ=255)
                                                                                   gCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgac                                                  >  1:17419/1‑141 (MQ=255)
                                                                                    cacacaAACATAAATTAATATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                 <  2:651799/143‑1 (MQ=255)
                                                                                       acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                              >  1:360642/1‑143 (MQ=255)
                                                                                            ccgCAATGAAAAGGGTGAGGCGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                         <  1:841554/141‑1 (MQ=255)
                                                                                             aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                        <  1:751466/143‑1 (MQ=255)
                                                                                               cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCGAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAtc                                      >  1:752665/1‑142 (MQ=255)
                                                                                                 aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATggg                                    >  2:100541/1‑143 (MQ=255)
                                                                                                    aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                 <  2:894439/143‑1 (MQ=255)
                                                                                                     aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCCGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCgct                                >  1:257175/1‑141 (MQ=255)
                                                                                                        tggtggGTGGGATAAGAATCCGGCTTCACGATA‑‑ACCCGGCCGCATGCCCCAGCAGCCACGACGGATAAAGGCCAATCAGGACCTCATGCCCGCCGAACTTCGCCTGTAAACGACTACGCGCCTGCTTAGCGCTGGGCGTTaaa                             <  2:854239/143‑1 (MQ=255)
                                                                                                             aGTGGGATAAGCATCAGTCTTGCCGATACCA‑‑CCGCGGTAGGTCCCAGCAGCAACGCCGGATCAGGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                        <  2:469565/143‑1 (MQ=255)
                                                                                                                   aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                  >  1:1136459/1‑143 (MQ=255)
                                                                                                                    tAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                 <  1:690293/143‑1 (MQ=255)
                                                                                                                     aaaaaTCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                <  1:397659/143‑1 (MQ=255)
                                                                                                                           cAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGAGGCGTCTaa          >  1:751441/1‑143 (MQ=255)
                                                                                                                                   ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt  <  1:21892/143‑1 (MQ=255)
                                                                                                                                            ||                                                                                                                                      
ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173221‑2173496

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: