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.4% 127.2 / 5.0 35pseudogene (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 . (15/17);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
*NC_0009132,173,3620C.94.1% 125.2 / 1.9 35pseudogene (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/0);  major base . (15/17);  minor base G (0/2);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACC  >  NC_000913/2173222‑2173493
                                                                                                                                           ||                                                                                                                                   
gttCCACGATGTACGGCGATCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGAACCGGAGTACCCGGCACACAAACAGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                                                                                                   <  1:244349/142‑3 (MQ=255)
   ccACGATGCGCGGCGACCGCCATCGCGCCGAAGAAGCCGATGGGGGCAAGATCGCGAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt                                                                                                                                <  2:204991/143‑6 (MQ=255)
             ggccGACCGCCTTCGCGACGATGAAGCCGATGGTGGCAAGGTCGCCAAACGGCAGCACCTGAGTACCCGGCACACGAACAGCGATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagc                                                                                                                    <  2:275768/140‑1 (MQ=255)
                   ccGCCATCGTCACGAAGAAGCTTATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                              <  1:213344/143‑1 (MQ=255)
                    cGCCACCCCCACGACGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCAGCTGATTACCCGGCACACCAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                                                                                             <  1:376398/143‑1 (MQ=255)
                           gccagcgagacgCCGATGGTGGCAAGATCGGCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAGTTAAAATGGTGAGTGGGATAAAAATCAAGCTTGCCGATACCA‑‑CCACCGTCTGTTCCTGCAGCAACGCCGGATc                                                                                                      <  2:702278/136‑1 (MQ=255)
                                      gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGGATGTCCCAGCAGCAACGCCCGATCAAGGCCAATCa                                                                                           >  2:324924/1‑143 (MQ=255)
                                       ccgagtgtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                          <  2:229725/137‑1 (MQ=255)
                                        cGATGGCGGGAAGATCGCCAAACGGCAGCACCCGCTGCGTCTTCACACAAACAGCAACTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                         <  1:485426/143‑1 (MQ=255)
                                         gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGACCAACGCCAACCAGGa                                                                                        >  2:516044/1‑143 (MQ=255)
                                               ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑ACGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTccg                                                                                  >  1:482237/1‑142 (MQ=255)
                                                gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                 >  1:419715/1‑143 (MQ=255)
                                                    gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                             <  1:471909/143‑1 (MQ=255)
                                                     aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAACAGCATAAAAAGCAGGCTTGCCGATACCA‑‑CCGCAATATGTCACAGCAGCAACGCCGGATCAAGGCCAAGCAGAAAATCCtcgcccc                                                                            >  2:118714/1‑137 (MQ=255)
                                                       cgccgggcGGCAGCGCTGGATTACCCGGCACAGAGCCGGCGAGTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                          <  2:2725/136‑1 (MQ=255)
                                                           gaaCGGCAGCACCTGATTACCCGGCACACAAGCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAGCGCCGGATCAAGGCCAATCAGTAACTCCTGACCGCCGAACt                                                                      <  2:88704/142‑1 (MQ=255)
                                                            aaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGCATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAACAGCAACGCCGGATCAAAGCCAATCAGAAACTCCTGACCGCCGAACtt                                                                     >  2:397972/1‑143 (MQ=255)
                                                             aCGGCAGCACAGGATGAGCCGGCACACAAACAGCAATTGAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTGCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                    <  2:248879/143‑1 (MQ=255)
                                                                  agcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg                                                               >  1:576298/1‑143 (MQ=255)
                                                                       cTGATTGCGCGCCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                          <  1:634439/143‑1 (MQ=255)
                                                                       cTGACTACCCGGCACACAAACAGCAATGAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                          <  2:23741/143‑1 (MQ=255)
                                                                                cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGGTAATCAGAAACCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGaa                                                 >  1:160675/1‑142 (MQ=255)
                                                                                       caAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCGTGCtt                                          >  1:72532/1‑143 (MQ=255)
                                                                                                aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGCTggg                                 >  1:105584/1‑143 (MQ=255)
                                                                                                 ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAACAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTCGCGATGGGc                                >  2:302788/1‑143 (MQ=255)
                                                                                                 gtgaagtGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                <  2:576298/137‑1 (MQ=255)
                                                                                                   aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                              >  1:199972/1‑143 (MQ=255)
                                                                                                        ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGCCGGGCGCTCAAc                         >  2:706103/1‑143 (MQ=255)
                                                                                                            aGTGGGATAAATATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                     <  1:382799/143‑1 (MQ=255)
                                                                                                             gtgtgggaaaacTCAGGCTTGCCGACACCA‑‑CCTCCGCATGTCCCAGCAGCAACGCCGGATCAAGGCCATTCAGGAACTCCTGACCGCCGAACTGCGCCTGTAATCGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcc                    <  2:511480/136‑1 (MQ=255)
                                                                                                             gTGGGATAAAAATCATACTTGCCGATCCCTCCGCCGTATGTCCCAGCTGCCACGCCGGATCAAGGCCAACCTGGATCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcc                    <  2:229418/143‑1 (MQ=255)
                                                                                                                   tAAACATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCATTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg              <  1:359686/143‑1 (MQ=255)
                                                                                                                    aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg             <  2:163271/143‑1 (MQ=255)
                                                                                                                       aaTCAGGCCGGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt          <  2:72532/143‑1 (MQ=255)
                                                                                                                           agccTTGCCGAGACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTGCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt      <  2:579703/140‑1 (MQ=255)
                                                                                                                             gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa    >  1:466831/1‑143 (MQ=255)
                                                                                                                               ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGATGGTACCTTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTCAACCATCCGTGAGGCGTTTGATCAcc  >  1:366077/1‑143 (MQ=255)
                                                                                                                                           ||                                                                                                                                   
TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACC  >  NC_000913/2173222‑2173493

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: