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.94.4% 68.0 / ‑3.5 18intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (12/5);  total (13/5)
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.93e-01
*NC_0009132,173,3640C.100.0% 74.3 / NA 18intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (13/5);  total (13/5)

CCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173120‑2173613
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                         
ccAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGACATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACACGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATg                                                                                                                                                                                                                                                   >  1:184201/1‑251 (MQ=255)
                  tGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGTTAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGAAACGCCgg                                                                                                                                                                                                                                 >  1:2052/1‑251 (MQ=255)
                                                atccacatggggttccTCAAAATGATGCCCCCCGAAATGAAGGGCCGGAACAGTTTCCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGAAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                                                                                                                                                   <  2:66544/236‑1 (MQ=255)
                                                         gTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                                                                                                                                              >  2:34632/1‑247 (MQ=255)
                                                         gTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                                                                                                                                              <  1:34632/247‑1 (MQ=255)
                                                            aTGCTCATCATGATGACAACCGAGATTACGGTGCGGAACAGATTTCCACGATCAACGGCGACCGCCATCGCCACGAAGAAGACGATGGTGGCCAGATCGCCAAACGGCAGCACATGATTACCCGGCACCCAACCAGCCATTCAAATGGCGCGTGGGATAACAATCAGGCTTGCCGATACCACC‑‑GCAGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAATCCTGACCGCCGACATTCGCCTg                                                                                                                                                                                       >  2:297521/1‑251 (MQ=255)
                                                                             acCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCAGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACGACGAACCt                                                                                                                                                                      >  2:168081/1‑251 (MQ=255)
                                                                                             cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCATGTTCGCCAAACGCCACCACCTGATAACCCGGAACACAAACAGCAATTAAAATGGAGAGCGGGATAAAAATCAGATTTGCCGATACCACAC‑‑CCGTATGTCCCAGCAGCAACGCCGGATACAGGCCAATCAGGACCTCCTGACCGCCGAACTACGCCTGTAAACGACAACGAACCTGCTTAGCGATGGGCGt                                                                                                                                                      >  1:278906/1‑251 (MQ=255)
                                                                                             cGGAACAAATTCCCACGATGCACGGCCACCGCCATCGCCACGAAGAACCCGATGGTTCCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGTTGAGTGTGATAAAAACCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTTTAAACGACTACGAGCCTGCTTAGCGATTGGCGt                                                                                                                                                      <  1:108344/251‑1 (MQ=255)
                                                                                              ggAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                                                                                                                                     >  2:84650/1‑251 (MQ=255)
                                                                                                   aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCGGTAACCGACTACGAACCTGCTACGCGATGGGCGTAAAAcc                                                                                                                                                >  2:215490/1‑251 (MQ=255)
                                                                                                                 cccggcgccccacaccGCCACAAATAAGCCGAGGGTGCCAAGATCGCAAAACGTCAGCACCTTATTACCCGGCACCCAAACAGAAATTAAAATGGTGAGGGGGAAAAAAAGGAGGTTGGCCTATACCACC‑‑TCCGTATGTCCCATAAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTTCTTAGCGATTGGCGTTAAACCATCCATGATGGGtt                                                                                                                                  <  2:257495/236‑1 (MQ=255)
                                                                                                                       gacagcaaTCGACACGAAGAAGCAGATGGGGGCCAGATCTCCAAACGGCAGCACCTGATTACCCGGCACACCAACCGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTACTGACCGCCGAACTTCGCCTGTAAACGACTCCAAGCCTGCTCAGCGATGGGCGTTCAACCATCCCTGATGGGATTAATcc                                                                                                                            >  2:273068/8‑250 (MQ=255)
                                                                                                                                                        aaGATAGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttt                                                                                           >  2:257920/1‑251 (MQ=255)
                                                                                                                                                             cGCCAAAGGGCAACACCTGAGTCCCCGCAACAAAAACACCAATAAAAATGGTGAGTGGAAAAAAAATCAGGCTTGCGGAGACCACC‑‑GCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATAAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTGAATCACCCGTGGCATGAGCAGCATCACTGCCGCCGTATTTAcc                                                                                      <  2:261957/251‑1 (MQ=255)
                                                                                                                                                                                                  aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTACCCcgcg                                                  >  1:160375/1‑250 (MQ=255)
                                                                                                                                                                                                                       gATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcccc                             >  1:47111/1‑250 (MQ=255)
                                                                                                                                                                                                                              aTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                                         >  1:55616/1‑151 (MQ=255)
                                                                                                                                                                                                                                               ccccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAgg      <  1:84650/249‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCg  >  1:234648/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ccgccgCATGGCCCAGCAGAAACGCGGGATCAAGGCCAACTGGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGGTAACCCATCCAGGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTTCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACTGTGACAGGCTCg  <  2:141803/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                         
CCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173120‑2173613

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: