Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑490/+918) gatD ← / ← gatB galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3630C.82.1% 82.2 / 14.4 28intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base . (12/11);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 5.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01
*NC_0009132,173,3640C.82.1% 84.4 / 12.3 28intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base . (12/11);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 5.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173120‑2173606
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
ccAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATg                                                                                                                                                                                                                                            >  1:92874/1‑251 (MQ=255)
           cgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGCTAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCagcagc                                                                                                                                                                                                                                  >  1:231052/1‑250 (MQ=255)
             ccggtggggTGTGAAAACCATGAGTTTGCGTCGCGACCCAAAGGGTGATGCTCATAAGGATGCCCCCCAAGATAAAGGTGCGGAACAGATTTCCACGATGCACGCCGCCCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAAc                                                                                                                                                                                                                               <  2:503529/245‑1 (MQ=255)
                gtggggtggccAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAAGGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCCCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGAGCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGGTTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                                                            <  1:24065/241‑1 (MQ=255)
                  ggggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                                                                                                                          <  1:338965/250‑1 (MQ=255)
                                         cGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                                                                                                                   >  1:15563/1‑251 (MQ=255)
                                           gcgcgATCCACAGGGTGATGCTAAAAAGGAGGACCCCCGATATTAAGGTGGGGAAAAGATTTCCACGATGCACGGCGACCGCCAACGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGGGGGATAAAAAGCAGGCTCGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                                                                                                                                                 <  2:142110/250‑1 (MQ=255)
                                               gATCCACAGGGTGATGCTCATAATGATGACACCCGAGCTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGCTAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCACCGCCGGATCAAGGCAAATCAGGAACTCCTGTccgc                                                                                                                                                                                             >  2:161494/1‑251 (MQ=255)
                                                                     ataatgaCCCCCGAGATTAAGGTGCGGAACAGATTCCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGGGGCAAGACCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATGAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCTATAC‑‑CGCCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                       <  1:336915/248‑1 (MQ=255)
                                                                        atgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                    >  1:446887/1‑251 (MQ=255)
                                                                          gaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                    >  1:204072/1‑249 (MQ=255)
                                                                             acCCGAGCTTAAGGTGCGGACCCGATTTCCACGATGCCCGGCGACCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGCCAGCACCTGATTACCCGGCACACAAACAGCAATTCAAATGGTGAGTGGGATAAAAATCAGGCTTTCCGATACCACC‑‑GCCGTATGTCCCAGCAGCACCGCCGGATCAAGGCCCATCAGCAACCCCTGACCGCCGACCTTCGCCTGTACACGCCTACGCGCCt                                                                                                                                                               >  2:212387/1‑251 (MQ=255)
                                                                                         ggTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                                                                                                                                   >  1:24397/1‑251 (MQ=255)
                                                                                            gcgaaaCAGTTTTCCACGTTGCACGGCGACCGCCTTCGCCACGAAGATGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGTATAAAAAGCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCGGGAACTCCTGACCGCCGAACTTCGCCTTTAAACGACTACGAGCCTGCTTAGCGATGGGc                                                                                                                                                 <  1:200689/246‑1 (MQ=255)
                                                                                                  aaaatgtcccggaggcagggCGCCCGCCATCGCCACAAAGAAGCCGAGGGGGGAAAGATCGCAAAGCGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGGGAGGGGGATAAAAATCGGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                                                                                                                                          <  2:92874/233‑1 (MQ=255)
                                                                                                          cccgaggcccGGCGCCCGCCATCGCCACGAAGAAGCCGAGGGTGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                                  <  1:350028/242‑1 (MQ=255)
                                                                                                            ggagcccgggCGACCCCCAGCGCCACGAAGAAGCCGAGGGGGGAAATATCGCCAACCGGCAGCACCTGATTCCCCGGCACAAAACCGGCAATTAAAAGGGTGAGGGGGATAAAAATCAGGCTGGCCTATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                                                                <  2:24397/243‑1 (MQ=255)
                                                                                                                                        agaagCCGATGGGGGAAAGACCGCAAAACGGCAGAACCTAATGACCCGGCACCCAAAAAGCAATAAAAAGGGGGAGGGGGATAAAAACCGGGCTGCCCGACACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTTACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcagc                                                                                                    <  1:212387/251‑1 (MQ=255)
                                                                                                                                                       cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAGACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                                   >  1:137900/1‑125 (MQ=255)
                                                                                                                                                       cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑ACCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTGCTTATCGATGAGCGTTCAACCATCCATGATGGGTTTAATCACCCGTGGCATTATCCGCATCACTGCCGCCGtt                                                                                     >  2:356606/1‑251 (MQ=255)
                                                                                                                                                             cGCCAAACGGCAGCACCTGTTTCCCCGGCACAAAAACAGCAATAAAAAGGGTGAGTGGGAAAAAAATCAGGCTGGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                               <  1:78870/251‑1 (MQ=255)
                                                                                                                                                                                aTTACCCGGCACAAAAACAGAATTTAAAGTGGTGAGGGGAATAAAAA‑CCGGCCTGGCGA‑ATCCCACCGCCGGATGCCCCAGCAGCAACGCGGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGTCATTAGCAGCATCACGGCCGCCGTTTTTACCGCCGGCTGCAATACAccgt                                                            <  2:69212/251‑3 (MQ=255)
                                                                                                                                                                                    cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGTCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc                                                        >  1:444652/1‑251 (MQ=255)
                                                                                                                                                                                              aaaCAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAccc                                              <  1:381078/251‑1 (MQ=255)
                                                                                                                                                                                                  aGCAATTAAAATGGCGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTACCCcgcgc                                          >  2:202459/1‑250 (MQ=255)
                                                                                                                                                                                                        aaaaaTGGGGAGTGGGAAAAAAAGCAGTCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGACCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGc                                    <  1:161494/250‑1 (MQ=255)
                                                                                                                                                                                                         aaaaGGGTGGGTGGGAAAAAAACCAGGCTTGCCGA‑ACC‑CACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGATCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCTTTAAACCATCCATGATGGGTTTAATCGCCCGTGGCATTAGCAGCATCACTGCCGCCGTTTGTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                   <  1:509680/251‑1 (MQ=255)
                                                                                                                                                                                                                             aaTCAGGCTTCCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATTCACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACaa               <  1:44505/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gatacacccGCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGACTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGCTGCCGATAATCAGCCCCATCACAAAACCCCCGGTGAc  >  2:79505/7‑251 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
CCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173120‑2173606

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑