Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,361 Δ2 bp 100% pseudogene (917‑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,3610C.88.1% 147.5 / 12.5 42pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base . (19/18);  total (21/21)
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
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*NC_0009132,173,3620C.88.3% 145.7 / 5.8 42pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base . (19/18);  minor base G (0/2);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173221‑2173497
                                                                                                                                            ||                                                                                                                                       
aTTTCCTCGCGGCGCGGCGGCCGCCATCGCCACTAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                        <  2:315471/143‑2 (MQ=255)
   tCCGCGATGCACAGCGACCGCCATCGCCACGAAGGAGCCGATGGTGGCCAGATCGCCAAACGGCAGCACCTGACTACCCCGCACACAAACAGCAATGAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccg                                                                                                                                     <  2:374190/143‑5 (MQ=255)
   tCCACGATGCGCGGCGACTGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCAGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccg                                                                                                                                     <  2:72947/143‑5 (MQ=255)
     cACGATGCACGGCGACCGCCATCGCCACGAAGAAACGGATTATGGCAAGATGGCCAAACGGCAGCACCTGATTACCCGGCACGCAAACAGCAATTGAAATGGTGAGTGGGATAAAAGTCAGGCTTGACGATaccaccgccgta                                                                                                                                   >  1:485343/1‑137 (MQ=255)
      aCGATGCACGGCGACCGCCATCGCCACCAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                                  <  1:802835/143‑8 (MQ=255)
        gctGCACGCCCACCGACATCGCCACGAAGCAGCCGAAGGTGGCAAGCTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                              <  2:6822/141‑1 (MQ=255)
             aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                         <  1:792870/143‑1 (MQ=255)
              cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagc                                                                                                                        <  2:828273/143‑1 (MQ=255)
              cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGCGAGCGGGATAACAATCAGGCTTGCCGACAC‑‑CCCCGCCGTACGCCCcagc                                                                                                                        >  1:383491/1‑143 (MQ=255)
              cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATAGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagc                                                                                                                        >  1:397172/1‑143 (MQ=255)
                gCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAGTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                      >  1:842761/1‑143 (MQ=255)
                 cGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                     <  1:327695/143‑1 (MQ=255)
                        cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                                                                                              >  1:453177/1‑143 (MQ=255)
                             ccACGAAGAAGCCGATGGTGGCAAGATCGCCAGACGGCAGCACCTGATTACCCGGCACACAAACAGCAATCAAAATGGTGAGAGGGATAAAAATCAGGCCGACCGCTA‑‑ACACCGCCATAGGTCACAGCACCAACGCCGGAtaa                                                                                                         >  2:491781/1‑141 (MQ=255)
                                cGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                      <  1:69587/143‑1 (MQ=255)
                                           aTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGAAAATAAAATGGTGACCGGACGAAGATTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTACCAGCAGCAACGCCGGCACAAGGCCGAGCAGaaa                                                                                           >  2:112293/1‑143 (MQ=255)
                                           aTGGTGGCAAGATCGCCAAACGGCAGCACCGGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAACAGTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                           <  2:205642/143‑1 (MQ=255)
                                               tggCAAGATCGCCAAACGGCAGCACCTGATTACCAGGCACACACACAGCAATTAAAATGGTGAGTGGGATACAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCCAGGAGGAACTcc                                                                                       >  2:157648/1‑143 (MQ=255)
                                                ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                      <  1:76645/143‑1 (MQ=255)
                                                     gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                 >  1:647494/1‑143 (MQ=255)
                                                      aTCGCCAAACGGCAGCACCTGATTGCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                <  1:629565/143‑1 (MQ=255)
                                                          ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                            <  1:69352/143‑1 (MQ=255)
                                                             caCGGCAGCACCTGACTACCCCGCACACCATGAGCAATTAAAGTGGTGCGTGGGATAAAAATCGGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                         <  2:385782/142‑1 (MQ=255)
                                                                ggCAGCACCTGATTACCCGACACACAAACAGCAATAAAAAGGGTGAGTGGGATAAAAATCAAGCTTGCCGATACCA‑‑CCGCCGTAAGACCCAGCAGCAACGCCGGCGCAAGGCCACTCACGACCAACTGACCGCCGAACtaccc                                                                      >  2:292442/1‑139 (MQ=255)
                                                                gccagcaCCTGATTACCCGGCACACAATCAGCAATTAAAATGGTGAGTGGGATACAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      <  2:132382/141‑1 (MQ=255)
                                                                 gcagcaCCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                     <  1:647364/143‑1 (MQ=255)
                                                                  cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGACCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                    >  2:634535/1‑143 (MQ=255)
                                                                          gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                            <  1:438505/143‑1 (MQ=255)
                                                                            ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      >  1:369681/1‑131 (MQ=255)
                                                                            ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      <  2:369681/131‑1 (MQ=255)
                                                                                ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                      <  1:730168/143‑1 (MQ=255)
                                                                                ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAACTCCTGACAGCCGAACTTCGCCCGTAAACGACTACGa                                                      >  2:58330/1‑143 (MQ=255)
                                                                                 cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAAGACTACGAg                                                     >  1:463100/1‑143 (MQ=255)
                                                                                      cacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCCATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGc                                                >  2:42589/1‑143 (MQ=255)
                                                                                                aagtAAGATGGTGAGTGGGATACAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCCGCAGCACCGCCGGATCAAGGCCAGTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGTGCCTGCTTAGCGATgg                                      <  2:519048/140‑1 (MQ=255)
                                                                                                 aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑GAGCCGTATGTCCCAGCAGCAACGACGGAGCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGGGCCTGCTTAGCGATggg                                     >  2:889281/1‑143 (MQ=255)
                                                                                                 aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACACCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCATGTAAACGACTACGAGCCTGGCTAGCGATggg                                     >  2:484489/1‑143 (MQ=255)
                                                                                                    aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCgc                                  >  2:643502/1‑142 (MQ=255)
                                                                                                     aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCCATCAGGAACTCCTGACCACCGAACTTCGCCTGTAACCGCCTCCGCGCCTGCTTAGCGATGGGCggt                                 >  1:368788/1‑141 (MQ=255)
                                                                                                      acgggTGAGCGCGATAAAAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                                <  1:69417/140‑1 (MQ=255)
                                                                                                            gAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                          >  1:861365/1‑143 (MQ=255)
                                                                                                             cttGGGATAAAAATCCGGCTTGCCGAGATCA‑‑CCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                         <  2:889678/141‑1 (MQ=255)
                                                                                                                ggcaTAAAAAGCACGCTTGCCGATAC‑‑CTCCGCCGGATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                      <  2:613244/140‑1 (MQ=255)
                                                                                                                        aaTCAGCCTAGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGTATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt              <  1:404492/143‑1 (MQ=255)
                                                                                                                         aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt             >  1:147475/1‑143 (MQ=255)
                                                                                                                                ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGACCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATcccc      >  1:623828/1‑143 (MQ=255)
                                                                                                                                    cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg  <  1:634535/143‑1 (MQ=255)
                                                                                                                                            ||                                                                                                                                       
ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173221‑2173497

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: