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.93.1% 111.5 / 1.7 29pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (13/14);  total (14/15)
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.93.1% 111.2 / 1.2 29pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (13/14);  total (14/15)
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

CACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACC  >  NC_000913/2173226‑2173493
                                                                                                                                       ||                                                                                                                                   
cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTACAATGGTGAGTGGGCTCAAAATCAGGCGTGCCGCTacccccgccgcg                                                                                                                               >  2:579669/1‑137 (MQ=255)
         cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagc                                                                                                                    >  1:216112/1‑143 (MQ=255)
         cGGCGACCGCCAGCGCCCCGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagc                                                                                                                    <  2:9336/143‑1 (MQ=255)
                 gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGCGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGGCCCAGCAGCAACGc                                                                                                            >  2:439471/1‑143 (MQ=255)
                       gCCACGAAGAGGCCGAGGGGAGCAAGATCGCCAAACGGCAGCACCTGAGTAGCAGACACACAAGCAGCAATTAAAATGGTGAGTGGGATAAAGATGAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                      <  2:54466/143‑1 (MQ=255)
                          aCGGAGAAGCCGATGGTGGCTAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                   <  1:670236/143‑1 (MQ=255)
                           cGAAGAAGCCGATGGTGGCAGGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                  <  2:534351/143‑1 (MQ=255)
                             atgaagCCGTACGTGGCAAGATCGCCAAACGGCCGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                <  2:406860/141‑1 (MQ=255)
                                  gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACAGCAATTCAACTGGTGCGTGGGATACAAATCAGGCTTCCCGATACCA‑‑CCGCCGTATGTCCACGCCGCACCGCCGGCTCAACGccgctcc                                                                                           >  1:199538/1‑138 (MQ=255)
                                                gATTGCCACACGTGAGCACCGGCGTGCCCGGCACACAAACAGCAGTTATAGTGGTGAGTGGGTTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                             <  2:191417/143‑1 (MQ=255)
                                                 aTCGCCAAACGGCAGCACCTAATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGAATAAAAAAAAAACTAGCAGAAACAA‑‑CCGACATATATCACAGCACCAACGCCGGATCAAAACCAATCAGGAACACCTGAacgc                                                                            >  2:709411/1‑143 (MQ=255)
                                                     ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                        <  1:53601/143‑1 (MQ=255)
                                                      cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTCGCCGATACCA‑‑CCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCGCATCAGGAACTCCTGACCGCCGAAc                                                                       >  2:533622/1‑143 (MQ=255)
                                                       aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                      >  1:604697/1‑143 (MQ=255)
                                                         aCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTCGCCGATACCACCC‑‑CCGTATGTCCCCGCCGCAACGCCTGATCAACGCCAATCGCGCACTCCTCACCGCCGAACTTc                                                                    >  2:216962/1‑143 (MQ=255)
                                                          cGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                   <  2:604697/143‑1 (MQ=255)
                                                            gcagcaCCTGACTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                 <  1:69254/143‑1 (MQ=255)
                                                             cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                <  1:674455/143‑1 (MQ=255)
                                                                   cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaacc                                                          >  1:296338/1‑141 (MQ=255)
                                                                      aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc                                                       <  1:414102/143‑1 (MQ=255)
                                                                                  acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                           >  2:356309/1‑143 (MQ=255)
                                                                                             aaagaaTGGGGAGGGGGAGAAAAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                <  1:62018/139‑1 (MQ=255)
                                                                                                aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCCAGGCCAATCAGGAACTCCTGACCCCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                             >  2:611836/1‑143 (MQ=255)
                                                                                                   gggTGAGTGGGATAAAAATCGGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                          <  1:451221/142‑1 (MQ=255)
                                                                                                            ggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTCGCGATGGGCGCTACAACATCCatg                 >  1:711728/1‑143 (MQ=255)
                                                                                                             gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAGCGCCGGATCAAGGCCAATAAGGAACTACTAACCGCCGAACTTCGCATGTAAAAGACTACGAGCCTGCATAGCGATGGGCGTTTAACCATCCatga                >  2:670687/1‑143 (MQ=255)
                                                                                                             gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTCTAAACGACTACGAGCCTGCTTAGCGATGCGCGTTAAACCATCCatga                >  1:707901/1‑143 (MQ=255)
                                                                                                                    aTCAGGCTTGCCGATAC‑‑CCCCGCCGGACGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTGCGCCTGTACACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt         <  1:370564/143‑1 (MQ=255)
                                                                                                                     tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTa        <  1:80182/143‑1 (MQ=255)
                                                                                                                           ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATCGGCGTTAAACCATCCATGATAGGTTTAATCAcc  >  2:208386/1‑143 (MQ=255)
                                                                                                                                       ||                                                                                                                                   
CACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACC  >  NC_000913/2173226‑2173493

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: