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.87.5% 76.1 / 5.4 24intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (15/6);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 4.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3640C.87.5% 75.9 / 5.6 24intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (15/6);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 4.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCT  >  NC_000913/2173125‑2173611
                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                                       
aTTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTccc                                                                                                                                                                                                                                             >  2:251803/1‑250 (MQ=255)
  tAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCcagc                                                                                                                                                                                                                                          >  2:6617/1‑251 (MQ=255)
         cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAAcc                                                                                                                                                                                                                                   >  2:9379/1‑250 (MQ=255)
                                                                               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)
                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                                       
ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCT  >  NC_000913/2173125‑2173611

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