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.90.0% 79.6 / 3.5 20pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (10/8);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.90.0% 78.2 / 1.3 20pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (10/8);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGC  >  NC_000913/2173297‑2173419
                                                                ||                                                         
ttACCCGGCACACAAACAGCAATTCAAATGGTGAGTGGGATAAAAAGCAGGCTTGCCGATaccac                                                            <  1:3957177/65‑1 (MQ=255)
 tACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccacc                                                           >  1:3597552/1‑65 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:3241510/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:2196781/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:2312307/1‑64 (MQ=255)
  aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                          >  1:3138618/1‑64 (MQ=255)
    ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                        >  1:111789/1‑62 (MQ=255)
    ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                        >  1:4202715/1‑62 (MQ=255)
    ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                        >  1:1926774/1‑62 (MQ=255)
         gcacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATg                                                 <  1:437924/64‑1 (MQ=255)
         acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATg                                                 >  1:3440335/1‑65 (MQ=255)
         acacaAACAGCAATTAAAATGGTGAGTGGAATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATg                                                 <  1:3309758/65‑1 (MQ=255)
         acacaAACAGAAATTAAAATGGTGAGTGGGATAAAAATCAGGCTGGCCGATACCA‑‑CCGCCGTATg                                                 <  1:3183425/65‑1 (MQ=255)
                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATTCCA‑‑CCGCCGTATGTCCcagca                                         <  1:4776170/65‑1 (MQ=255)
                   cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                       >  1:2540656/1‑65 (MQ=255)
                       tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                   >  1:784037/1‑65 (MQ=255)
                            tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                              <  1:2715479/65‑1 (MQ=255)
                            tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                              <  1:938939/65‑1 (MQ=255)
                            tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                              <  1:4071536/65‑1 (MQ=255)
                            tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                              >  1:4514614/1‑65 (MQ=255)
                             ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                             <  1:1035449/65‑1 (MQ=255)
                             ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                             >  1:1330378/1‑65 (MQ=255)
                                 aGTGGGATAAAATTCGGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                         <  1:654888/65‑1 (MQ=255)
                                 aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                         <  1:2869601/65‑1 (MQ=255)
                                        tAAAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                  <  1:945542/65‑1 (MQ=255)
                                         aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                 <  1:3027174/65‑1 (MQ=255)
                                           aaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg               >  1:1183871/1‑65 (MQ=255)
                                            aaTCAGGCTTGCCGATACCACCC‑‑CCGTATGTCCCAGCAGCAAAGCGGGATTAAGGCCAATCAggg              >  1:3804018/1‑64 (MQ=255)
                                            aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa              >  1:4161158/1‑65 (MQ=255)
                                            aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa              >  1:1499667/1‑65 (MQ=255)
                                            aaTCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCGATCAGGa              <  1:4579528/65‑1 (MQ=255)
                                             aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa             >  1:1471935/1‑65 (MQ=255)
                                               cAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAGCt           <  1:4550974/65‑1 (MQ=255)
                                                 ggCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc         >  1:1958299/1‑65 (MQ=255)
                                                 cgCTTGCGGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc         <  1:921966/64‑1 (MQ=255)
                                                        cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc  >  1:4749149/1‑65 (MQ=255)
                                                                ||                                                         
TTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGC  >  NC_000913/2173297‑2173419

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