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.94.4% 76.6 / ‑1.3 18pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (11/6);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.94.4% 76.2 / ‑1.8 18pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (11/6);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTG  >  NC_000913/2173299‑2173414
                                                              ||                                                    
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                     >  1:1487363/1‑64 (MQ=255)
aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                     >  1:1506731/1‑64 (MQ=255)
 cccGGCACACAAACAGCAATTAAAATGGTGGGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                    <  1:2451212/65‑3 (MQ=255)
  ccGGCACACAAACAGGAATAAAAATGGTGAGTGGGAAAAAAATCAGGCTTGCCGATacccccgcc                                                   <  1:1925136/65‑4 (MQ=255)
  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                   <  1:3205595/65‑4 (MQ=255)
  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                   >  1:818177/1‑62 (MQ=255)
  ccGGCACACAAACAGCAATTAAAATGGTGAGTGCGATAAAAATCAGGCTTGCCGATaccaccgcc                                                   >  1:1742675/1‑62 (MQ=255)
  ccGGCACAAAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                   <  1:2581860/65‑4 (MQ=255)
    ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGt                                               <  1:89957/65‑1 (MQ=255)
       acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATg                                            >  1:1138832/1‑65 (MQ=255)
          caAACAGCAATAAAAAGGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcc                                         <  1:1760730/65‑1 (MQ=255)
           aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                        >  1:3586502/1‑65 (MQ=255)
           aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGAccc                                        >  1:2206389/1‑65 (MQ=255)
            aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                       >  1:1650975/1‑65 (MQ=255)
            aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                       >  1:1423509/1‑65 (MQ=255)
            aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                       >  1:3170875/1‑65 (MQ=255)
                 cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                  >  1:1596194/1‑65 (MQ=255)
                    ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                               <  1:2147436/65‑1 (MQ=255)
                    ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                               <  1:3444922/65‑1 (MQ=255)
                     tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                              <  1:1615010/65‑1 (MQ=255)
                      aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                             <  1:2113284/65‑1 (MQ=255)
                          tGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                         <  1:2504111/65‑1 (MQ=255)
                          tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                         >  1:1930123/1‑65 (MQ=255)
                             tGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCa                      >  1:1898087/1‑65 (MQ=255)
                             tGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCa                      >  1:2737918/1‑65 (MQ=255)
                                gTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                   <  1:2274080/65‑1 (MQ=255)
                                  gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCCAGGcc                 >  1:222481/1‑65 (MQ=255)
                                  gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                 <  1:2466099/65‑1 (MQ=255)
                                          aaTCAGGCTTGCCGATACCA‑‑CCGCCGCATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa         >  1:778644/1‑65 (MQ=255)
                                            tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc       <  1:123359/65‑1 (MQ=255)
                                            tCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc       <  1:2902389/65‑1 (MQ=255)
                                                 cTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg  >  1:318139/1‑65 (MQ=255)
                                                              ||                                                    
ACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTG  >  NC_000913/2173299‑2173414

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: