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.1% 70.9 / ‑2.2 17pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (8/8);  total (9/8)
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.1% 70.0 / ‑2.2 17pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (8/8);  minor base A (1/0);  total (9/8)
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

TTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAA  >  NC_000913/2173297‑2173423
                                                                ||                                                             
ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                                >  1:2211661/1‑65 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                              >  1:620029/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                              >  1:558743/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                              >  1:3264734/1‑64 (MQ=255)
   cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                             >  1:1254841/1‑63 (MQ=255)
    ccGGCACACAAACAGCTATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                            >  1:827966/1‑62 (MQ=255)
    ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                            <  1:1246416/65‑4 (MQ=255)
    ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                            <  1:766318/65‑4 (MQ=255)
    ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                            <  1:2962185/65‑4 (MQ=255)
           gaanaCAGCAATTAAAATGGTGATGGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTc                                                   <  1:2141557/63‑1 (MQ=255)
               aCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                               >  1:1586517/1‑65 (MQ=255)
                 aacAATTAAAATGGTGAGTGGAATCAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                             <  1:2789298/63‑1 (MQ=255)
                            tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                  <  1:1267886/65‑1 (MQ=255)
                            tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                  >  1:3104331/1‑65 (MQ=255)
                            tGGTGAGGGGGATAAAAATCAGGCTTGCCGGTAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGa                                  <  1:2462742/65‑1 (MQ=255)
                             ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                 >  1:2115078/1‑65 (MQ=255)
                             ggTGAGTGGGATAAAAATCAGGCTTGCCGAT‑‑CGCACCGCCGTATGTCCCAGCAGCAACGCCGGAt                                 >  1:3221618/1‑65 (MQ=255)
                                  gTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                            <  1:2503983/65‑1 (MQ=255)
                                    gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                          >  1:2728554/1‑65 (MQ=255)
                                       aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                       <  1:2553730/65‑1 (MQ=255)
                                       aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATATCCCAGCAGCAACGCCGGATCAAGGCCAAt                       <  1:1636494/65‑1 (MQ=255)
                                            aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCCATCAGGa                  >  1:2120247/1‑65 (MQ=255)
                                            aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa                  >  1:770306/1‑65 (MQ=255)
                                            aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAACAGGa                  >  1:567019/1‑65 (MQ=255)
                                               cAGNCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt               <  1:2145565/65‑1 (MQ=255)
                                               cAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt               <  1:3283640/65‑1 (MQ=255)
                                               cAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt               <  1:1305972/65‑1 (MQ=255)
                                                 ggCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc             <  1:1994657/65‑1 (MQ=255)
                                                   cTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg           >  1:1004277/1‑65 (MQ=255)
                                                        cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc      >  1:1584676/1‑65 (MQ=255)
                                                              ccccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa  <  1:3370148/63‑1 (MQ=255)
                                                                ||                                                             
TTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAA  >  NC_000913/2173297‑2173423

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: