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.95.7% 190.0 / 0.0 46pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (20/24);  total (21/25)
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.98e-01
*NC_0009132,173,3620C.97.8% 188.1 / ‑0.8 46pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (20/24);  minor base A (1/0);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACC  >  NC_000913/2173224‑2173493
                                                                                                                                         ||                                                                                                                                   
tCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGTTGGTGGCAAGCTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGCTAAAAAGCAGGCTTGCCGATaccaccgccg                                                                                                                                 >  2:976692/1‑139 (MQ=255)
     gATGCACGGCGACCGCCATCGCCACGACGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                          <  2:560009/143‑1 (MQ=255)
       tGCACGGCCACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACTAGCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcc                                                                                                                        <  1:929086/143‑1 (MQ=255)
         cACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                      <  1:162606/143‑1 (MQ=255)
            ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAAGTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCGGAGACCA‑‑CCGCAGTATGTcccgcca                                                                                                                   >  2:493567/1‑139 (MQ=255)
            cgccaCCGCCATCGCCACGAATAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATTGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                   <  1:452189/139‑1 (MQ=255)
                     cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                                                                                          >  2:393070/1‑143 (MQ=255)
                       tCGCCACGAAGAAGGCCATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                        <  1:215654/143‑1 (MQ=255)
                       tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTCAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                        >  2:760182/1‑143 (MQ=255)
                                    gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATACAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCa                                                                                           >  2:597319/1‑143 (MQ=255)
                                    gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                           <  1:629192/143‑1 (MQ=255)
                                     ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                          <  1:44994/143‑1 (MQ=255)
                                      cGATGGGGGCAAGATCGCCAAACGGCAGCGCCAGATGACCCGGTCCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                         <  2:312950/143‑1 (MQ=255)
                                       gATGGTGGCAAGATCACCAAAAAACAGAACCTGATTACCCGGCACACAAACAGCAATTAAAATAATGAGAAAGAAAAAAAACAGGATTGCCAATACAA‑‑CCGCAATATATCCCAGCAACAACGACAGACCAAAGCCAAACaaaa                                                                                        >  2:160864/1‑140 (MQ=255)
                                        aTGGTGGCAAGAGCGCGAAGCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAGCGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                       <  1:522253/143‑1 (MQ=255)
                                           gtggCTAGATCGCCAGACGGCAGCACCTTGGTTCCCGGCCCACAAACTGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                    <  2:108641/143‑1 (MQ=255)
                                           gtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                    >  1:548732/1‑143 (MQ=255)
                                             ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGAGAGCGGGATAACAATCAGGCTTGACGTTA‑‑ACACCGCCGTATGGCCAAGCAGCACAGCCGGATGAAGGCCAATCAGGCACGCCt                                                                                  >  2:274301/1‑143 (MQ=255)
                                             ggCAAGACCGCCAAACTGCAGCACCTGATTACCCGGCACACAAGCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                  <  2:803408/143‑1 (MQ=255)
                                               cAAGATCGCCAAACGGCAGCACCTGATGACCCGGCACACAAACAGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                                <  1:330931/143‑1 (MQ=255)
                                                aaGATCGCCAAACGGCCCCACCTGATTCCCCGTCACACAAACAGCACTCAAAACGGTGCGTCGCATAAACATCAGGCTCGCCGATACCA‑‑CCGCCGTATCTCCCAGCAGCAACGCCGGATCCATGCCACTCCGCGACTCCTgcc                                                                               >  2:468688/1‑141 (MQ=255)
                                                aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCACTCAGGAACTCCTGAc                                                                               >  2:476577/1‑143 (MQ=255)
                                                         aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCACCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCccgcgccacc                                                                      >  2:719542/1‑137 (MQ=255)
                                                           aCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                    <  1:967944/143‑1 (MQ=255)
                                                            cGGCAGCACCTGATTACCCGGCACACAACGAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                   <  2:802213/143‑1 (MQ=255)
                                                            cGGCAGCACCGGACCACACGGCACACAAACAGCAATTAATCTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTTCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                   <  1:806655/143‑1 (MQ=255)
                                                             ggCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                  <  1:891439/143‑1 (MQ=255)
                                                              gcagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCAACGCCAGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGcc                                                                 >  2:70701/1‑143 (MQ=255)
                                                               cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGGCt                                                                >  1:267572/1‑143 (MQ=255)
                                                                 gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGACTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                              >  1:120532/1‑143 (MQ=255)
                                                                     cTGATTACCCGGCACACAAACAGCAATTAAACTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCGGTAAAc                                                          >  2:246537/1‑143 (MQ=255)
                                                                     cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                          <  1:288837/143‑1 (MQ=255)
                                                                              cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                 >  1:610926/1‑143 (MQ=255)
                                                                               ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCCGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                <  2:520213/143‑1 (MQ=255)
                                                                                  acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACCA‑‑CCGCCGTATGTCGCAGCAGCAACGCCGGATCAAGGCCGATCAGGCACTCCTGACCGCCGAACTGCGCCCGTACacgacgacgcacccg                                             >  2:857437/1‑137 (MQ=255)
                                                                                   cacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAGCTCCTAACCGCCCAACTACGCCTGTAAACGACGACGAGCCTGc                                            >  2:1017358/1‑143 (MQ=255)
                                                                                      aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CAGCCGTATGTCACAGCAGCAACGGCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCCTa                                         >  2:169988/1‑143 (MQ=255)
                                                                                           gCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                    <  1:713020/143‑1 (MQ=255)
                                                                                             aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTATAAACGACTACGAGCCTGCTTAGCGATgg                                  <  1:282202/143‑1 (MQ=255)
                                                                                                          aGTTGGATAAAAATCACGCTTGCCTATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCACGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                     <  1:94979/143‑1 (MQ=255)
                                                                                                          aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                     >  1:210864/1‑143 (MQ=255)
                                                                                                          aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                     <  1:63367/143‑1 (MQ=255)
                                                                                                            tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCCTGCTTAGCGATGGGCGTTAAAgcatcca                   >  1:969721/1‑143 (MQ=255)
                                                                                                                       tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTa        <  1:338477/143‑1 (MQ=255)
                                                                                                                       gcAGGCTTGCCGAGAGCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTa        <  2:515300/142‑1 (MQ=255)
                                                                                                                        gaGGCTTGCCGAGAC‑‑CGCCGCCGTACGTCCCCGCAGCAGCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa       <  2:587661/142‑1 (MQ=255)
                                                                                                                             ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc  >  1:150792/1‑143 (MQ=255)
                                                                                                                                         ||                                                                                                                                   
TCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACC  >  NC_000913/2173224‑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: