Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,448 C→T intergenic (‑575/+834) 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,4480CT96.4% 81.5 / ‑6.8 28intergenic (‑575/+834)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (14/13);  total (15/13)
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.09e-01

ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGAT  >  NC_000913/2173204‑2173687
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
aTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGc                                                                                                                                                                                                                                         >  1:182645/1‑251 (MQ=255)
  tAAGGGGCGGAACAGATTCCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCg                                                                                                                                                                                                                                        <  1:251803/250‑1 (MQ=255)
    aGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATg                                                                                                                                                                                                                                     >  1:446987/1‑251 (MQ=255)
                  tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACAAGCTTGCTTAGCGATGAGCGTTAAAccatc                                                                                                                                                                                                                       >  1:20852/1‑251 (MQ=255)
                                   gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa                                                                                                                                                                                                      <  1:343933/251‑1 (MQ=255)
                                             cccACGAAGAAGCCGATGGGGGCAAGACCCCCAAACGGCAGCCCATGATTCCCCGGCACACAAACAGCAATAAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAt                                                                                                                                                                                            <  2:182645/250‑1 (MQ=255)
                                                   aagaatCCGGTGGTGGAAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGAAATAAAAATGGTGAGGGGGAAAAAAACCAGGCTTGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcag                                                                                                                                                                                      <  2:43879/251‑1 (MQ=255)
                                                      aagCCGATGGTGGCACGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCATGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGAAGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCAt                                                                                                                                                                                   >  1:417494/1‑251 (MQ=255)
                                                          cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTCATCACCCGTGGCATTAGCAGCATCACt                                                                                                                                                                               >  1:338762/1‑251 (MQ=255)
                                                              ggtgccacggtcgcccaccggctccaCCTGATTCCCCGGCCCACAACCAGCAATTAAAATGGTGATTGGGATAAAAATCATGCTTTCCGATCCCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCATTCAGGAACTCCTTACCGCCGAACTTCGCCTGTTTACGACTTCGAGCTTTCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgcc                                                                                                                                                                            <  1:46525/226‑1 (MQ=255)
                                                                tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAACC‑‑GCCGTATGTCCCAGCAGCAACGCCGTATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                                                                                                         >  1:289535/1‑251 (MQ=255)
                                                                                     gcaCCTGAATACCAGGCCCACAAACAGAAATATAAATGGTGAGTGGGATAAAAATCAGGCTTGACGATACAACC‑‑GCCGTATGTCAAAGCAGCAACGCCGGATAAAGGCAAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGCCTACGAGCTTGCTTAGCGATGGGCGTTAAACAATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGACGTTTTTACAGCAAGCTGCat                                                                                                                                                     >  2:410179/1‑249 (MQ=255)
                                                                                                   cgcaccaaaaCACCAATTAAATGGGTGATTGGAATAAAAACCATGTTTGCCGTTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCTCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTATCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATc                                                                                                                                      <  1:108070/244‑1 (MQ=255)
                                                                                                      acacaAACAGCAATTAAAATGGTGCGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGAAGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGATGCAATACACCTTTGACATCGt                                                                                                                                    >  1:329807/1‑250 (MQ=255)
                                                                                                        acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACA‑‑GCCGCATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCAGCACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATAACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGtcac                                                                                                                                 >  1:56418/1‑248 (MQ=255)
                                                                                                          aaaCAGCACTTAAAATGGTGCGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACACCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCAATGGGCGTTAAACCATCCCTGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAcac                                                                                                                               >  2:271068/1‑249 (MQ=255)
                                                                                                             aagcaatgaaatgggtaaggggaaaaaaaaCCAGGCTTGCCGACAC‑‑CCCCGCCGTATGCCCCAGCAGCAACGCCGGATCAAGGCCTATCAGGAACTCCTGACCGCCGAACTTCTCCTGTAAACGACGACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcg                                                                                                                            <  1:91947/226‑1 (MQ=255)
                                                                                                                     aaaaTGGTGAGTGGGATAAAAACCAGGCTGCCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                                                                                                    <  1:110256/251‑1 (MQ=255)
                                                                                                                             gagggggATAAAAAAAAGGCTTGCGGATACAACC‑‑GCCGTATGCCCCGGCAGAAACGCCGGATCAAGGCAAAGGAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCAGGATGGGTTTAATCACCGGGGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCa                                                                                                            <  1:196628/247‑1 (MQ=255)
                                                                                                                              aGTGGGATAAAACTCCTGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCCAGGCCAATCCGGACCTCCTGACCGCCGAACTTCGCCTGTCAACGACTACTAGCCTGCTTAGCGTTGGGCGTTCAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTTCCGCCGTTTTTACCGCCAGCTTCACTACACCTTTGACATCGTACCCCGCGTGGATGCCGCTAATCa                                                                                                            >  2:196628/1‑250 (MQ=255)
                                                                                                                                                     cGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTg                                                                                     >  2:136147/1‑250 (MQ=255)
                                                                                                                                                             caccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACACGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCt                                                                              >  1:370309/3‑251 (MQ=255)
                                                                                                                                                                cgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCa                                                                                                                                                                                 <  2:107414/149‑1 (MQ=255)
                                                                                                                                                                             agcagcAACGCCGGAGCAAGCCCAATCAGGAACACCTGACCGCCGAACGTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCaaa                                                              <  2:373687/251‑1 (MQ=255)
                                                                                                                                                                                    aCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTgcgc                                                       >  1:73684/1‑251 (MQ=255)
                                                                                                                                                                                       ccGGATCAAGGCCAATCAGGAACTCCTACCCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGTCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCGATCTCAAAACCCACGGTGACAGGCTCGCCAAATTGACCAAAGTTGCGCTGa                                                    <  2:30213/251‑1 (MQ=255)
                                                                                                                                                                                                                       ccgAACTTCGCCTGAAAACGACTCCGAGCTTGCTTCGCGATGGGCGTTAAACCATCAATGATGGGTTTAATAACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGAACCAAATTTGCGCTGAATATAGTCGGCGCTAACTTTAATTCGGTTAAc                    >  2:317412/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cGACTACGGGCTTGCTTAGCGAAGGGGGTTAAACCATCCACGATGGGTGTAATCACCCGTGGCATTATCCTCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCAATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATTTCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCgac  <  2:77696/251‑2 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGAT  >  NC_000913/2173204‑2173687

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: