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.86.7% 91.4 / 12.5 30intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base . (20/6);  total (21/9)
Fisher's exact test for biased strand distribution p-value = 6.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01
*NC_0009132,173,3640C.86.7% 93.9 / 10.0 30intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base . (20/6);  total (21/9)
Fisher's exact test for biased strand distribution p-value = 6.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGC‑GGAACAG‑ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173125‑2173606
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                  
atgaGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGC‑GGAACAG‑ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCca                                                                                                                                                                                                                                       <  1:172469/248‑1 (MQ=255)
                     cAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGC‑GGAACAG‑ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                                                  >  2:155660/1‑251 (MQ=255)
                     cAAACCAATAGCTTTCGTCGCGATACAAAGGGTGAAGATCATAATGATGACACCCGAGATTAAGGTGA‑GGAACAG‑ATTTCCCCGATGCACGACCTCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGCTCGCCAACCGGCAGCCCCTGATTACCCGGCACAAACACAGCAATTACAATGG‑TGAGGGGGATAAAAATCAAGCTTGCCGATACCACC‑‑GCCATCTGTCCCAGCACCAACGCCGGCTCAAGGc                                                                                                                                                                                                                  >  2:252062/1‑251 (MQ=255)
                         ccAATAGTTGGGGTCGCGACCCACAGGGTGATGCTCATAATGATAACGCCCAGGATTAAGGTG‑‑GGAACAGGATTCCCAAGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCGCCTGATTACCCGCCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAAACAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                                                                                                                              <  2:129116/251‑1 (MQ=255)
                              aGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGC‑GGAACAG‑ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                                                          >  2:240047/1‑250 (MQ=255)
                                  ggcgtctggttCCACAGGGTGAGGCTCATAATGAAGACCCCCGAGATAAAGGTGC‑GAAACAG‑ATTTCCACGTTGCACGGCGCCCGCCACCGCAACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAGGG‑TGAGTGGGATAAAAATCAGGCTAGCCGATTCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                                                                                                                                      <  1:155660/240‑1 (MQ=255)
                                                caGGGGGATGCTCATAATGATGACACCCGAGATTAAGGTGC‑GGAACAG‑ATTTCAAAGATGCACGGCGACCGCCATCGAAACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGACAGCACCTGATTAACCGGCACATCAACAGCAATTAAAATGG‑TGAGTGGGATACGAATCAGGATTGAAGATA‑‑ACACCGCCGTATGTCCCAGCTGCAACGCCGGATCAACGCCAATCAGGAACTCCTGAccgcagcact                                                                                                                                                                                       >  2:99390/1‑245 (MQ=255)
                                                  gggTGATGCTCATAATGATGACACCCGAGATTAAGGTGC‑GGAACAG‑ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                                     >  2:180228/1‑251 (MQ=255)
                                                  gggTGATGCTCATAATGATGACACCCGAGATTAAGGTGC‑GGAACAG‑ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                                     >  2:134259/1‑251 (MQ=255)
                                                    gTGATGCTCATAATGATGACACCCGAGATTAAGGTGC‑GGAACAG‑ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                                                                                                                                   >  2:283755/1‑251 (MQ=255)
                                                            aaaaattataacccccGAGATTAGGGTGG‑GGAACAG‑ATTCCCAGGTTTCAGGCCGCCCGCCATCGCCACGAAGAAGCCGAGGGGGGCAAGAGCGCAAAACGGCGGCACCTGATTCCCCGGCCCACAACCAGCAAT‑AAAATGGGTGGGTGGGAAAAAAATCAGGTGGGCCGCTAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                                           <  2:138746/238‑1 (MQ=255)
                                                                         ccagagaTTAAGG‑GCAGGAACAA‑ATTTCCACCATGCACGCCGCCCCCCATCCCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCACCACCTGATTACCCGCAACACAACCAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGCCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACTAGCCTg                                                                                                                                                              <  1:7620/248‑1 (MQ=255)
                                                                                    ggTGC‑GGAACAA‑ATTTCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                                                                                                                                   <  1:374957/251‑1 (MQ=255)
                                                                                                    ttCCACGATGCACGGCGACCGCCATCGCCACGAAGAATCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAAAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcc                                                                                                                                     >  2:60663/1‑251 (MQ=255)
                                                                                                       cACGATGCACGGCACCCGCCAGCTCCACGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTGCCCGCCACACAAACACCAATTAAAATGG‑TAAGTGGAAAAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACTACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                                  <  1:225883/251‑1 (MQ=255)
                                                                                                          gATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                                                               >  1:115736/1‑251 (MQ=255)
                                                                                                              cgggccacccccctccGCCACGAAGAAGCCGATGGGGCCAAAACCGCCAAACGGCAGCACCTAATTACCCGGCACACAAACAGCAATTAAAATGG‑TAAGTGGGATAAAAATCAGGCTTCCCAAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt                                                                                                                           <  1:62426/236‑1 (MQ=255)
                                                                                                                  gCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACAGCCGCACTTCGCCTGTAACCGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                                                       >  1:172884/1‑251 (MQ=255)
                                                                                                                    cacccccacccccACGAAGAAGCGGATGGGGGAAAGAGCGCCAAACGGCACCACTGGATAACCCGGCACACAAACAGCAATAAAAAGGG‑TAAGTGTGATAAAATCCAGGCTTGCCG‑‑ACCCCACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGATCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa                                                                                                                     <  2:241969/240‑1 (MQ=255)
                                                                                                                                            gATGGTGGCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑AGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACGACGAGCCTGCTTAGCGGTGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTg                                                                                             >  2:201601/1‑251 (MQ=255)
                                                                                                                                                                  ggCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                                       >  1:2213/1‑251 (MQ=255)
                                                                                                                                                                                   cGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAt                                                      >  1:167703/1‑251 (MQ=255)
                                                                                                                                                                                         acaAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAc                                                >  2:175115/1‑251 (MQ=255)
                                                                                                                                                                                              cAGCAATTCAAATCG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGACAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGCGGGGTTTAATCACCCGTGGCATTAGCCGCATCACTGCCGCCGTCTTTACCGCCAGCTGCAATACACCTTTCGCATCGTAACCcgcg                                           >  2:185056/1‑251 (MQ=255)
                                                                                                                                                                                                   aTTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGAt                                      >  1:206807/1‑251 (MQ=255)
                                                                                                                                                                                                       aaaTGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGCTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGCTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCg                                  >  1:293325/1‑251 (MQ=255)
                                                                                                                                                                                                              tGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc                            >  2:238840/1‑251 (MQ=255)
                                                                                                                                                                                                                        aaaaaTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc                   >  2:79823/1‑250 (MQ=255)
                                                                                                                                                                                                                         aaaaTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac                 >  2:10238/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ttGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATAAGCCCCATCACAAAACCCAcc       >  2:268218/1‑250 (MQ=255)
                                                                                                                                                                                                                                      tgccgatccccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACg       <  1:283755/241‑1 (MQ=255)
                                                                                                                                                                                                                                          gataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGTGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGCATGCCGATAATCAGCACCATCACAAAACCCACGGTGAc  >  2:81142/8‑251 (MQ=255)
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                  
ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGC‑GGAACAG‑ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173125‑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: ‑