Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑1/+1) gatC ← / ← gatC pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.97.5% 166.2 / ‑3.1 40pseudogene (2/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (16/23);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01
*NC_0009132,173,3620C.97.5% 165.6 / ‑3.3 40pseudogene (1/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (16/23);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173227‑2173497
                                                                                                                                      ||                                                                                                                                       
aCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCCGGCTTGCCGATaccaccgccgtat                                                                                                                                  >  1:292558/1‑136 (MQ=255)
    tGCATGGCGACCGCCCTCGCCCCGAAGAAGCCGATGGCAGCAAGACCGCCAAACGGCAGCACCTGACGACCCGGCACACAAACAGCCATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcc                                                                                                                            <  1:642111/143‑1 (MQ=255)
     ccACGGCGACCGCCATCCCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                                                                                                           <  1:347212/142‑1 (MQ=255)
       aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                         >  1:811507/1‑143 (MQ=255)
       aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCAGGCACACAAACAACAATTAAAATGGTGAGTGGGATAAAAATCAGACTCGCCGCTACCA‑‑CCGCCATATGTACcag                                                                                                                         >  2:370025/1‑143 (MQ=255)
            gcccGCCAGCGCAACGAAGAGGCCGATGGGGGCAAGATCGCCAAACGGCAGCAGCTTATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                                                                                                    <  1:119041/141‑1 (MQ=255)
               cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                                                                                                 >  2:730620/1‑143 (MQ=255)
                         aCGTAGAAGCCGATGGTGGCAAGCCCGCCAAATGGCAGCGGCTGATTACTCGTCTGACACGCAGCAATTAAACTCGTGAGTGGGATAAAAATGAGGCCTGCCGATACCA‑‑CGGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                       <  2:139520/143‑1 (MQ=255)
                         aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                       >  1:698199/1‑143 (MQ=255)
                          cGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                      <  2:244159/143‑1 (MQ=255)
                               acgCCGATGGTGGCCAGATCGCCAAACCCCAGCACCAAATTCCCCGGCCCACAAACAGCAATTAAAATCGCGCGTCGGATAAAAATCAGGCTTGCCTATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                 <  1:418658/141‑1 (MQ=255)
                                     aTGGGGGCAAGATCGCCAAACGGCAGCACCTGATTAGCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                           <  2:566213/143‑1 (MQ=255)
                                         tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc                                                                                       >  1:276039/1‑143 (MQ=255)
                                              aGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                  >  1:687778/1‑143 (MQ=255)
                                                     cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGAACAGAAGGTGCTTTGTCTATCAGGAACTCCTGACCGCCGAAc                                                                           <  1:817029/143‑1 (MQ=255)
                                                         cGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                       <  1:730620/143‑1 (MQ=255)
                                                             agcaCCTGATCACCCGCCACACAACCAGCAATTAAAACTGTCAGTGGCATAAAAAACAGGCTTGCCGATACCA‑‑CCGCCGTAAGTCCCAGCCGCACCGCCGGATAACGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTg                                                                   <  1:838520/143‑1 (MQ=255)
                                                                    gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                            <  1:509580/143‑1 (MQ=255)
                                                                      tcaCCCGGCACACAAACAGCAATTAAAATTGTGAGTGGGATAAAAATTAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                          <  2:418377/141‑1 (MQ=255)
                                                                       tACCCGGCCCACAAACAGCAATTAAGATGGTGAGTGGGATCAACATCAGGCTTGCCGCTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                         <  2:672268/143‑1 (MQ=255)
                                                                        aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACGGCCGAACTGCGCCTAGACACGACTAc                                                        >  2:454273/1‑143 (MQ=255)
                                                                        aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCACAGCACCAACGCCGGATCCAGGCCCATCAGGCACCCCTGACCGCCCAACTTCGCCTGCACACGACTAc                                                        >  2:244830/1‑143 (MQ=255)
                                                                          ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTGCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAACACCTGACCGCCGAACTTCGCCTGTAACCGGCTACgc                                                      >  2:553127/1‑142 (MQ=255)
                                                                           cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAACATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTCCgcg                                                     >  2:582035/1‑141 (MQ=255)
                                                                            ggCACACAAACAGCAATTAAAATGGAGAGTGGGATAAACATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCGCAGCAGCAACGCCGGATCAAGGCCAATAAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGCCTACGacc                                                    >  2:86767/1‑141 (MQ=255)
                                                                                   aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCGATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACCACGAGTGTGCtgg                                             >  1:187849/1‑141 (MQ=255)
                                                                                            ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                    >  1:668195/1‑143 (MQ=255)
                                                                                            tcaaaaTGGTGAGCGGGCTAACACTGCGGCGTGCCGCTAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                    <  1:514734/141‑1 (MQ=255)
                                                                                            atAAAAAGGTGAGGGGGATAAAAAACAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGACCAGATAGAGCGGAGGTTTTCGGGTACCCCTAACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                    <  1:396726/142‑1 (MQ=255)
                                                                                              aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGGCTACGAGCCTGCTTAGCGATGGGCGt                                  <  2:208957/143‑1 (MQ=255)
                                                                                                        gttaaacaaaaaTCTGGCGGGCCCAAACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcc                        <  1:389053/136‑1 (MQ=255)
                                                                                                          ggtaTAAAAACCCGGCTCTCCGATACCA‑‑CCGCCTTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                      <  1:53362/140‑1 (MQ=255)
                                                                                                           ggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCCGTAAACAATCCatg                     >  1:530040/1‑143 (MQ=255)
                                                                                                                 agagcgggCTTGGCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAGCGCCGGATCATGGTCAGTCAGGTACTCCCGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt               <  2:410825/137‑1 (MQ=255)
                                                                                                                 aaaGCAGGGTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTATAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt               <  1:309078/143‑1 (MQ=255)
                                                                                                                  agtCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt              <  2:272124/141‑1 (MQ=255)
                                                                                                                  aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTAGCCCATAACCAACTACGCGCCTGCCTAGCCATGGGCGTTAACCCATCCACGGTgggca              >  2:558962/1‑141 (MQ=255)
                                                                                                                    tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTa            <  1:701777/143‑1 (MQ=255)
                                                                                                                         gctGCCGAGACCA‑‑CCGCCGGCTGTGCCGGCGGCTACGCCGGATCAAGGCCAATCAGCAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGGGCCTGTTAGGCAATGGTCGGTAAACCATCCATGATGGGTTTAATCAc       <  2:62819/141‑1 (MQ=255)
                                                                                                                          ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGACCAAGGCCAACCAGGAACTCCTGACCGCCGACGTTCGCCTGTAAACCACTCCGAGCCTGCTTAGCGATGGGCGTTAAACCATCCACGATAGGTTTGATCaac      >  2:684835/1‑141 (MQ=255)
                                                                                                                              cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg  <  2:687778/143‑1 (MQ=255)
                                                                                                                                      ||                                                                                                                                       
ACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173227‑2173497

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: