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.91.7% 137.8 / 4.7 36pseudogene (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/3);  new base . (17/16);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
*NC_0009132,173,3620C.94.3% 135.2 / 2.4 36pseudogene (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/2);  new base . (17/16);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTA  >  NC_000913/2173219‑2173487
                                                                                                                                              ||                                                                                                                               
aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                                                                                                  >  1:335180/1‑143 (MQ=255)
       cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCATGCCGATaccaccgccgta                                                                                                                           >  1:325229/1‑137 (MQ=255)
               aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACCA‑‑CCGCCGTATGTCCcag                                                                                                                 >  1:315973/1‑143 (MQ=255)
               aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                 <  1:270215/143‑1 (MQ=255)
                cGGCGCCCGCCACCGCCACGAAGGCGTCAATGGGGGCGGGGTCGTCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagc                                                                                                                <  1:160116/143‑1 (MQ=255)
                cGGCGACCGCCATCGCCGCGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATCAACATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagc                                                                                                                <  1:168760/143‑1 (MQ=255)
                               ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGAACa                                                                                               >  2:638586/1‑143 (MQ=255)
                                  cGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGGGA‑‑CGCCGGACCACgg                                                                                            >  2:410638/1‑143 (MQ=255)
                                    aagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGcc                                                                                          >  1:680308/1‑143 (MQ=255)
                                         gCCGATGGTGGAAGGAGAGCGAAGCGGCAGCACCTGATTACCCGGCACACAAACAGCAATGAAAATGGCGAGTGGGATAAAAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCa                                                                                     <  1:391041/143‑1 (MQ=255)
                                          ccGATGGGGGCAAGCTCGCCGTACGGCAGCACCAGCCTCAACAGCGCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAg                                                                                    <  1:401053/143‑1 (MQ=255)
                                           cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCAGTATGTCCCAGAAGCAA‑‑CGCCGGATCAAGGCCAATCAgg                                                                                   >  1:751750/1‑143 (MQ=255)
                                           cGAGGGTGGCAAGATCGTCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAgg                                                                                   <  2:237908/143‑1 (MQ=255)
                                             gtggtgtGAAGATCGCCAAACCGCAGCACCGGATTCCCCGGCACACAAACAGCAATTAAAATGGCGAGTGGGATAAAATTCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGaa                                                                                 <  2:348615/142‑1 (MQ=255)
                                               ggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCACTCAGGAACt                                                                               >  2:321604/1‑143 (MQ=255)
                                                gtgaCACGATGGCCAACCGGAAGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCAGATCAAGGCCAATCAGGAACTc                                                                              <  2:869597/143‑1 (MQ=255)
                                                 tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACCCAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTcc                                                                             >  1:141294/1‑143 (MQ=255)
                                                  ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTTGTGTGCCGAGAGCGAGAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCt                                                                            <  1:875020/143‑1 (MQ=255)
                                                       gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                       >  2:388067/1‑143 (MQ=255)
                                                        aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCAGGATCAAGGCGAGTGAGCAACTCCTGAcggc                                                                      >  2:757520/1‑143 (MQ=255)
                                                                  gtcagcaCCCGACTCCCCGGCACACAAACAGCAAGTAAAATGGTGAGTGGGATCAAAATCACGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                            <  1:284760/141‑1 (MQ=255)
                                                                    cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCCGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCCATCAGGAACTCCTGACCGCCGCACTTCGccc                                                          >  2:845106/1‑142 (MQ=255)
                                                                    cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGACTAGCCGATACAA‑‑CCACTGTATGTAAAGCGCATTTCCGCGGATCCAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGccc                                                          >  2:857589/1‑142 (MQ=255)
                                                                         cgcgATTACCCGGAACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                     <  1:464224/140‑1 (MQ=255)
                                                                            gATTCCCCTGCACACGAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAGTTCGCCTGTAAACGa                                                  <  2:742777/143‑1 (MQ=255)
                                                                               tACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTCTGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGCCTa                                               >  1:306612/1‑143 (MQ=255)
                                                                                  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTCAACGACTACGa                                            >  1:621466/1‑143 (MQ=255)
                                                                                  ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGCATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                            >  2:845678/1‑143 (MQ=255)
                                                                                    ggccgccaAACAGCACGTAAAAGGGTGAGTGGGATAAACGTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                          <  2:670900/137‑1 (MQ=255)
                                                                                             aCAGCAATTAAAATGGTGAGTGGGAGAAAAATCAGGCTTGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                 <  1:314782/143‑1 (MQ=255)
                                                                                               aGCAATTAAAATGGTGAGTGGGCTCACAATCAGGCTTGCCTATACCA‑‑CCGCCATATGTCCAAACAACAA‑‑AGACGAATAAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTCGCGa                               >  1:611735/1‑143 (MQ=255)
                                                                                                 cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTACTCGACAACGAGCCTCCTTAGCGGTg                             >  2:415551/1‑143 (MQ=255)
                                                                                                               aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc               <  1:590217/143‑1 (MQ=255)
                                                                                                                      tAAAAACCAGGCTCGCCGCTACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGACGGACCAAGGCCAAGCCGTAACCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg        <  2:478067/143‑1 (MQ=255)
                                                                                                                       aaaaaTCAGGCTCGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg       <  1:757520/143‑1 (MQ=255)
                                                                                                                         aaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCACGGCCAATCCGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGACGGACGTTAAACCATCCACGATggat     >  2:693448/1‑141 (MQ=255)
                                                                                                                          aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt    <  1:213514/143‑1 (MQ=255)
                                                                                                                            tCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTa  <  1:415453/143‑1 (MQ=255)
                                                                                                                                              ||                                                                                                                               
AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAA‑‑CGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTA  >  NC_000913/2173219‑2173487

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: