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.89.7% 100.9 / 4.8 29intergenic (‑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 . (21/5);  total (21/8)
Fisher's exact test for biased strand distribution p-value = 1.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
*NC_0009132,173,3640C.89.7% 100.2 / 5.5 29intergenic (‑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 . (21/5);  total (21/8)
Fisher's exact test for biased strand distribution p-value = 1.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

AGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTG‑CGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173122‑2173606
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
aGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTG‑CGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTc                                                                                                                                                                                                                                          >  2:347585/1‑251 (MQ=255)
   aTTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTG‑CGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTccc                                                                                                                                                                                                                                        >  2:332700/1‑250 (MQ=255)
                tGGGTGTGCAAACCAATAGTTTG‑CGTCGCGATCCACAGGGCGATCACCAACATCATGACACCCGAGATTACGGAGCCACACAGATTTCCACGATGCACGGCAACCGCCATAGCCACGAAGAAGCCGATGGAGGCAAGATCGCCAACC‑GGCAGCACCTGATTAACCGGCAAACAAACAGCAATTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                                                                                                                          >  1:433539/1‑251 (MQ=255)
                      gcaaaaCCAAT‑GTTTGGCGGCGCGACCAACAGGGTGATGCTCAAAATGAGAACACCCGAAATTAGGGGGCGGAACAGATTTCCACGATGCCCGGCGACCGCCACCGCAACGAAGAAGCCGATGGTGGCAAGAGCGCCAAAC‑GGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAAGGGTGAGGGGGATAAAAATCAGGCGTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                                                                                                                                    <  2:427812/248‑1 (MQ=255)
                         aaaCCAATAGTTTG‑CGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                                 <  1:331321/251‑1 (MQ=255)
                              aaTAGTTTG‑CGTCGCGATCCACAGGGTGATGCTCATAATGAGGACACCCGAGCTTAAGGTACGGAACAGATTTCCACGATGCACGGAGACCGCCATCGACACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCAGGCACACAAAACGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCCACGCCGGATCACGGCCAATcc                                                                                                                                                                                                            >  1:425065/1‑250 (MQ=255)
                              aaTAGGTTG‑CGTCGCGATAAAGAGGGTGATGCTCATAGTGATGACACCCGAGATTAAGGTGCGGAACAGATTTCAACGATGCACGGCGACCGCAGTCGCCACGAAGCAGCCGATGGTGGCAAGATCGCCTAAC‑GGCAGAACCTGCTTACCCGGCCCACAAACAGCAATTAAAATGGCGAGTGGGATAAAAGTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGactcaggcccatc                                                                                                                                                                                                             >  2:417201/1‑238 (MQ=255)
                                         gTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCCCGAAGAAGCCGATGGTGGCAAGATCGCCACAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                                                  >  1:58326/1‑251 (MQ=255)
                                                   acaGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                                                                                                                                        >  1:17669/1‑251 (MQ=255)
                                                      gggTGCTGCCCATAATGATGACAACCGAGATTAAGGTGCGGAACAGATTTCAACGATGAACGGCGACCGCCATCGCAACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGAATACACGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATAAGGCTTGCCGATACCACC‑‑GCCGTCTGTCCCAGCAGACACGCCGGATCAAGGCCAATCTGGACCTCCGTACCGCCGACCTTc                                                                                                                                                                                     >  2:424387/1‑251 (MQ=255)
                                                                  tAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACg                                                                                                                                                                          >  1:253913/1‑250 (MQ=255)
                                                                                  gaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                                           >  1:5865/1‑249 (MQ=255)
                                                                                    ttAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                       >  1:29768/1‑251 (MQ=255)
                                                                                         ggtggaaacgatttcccccaatccACGGCGACCCCCATCGCCACAAAGAAGCCGAGGGGGGCAAGATCCCAAAAC‑GCCAGCACCTGATTACCCGGCAAAAAAACAGAAATAAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTGCGCCTGTAAACGACTACGAGCCTGCTTAGCGATggg                                                                                                                                                  <  2:435403/228‑1 (MQ=255)
                                                                                             ggAACAGATTTCCACGATGTACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGCTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGCTAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                                                                                                                              >  2:406720/1‑251 (MQ=255)
                                                                                                        ccacaagccACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                   <  1:406720/243‑1 (MQ=255)
                                                                                                        ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAGGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                   >  1:313944/1‑251 (MQ=255)
                                                                                                                   ggcgccccccacccccACGAAAAAGCCGATGGTGGCAAGATC‑CCAAACCGGCAGCCCCTGATTACCCGGCACCCAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGCCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCTTCCATGATGGGTTTaa                                                                                                                        <  1:332700/237‑1 (MQ=255)
                                                                                                                      gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACC‑GGCAGCAACTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGCGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc                                                                                                                      >  1:233725/1‑250 (MQ=255)
                                                                                                                            cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCCGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg                                                                                                               >  1:56820/1‑251 (MQ=255)
                                                                                                                                                             gCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAACATCAGGCTTGCCGATACAACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGACCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGATTTTACCg                                                                              >  1:172561/1‑251 (MQ=255)
                                                                                                                                                               ccaacggggcgccccggATTACCCGGCACAAAAACAGCAATTAAATGGGTGAGTGGGATAAAAATCAGCCTT‑CCG‑GACCCCACCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCTTTAAACCATCCAGGATGGGTTTAATCACCCTTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGc                                                                             <  1:44075/235‑1 (MQ=255)
                                                                                                                                                                                 ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGATTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                                                           >  1:464867/1‑251 (MQ=255)
                                                                                                                                                                                       ggCACAAAAACAGCAATAAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGATCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATc                                                     <  1:206576/251‑1 (MQ=255)
                                                                                                                                                                                               aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTACCCcg                                             >  2:37917/1‑251 (MQ=255)
                                                                                                                                                                                                         aaaagggtggggggaaaaaaaaCCAGGCGTGCCGACAC‑‑CCCCGCGGTATGCCCCAGCAGCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                   <  2:5865/234‑1 (MQ=255)
                                                                                                                                                                                                            aTGGTGAGTGGGATAAAAATAAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                >  2:384708/1‑251 (MQ=255)
                                                                                                                                                                                                                   gTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCAGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATCCACCTTTGACATCGTAACCCGCCAGGATGCCGATAATCAGc                         >  2:331239/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGCTCCCGACCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                   >  1:449329/1‑128 (MQ=255)
                                                                                                                                                                                                                                            gataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAc  >  2:167065/8‑251 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
AGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTG‑CGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173122‑2173606

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