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.91.3% 165.5 / 9.2 46pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base . (21/21);  total (23/23)
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.84e-01
*NC_0009132,173,3620C.97.7% 170.1 / 1.8 45pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base . (21/21);  minor base A (1/0);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCG  >  NC_000913/2173221‑2173495
                                                                                                                                            ||                                                                                                                                     
aTTTCCACGATGCACGGCGACTGCGTTTGCCATTAATAAGCTGTTTGTTGTAAGATCGCCAAACGGCAGCACCTGAGAACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                      <  1:216204/143‑2 (MQ=255)
  ttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGACAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAAGGTGAGTGGGATAAAAATCCAGCTTACCGATACCACcgcg                                                                                                                                    >  2:137388/1‑143 (MQ=255)
        gATTCACGTCGGCCTCCTTCGCGACGAATAAGTCTGTGGTGGCAAGATCGCCAAACGTCAGCACATGATTGCCCGGCGCACAAACAGCTATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                            <  2:67619/143‑1 (MQ=255)
               ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                     <  1:454203/143‑1 (MQ=255)
                    ccGCCATCGCCACGAAGAAGCGCATGTTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                                <  2:694779/143‑1 (MQ=255)
                       ccATCGCAACGAAGAAGCCGATGGTGGCAAGATCGCAAAACGGCAGCTTCTGTTTACTTGTTACATAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                             <  1:305317/143‑1 (MQ=255)
                         attgCCACGAAGAAGTCGATGGTGGTTAGATCGCCAAACGGCAGCATCTGATTACCAGGCACACAAACAGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                           <  1:537220/140‑1 (MQ=255)
                            gCCACGAAGGAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                        <  2:390110/143‑1 (MQ=255)
                             ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGGGTGCGACAAAAATCCGGCTTGCCGATACCA‑‑CCGCCCTATGTCCCAGCAGCAACGCCGGATCa                                                                                                       >  1:827759/1‑143 (MQ=255)
                             ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAtaa                                                                                                       >  2:30287/1‑141 (MQ=255)
                              cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                      >  2:449323/1‑143 (MQ=255)
                              cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTAG‑‑AAAAAAACACAGACGAATATCCCAGCAGCAACGCCGGATCaa                                                                                                      >  1:395999/1‑143 (MQ=255)
                               aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATAAGACTTGCCGATACCA‑‑CCGCCGAATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                     >  1:107038/1‑143 (MQ=255)
                                cGAAGAATCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCAGACAAACAGCAATTAAAAAGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGAAACGACGGATCAAgg                                                                                                    >  2:545447/1‑143 (MQ=255)
                                cGAAGAAGCCGTTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAGCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                    <  2:557008/143‑1 (MQ=255)
                                        ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAA‑‑CCGCCGTATGTCACAGCAGCAACGCAGGATCAAGGCCAATCAg                                                                                            >  1:266401/1‑143 (MQ=255)
                                         cggtggcggCAACAGCGCCAAACCGCAGCACCTGATTCCCCGGCACACAACCAGCAATGAAAATGGTGAGTGGTATAAAAATCAGGCTTGCCTTTACCA‑‑CCGCCGTATGTCCCAGCAGCATCGCCGGATCAAGGCCAATCAgg                                                                                           <  2:782407/140‑1 (MQ=255)
                                          gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGATTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATAAAGGCCAATCAGGa                                                                                          >  2:626974/1‑143 (MQ=255)
                                            tggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCAGATCAAGGCCAGTCAGGAAc                                                                                        >  1:23103/1‑143 (MQ=255)
                                                  cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCAAACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAAGAGCAACGCCGGAACAAGGACAAGCAGGAACACCCGa                                                                                  >  2:98332/1‑143 (MQ=255)
                                                      attgCCAAACGGCAGCACCTGATTACCCGGCACACAAACTGCTATTTATATTGTTATTGGTATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                              <  1:714214/140‑1 (MQ=255)
                                                         gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCACTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCATGACCGCCGa                                                                           >  1:722554/1‑143 (MQ=255)
                                                         gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAACAGGAACGCCTGACCGCCGa                                                                           >  2:613096/1‑143 (MQ=255)
                                                          ccACACGGCATCACCTTATTAGCCGGCACACAATCAGCGATTAACATGGTGAGTGGGATCAAAATCAGGCTTGCCGATACCA‑‑CCTCCGTATGTCCCCGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                          <  2:277993/143‑1 (MQ=255)
                                                          ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGACAATCAGGAACTCCTGACCGCCGaa                                                                          >  1:97459/1‑143 (MQ=255)
                                                            aaaCGGCACCACCTGATTACCAAGCAAACAAACAGCAATTAAAACGGTGAGTGCGATAAAAATCAAGCTTGCCAATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCAGGATCAAGGCCAATCAGGAACCCCTCACCGCCGAACt                                                                        >  2:188480/1‑143 (MQ=255)
                                                                ggCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAAGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGTAACTCCTGACCGCCGAACTTCGc                                                                    <  2:273751/143‑1 (MQ=255)
                                                                 tcatcatcTGATTCCCCGGCACACAAACAGCAGTTAAAATGTTTACTGGGCTATCGCCCAAGGAGACCGATAC‑‑CCGCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                   <  1:368201/136‑1 (MQ=255)
                                                                        cTGATGACCCGGGAAACAAACAGTAATTAAAGTGGTGAGTTGGATAAAAATCAGGCATGCCGATACCA‑‑CCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                            <  1:715219/143‑1 (MQ=255)
                                                                               cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAACAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCATGTAAAAGACTACg                                                     >  1:130311/1‑143 (MQ=255)
                                                                                ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCAGAACAAGGCCAAACAGGAACTACTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                    >  1:725738/1‑143 (MQ=255)
                                                                                 cGGCACACAAACAGCAATTAAAATGGTGAGCGGGATAAAAAGCATGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGTATCAAGGCCAAGCAGGAACGCCTGACCGCCGGACTTCGCCTGTAAACGACTACGAg                                                   >  2:280248/1‑143 (MQ=255)
                                                                                     acacaAACAGCATTTAAACTGGCGAGTGGGATACAAAGCAGGCTTGCCGACACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTCACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTg                                               <  2:503737/143‑1 (MQ=255)
                                                                                     acacaAACAGCAATTAAAATTGTGTGTGGGATAAAAATCAGGCTTGTCGATACCA‑‑CCGTCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTg                                               <  1:492405/143‑1 (MQ=255)
                                                                                       acaACCAGCACTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                             <  2:127322/143‑1 (MQ=255)
                                                                                         aaaCAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCAGATAAAAGCCAATCAGAAACACCTGACCGCCGAAATTCGCCTGTAAACGACTACGAGCCTGCTTa                                           >  1:616471/1‑143 (MQ=255)
                                                                                           ccAGCAAGTAATCTTGTGATTGGGATAAAAATCAGGCTTTCCGATTCCA‑‑CCGCCGTATGTCCCAGCAGCACCGCCGGATCATGGCCAATCAGGAACTCCTGACCTCCGAACTTCGCCTGTATACGACTACGAGCCTGCTTAGc                                         <  2:564293/142‑1 (MQ=255)
                                                                                                  ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGAAACAAAGCAAAACAGGAAATCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                  >  1:676582/1‑143 (MQ=255)
                                                                                                        gtgTGAGGGGGATAAAACTGGGGGCTCGCGATACCA‑‑CCTCTGTATGTGCCAGCTGCCACGGCGGATCAGGGCCATCCAGGAGCTCCTGGCCGGTGTAGCTGTTATCTGAACGGTTACGAGCCTGCTTAGCGATGGGCGTTaaa                            <  2:168879/141‑1 (MQ=255)
                                                                                                         ggTGAGTGGGATAAAAATCATGCTTGCCGATACCA‑‑CCGCCTTATGTTCCATCATCAACGCTTGATCAAGGCCAATCAGGAACGCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                           <  1:244944/143‑1 (MQ=255)
                                                                                                             aGTGGGATAATAATCAGGCTTGCCGTTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAGTCAGGAACTCCTGACCGCCGAACATCCCATGTATTAAGTTACTTGCCTGCTTAGCGATGGGCGTTAAAccatc                       <  1:791219/143‑1 (MQ=255)
                                                                                                                 ggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCACCTGTAAACGAATACGAGCATACTTAACGATGGGCGTTAAAACTTCCatg                   >  1:136790/1‑143 (MQ=255)
                                                                                                                    tacaaaTCAGGCTCGCTGACACCA‑‑CCGCCGTATGTCCCATCAGCAACTTTTGATCAAGTCCAATCAGGAACTCCGGACAGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                <  1:484709/140‑1 (MQ=255)
                                                                                                                     aaaaaTCAGGCTGGCCGATAC‑‑CGCCGCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg               <  1:189185/143‑1 (MQ=255)
                                                                                                                                ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTCAGCGATGGGCGCTAAACCATCCAAGATGGGATTAATCAcc    >  2:468765/1‑143 (MQ=255)
                                                                                                                                gtGCCGATACCA‑‑CCGGCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTGCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc    <  2:719118/142‑1 (MQ=255)
                                                                                                                                  gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTGGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGGTAAACCATACATGATGGGTTATATCACACg  >  2:41131/1‑143 (MQ=255)
                                                                                                                                            ||                                                                                                                                     
ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCG  >  NC_000913/2173221‑2173495

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑