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.0% 71.8 / 8.0 23intergenic (‑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 . (13/7);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 6.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01
*NC_0009132,173,3640C.87.0% 72.0 / 7.9 23intergenic (‑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 . (13/7);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 6.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

GCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173132‑2173613
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                         
gccAGTTGGGTGTACAAACAAATAGATTGCGTCGCGATCCACAGGGTGATGCTCCTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCCG‑ATTACCCGGCCCACAAACAGCAATTACAATGGTG‑CGTGGGCTCAAAATCAGG‑CTTGCCGACACCACC‑‑GCCGTAAGTCCcagcacca                                                                                                                                                                                                                                        >  2:379228/1‑247 (MQ=255)
     ttGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCg                                                                                                                                                                                                                                  >  1:209222/1‑251 (MQ=255)
         gtgtGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                              >  1:117811/1‑251 (MQ=255)
           gtgCAAAGCATTAGTTTGCGTCGCGATCCACAGGGTGATGCTCCTAATGATGACACCCGTGCTTAAGGTGCGGCACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACCAACAGCAATTAAAATGGTG‑GGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTCTGTCCCAGCAGACACGCCgggccc                                                                                                                                                                                                                             >  2:172686/1‑246 (MQ=255)
              cAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                                                         >  1:295877/1‑251 (MQ=255)
                aaCCAATAGTTTGCGTCGCGACCCACAGGGTGAGTCGCATAATAATAACCCCCGAGATTAAGGTGCGGAACAAATTTCCACGATGCACGGCGACCGCCAACGCCACGAAGAAGCCGATGGTGTCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACGGCATTTAAAATTGTG‑AGTGGGATAAAAATCAAG‑CTTGCCAATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                                                                                                                                       <  1:189447/251‑1 (MQ=255)
                             cGTCGCGACCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                                                                                                                          <  1:73859/251‑1 (MQ=255)
                                        acaGGGAGGTGCTCATAATGATGACACCCGAGATTAAGGTACGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCAGATGGTGGCAAGATCGCCAAACGACAGCACCTG‑ATTACCCGGCACACAAACAGCAATTCAAATGGTG‑AGTGGGATACAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGTATCAAGTCCAATCATGAACTCCTGAccgcacaac                                                                                                                                                                                               >  2:35501/1‑246 (MQ=255)
                                           gggTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                                            >  1:58176/1‑251 (MQ=255)
                                                         atgatgACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                              >  1:86093/1‑251 (MQ=255)
                                                              gaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                           >  1:212442/1‑249 (MQ=255)
                                                                            aggggcgaaaaaGATTCCCCCGATGCACGGCGCCCCCCACCGCCACGAAGAAGCCAATGGTGGCAAGATCGCAAAACGGCAGCACCTA‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATg                                                                                                                                                           <  2:18601/240‑1 (MQ=255)
                                                                                     aaagattcccgcgaggcccGGCGACCGCCATCGCCCCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTTCCCGGCCCAAAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTaaa                                                                                                                                                  <  1:362455/233‑1 (MQ=255)
                                                                                              catgATGCACGGCGACCGCCATCGCCACAAAAAAACCCATTTTGCCAAAACCGCCAAACGTCAGCACCGG‑TTTCCCCGGCACACAAACAGCAATAAAAATGTGG‑ATTGGGATAAAAATCAGG‑TTTGCCGATTC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCTGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                                         <  1:42128/248‑1 (MQ=255)
                                                                                                             ccGCCATCCCCCCGAAGAAGCAGAGGGGGGCAAAAACGCCAAAGGGCAG‑ACCTGGATTACCCGGCACACAAACAGCAATTAAAAGGGTG‑AGTGGGATAAAAAACAGG‑CTTGCCGA‑ACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc                                                                                                                          <  2:212442/251‑1 (MQ=255)
                                                                                                                      ccACGAAGAAGCCGATGGTGGAAAGACCGCAAACCGCCAGCCCCTG‑ATAACCCGGCACAAAAACAGCAAT‑AAAAAGGTGGAGTGGGAAAAAAACCAGG‑CTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                                                 <  2:117811/251‑1 (MQ=255)
                                                                                                                              aagCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCAt                                                                                                         <  1:147292/251‑1 (MQ=255)
                                                                                                                                                 cgcaaaaCGGCAACACCTA‑ATTCCCCGCCACAAAAAAAGCAATTAAAATGGTG‑A‑TGGGGAAAAAAACAGGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                                      <  2:295877/247‑1 (MQ=255)
                                                                                                                                                          gcagcaCCTG‑ATTACCCGGTACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTgg                                                                                                                     >  1:202652/1‑211 (MQ=255)
                                                                                                                                                            agcaCCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCa                                                                            >  1:140730/1‑250 (MQ=255)
                                                                                                                                                                          ccGGCACAAAAACAGCAATAAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACa                                                              <  1:89213/251‑1 (MQ=255)
                                                                                                                                                                                       aGCAAATAAAATGGTG‑AGTGCGATACAACTCAGG‑CTCGCCGCTACCACC‑‑GCCATATCTCCCAGCAGCACCGCCGGATCAAGGCCCATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAAGCGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTCGCATTAGCAGCATCACTCCCGCCGTTTTTACAGCCAGCTGCACTACACCTTTGACATCGTAACCcgcg                                                  >  2:205418/1‑250 (MQ=255)
                                                                                                                                                                                                 aTGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                       >  1:256431/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCACCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAAGCAGCCCCATCACAAAACCCACGGTGACAgg      <  1:379228/247‑1 (MQ=255)
                                                                                                                                                                                                                                       cccccGCCGTATGCCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGACCGCCTGACCGCCGAACTTCGCCTGTAAACGACTACTAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTGAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTGTGACATCGTAGCCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGc     <  2:256431/250‑1 (MQ=255)
                                                                                                                                                                                                                                          ccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCGTTGACATCGTAACCCTCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCg  <  2:58176/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                         
GCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAATCAGG‑CTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173132‑2173613

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