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.95.7% 101.9 / ‑1.1 23pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (17/5);  total (17/6)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
*NC_0009132,173,3620C.95.7% 101.1 / ‑1.9 23pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (17/5);  total (17/6)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGC  >  NC_000913/2173297‑2173419
                                                                ||                                                         
ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                            >  1:1497080/1‑65 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:2973569/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:2599283/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccacag                                                          >  1:3679423/1‑63 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGCGGGATACAAATCAGGCTTGCCGATaccaccg                                                          >  1:3956877/1‑64 (MQ=255)
         acccaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATg                                                 <  1:321684/62‑1 (MQ=255)
         acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATg                                                 <  1:2151481/65‑1 (MQ=255)
           acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTc                                               >  1:3230141/1‑65 (MQ=255)
           acaAACAGCAATTAAAATGGTGAGTGGAATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTc                                               <  1:684889/65‑1 (MQ=255)
            caAACAGAAAGAAAAATGTTAATGGGGATAAAAACCGGTTTTCCCGATAC‑‑CCCCGCCGTATGTcc                                              <  1:772889/65‑1 (MQ=255)
             aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                             >  1:1532535/1‑65 (MQ=255)
             aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                             >  1:2912473/1‑65 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCGGATACCA‑‑CCGCCGTATGTCCcagaa                                         >  1:436986/1‑63 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                         >  1:2083716/1‑65 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                         >  1:2581389/1‑65 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                         >  1:30757/1‑65 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                         <  1:1704142/65‑1 (MQ=255)
                       tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                   <  1:3484770/65‑1 (MQ=255)
                        aaaaTGGTAAGTGGAATAAAAGTCAGTCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGc                                  <  1:1700927/65‑1 (MQ=255)
                            tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                              >  1:3205831/1‑65 (MQ=255)
                             ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                             >  1:2926023/1‑65 (MQ=255)
                                  gTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                        >  1:3465074/1‑65 (MQ=255)
                                   tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                       >  1:2539600/1‑65 (MQ=255)
                                    gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                      >  1:1368838/1‑65 (MQ=255)
                                         aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATCa                 >  1:3757482/1‑65 (MQ=255)
                                         aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                 <  1:1758956/65‑1 (MQ=255)
                                         aaaaaTCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                 <  1:1480360/65‑1 (MQ=255)
                                            aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa              <  1:3041195/65‑1 (MQ=255)
                                            aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa              >  1:1461576/1‑65 (MQ=255)
                                              tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc            <  1:3889247/65‑1 (MQ=255)
                                               cAGCCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt           <  1:2424298/65‑1 (MQ=255)
                                                  gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt        >  1:3518344/1‑65 (MQ=255)
                                                  gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt        >  1:3610271/1‑65 (MQ=255)
                                                  gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAACGCCAATCAGGAACTCCt        >  1:3509356/1‑65 (MQ=255)
                                                  gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAACGCCAATCAGGAACTCCt        >  1:2722066/1‑65 (MQ=255)
                                                        cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc  >  1:1785527/1‑65 (MQ=255)
                                                                ||                                                         
TTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGC  >  NC_000913/2173297‑2173419

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