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.92.3% 147.4 / 4.0 39pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base . (25/11);  total (27/12)
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.94.7% 145.9 / 1.8 39pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (25/11);  total (27/12)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAC  >  NC_000913/2173221‑2173492
                                                                                                                                            ||                                                                                                                                  
ctttCCACGCTGCACGCCGGCCGCCATCGCCACGAAGCAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAGACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                   <  2:665053/142‑2 (MQ=255)
aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                   <  1:132652/143‑2 (MQ=255)
 ttcccACGATGCACGTCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                                                                                                  <  2:109346/140‑3 (MQ=255)
  tgccACGAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                                                                                                 <  2:515398/141‑4 (MQ=255)
        gATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGGGCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAGTTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGGTAC‑‑CCCCGCCGTATGt                                                                                                                         >  1:203680/1‑143 (MQ=255)
                gCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCTCCCAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTACCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                 >  1:79394/1‑143 (MQ=255)
                   aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCCCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagcca                                                                                                              >  1:715176/1‑141 (MQ=255)
                          tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                       >  1:301080/1‑143 (MQ=255)
                              cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                   <  1:36997/143‑1 (MQ=255)
                               gcGAAGAAGCCGAGGGTGGCGAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                  <  1:548788/142‑1 (MQ=255)
                               aCGAAGAGGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACTGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                  <  1:401165/143‑1 (MQ=255)
                                     aagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACCAACAGCAATTAAAATGGTGAGTGGGACAAAAATCAGGCTTGCCCATACCA‑‑CCGCCGTATGGCCCAGCAGCCACGCCGGGTCGCGGcccac                                                                                            >  2:167248/1‑140 (MQ=255)
                                      acccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                           <  1:48793/141‑1 (MQ=255)
                                       gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCACTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATA‑‑ACACCGCCGTATGTCCCAGCAACAACGCCGGATCAAGACCAATCa                                                                                          >  1:517798/1‑143 (MQ=255)
                                        ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCTCACGGCCAAGCAg                                                                                         >  2:512746/1‑143 (MQ=255)
                                           aTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACACGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCGCAGCAGCAACGCCGGATGACGGCCAATCAGGaa                                                                                      >  1:560417/1‑143 (MQ=255)
                                                  cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTgc                                                                               >  1:323319/1‑142 (MQ=255)
                                                    aGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                             <  1:512746/143‑1 (MQ=255)
                                                     gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGACGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                            >  2:704041/1‑143 (MQ=255)
                                                     gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGAGCTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGAACAAGTTGAGTCGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                            >  2:164607/1‑143 (MQ=255)
                                                       tCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAGAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCGGCAACGCCGGACCAAGGCCAATCAGGAGCTCCTGAccgcc                                                                          >  1:107447/1‑143 (MQ=255)
                                                           cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGGTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                      >  1:716871/1‑143 (MQ=255)
                                                            aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCCGATCAAGGCCAATCAGAAACAACTGACCGCCGAAcc                                                                     >  2:624976/1‑142 (MQ=255)
                                                            aaaCGGCAGCACCTGATTACCCAGCACACAAACAGCAATTAAAATGATGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCACACCAGCAACGCCGGAACAAGGCCAATCAGGACCACCAGACAGCCAAAcc                                                                     >  2:464484/1‑142 (MQ=255)
                                                             aaCGGCAGCACCTGATTACCCGGCACCCAAACAGCAATTAAAATGGCGAGTGGGATACAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCTCACGGCCAATCCGGAACTCCTGAccgccgccctt                                                                    >  2:294153/1‑138 (MQ=255)
                                                                     caCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCACGCCCAATCAAGAACTCCTGACCGACGCGCTTCGccccca                                                            >  2:470238/1‑139 (MQ=255)
                                                                     caCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAGCTGCggcggga                                                            >  1:193669/1‑137 (MQ=255)
                                                                        cggATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCTAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                         <  1:548395/141‑1 (MQ=255)
                                                                        cTGATAATCCCGCACACATCCAGCAATTAAAATGGTGATTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                         <  2:445921/143‑1 (MQ=255)
                                                                                 cGGCATACAAACAGCAATTAATATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTTCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                <  2:502104/143‑1 (MQ=255)
                                                                                 cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCCGACCGCCGAACTTCGCCCGCACCCGACTACgcg                                                >  2:276227/1‑141 (MQ=255)
                                                                                  ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGCCTACGAGc                                               >  1:429142/1‑143 (MQ=255)
                                                                                        caAACAGCAATTAAAATGGTGAGTGGGATAAAAATAAGGCTTGCCGATACCA‑‑CCGACGTATGTCCCAGCAGAAACGCAGGATCAAAGCCGATCAGGCACTACTGACCGCCGAACTTCGACTGTAAACGACTCCGAGCGTGCta                                         >  2:8055/1‑142 (MQ=255)
                                                                                                aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACCACGAGCCTGCTTAGCGATgg                                 >  2:648998/1‑143 (MQ=255)
                                                                                                  ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                               <  1:658600/143‑1 (MQ=255)
                                                                                                    agaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                             <  2:134077/141‑1 (MQ=255)
                                                                                                        tGGGGAGTGGGATAAAAATCAGGCTTGCCGAGACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                         <  2:796486/143‑1 (MQ=255)
                                                                                                            gAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                     >  1:681888/1‑143 (MQ=255)
                                                                                                                     aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg            >  1:624669/1‑143 (MQ=255)
                                                                                                                       aaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGTACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCCATTGGCGTTACACTATCCATCATGCGt          >  1:430549/1‑143 (MQ=255)
                                                                                                                              gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa   >  1:321357/1‑143 (MQ=255)
                                                                                                                              gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTACTGACCGCCGAACTTCGGCTGTAAACGACTACGAGCCTGCTTAGCGATGGTCGATAAACCAACCCTGAGAGGTTAAATCa   >  2:479682/1‑143 (MQ=255)
                                                                                                                               cctacctaTACCACCC‑‑CCGCATGCCCCAGCAGCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGACGAACTTCGCCTTTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc  <  2:193669/136‑1 (MQ=255)
                                                                                                                                            ||                                                                                                                                  
ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAC  >  NC_000913/2173221‑2173492

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: