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.84.0% 74.0 / 13.5 25intergenic (‑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 . (13/8);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01
*NC_0009132,173,3640C.88.0% 82.8 / 8.9 25intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (13/9);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

CAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGA‑GATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAG  >  NC_000913/2173121‑2173608
                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                                    
cAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGA‑GATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGt                                                                                                                                                                                                                                             >  1:231181/1‑251 (MQ=255)
      tAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGA‑GATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCcag                                                                                                                                                                                                                                        <  1:60967/250‑1 (MQ=255)
            ccAGTTGGGTGGGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGA‑GATTAAGGTGCGGAACAGATTCCCACGATGCACGGCGACCGCCAGCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAAc                                                                                                                                                                                                                                 <  1:19211/251‑1 (MQ=255)
                         cAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGA‑GATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                                                    >  1:280676/1‑251 (MQ=255)
                                aTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGA‑GATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                               >  1:151326/1‑249 (MQ=255)
                                                          gATGCTAAAAATAATGAC‑CCCGAAGATTAAGGGGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACAAAGAAGCCGATGGTGGCAATATCGCCAAAGGGCAGCACCTGATGACCCGGCACACAAACGGCAATTAAAATGGGGGGGGGGA‑AAAAATCAGGGCTTGGCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGATCTCCTGACCGCCGAACTTCGCCt                                                                                                                                                                                   <  2:167941/251‑1 (MQ=255)
                                                                   aaTGATGACACCCGA‑GATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGAc                                                                                                                                                                          >  1:280628/1‑251 (MQ=255)
                                                                           cacCCGA‑GATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGct                                                                                                                                                                  >  1:70153/1‑250 (MQ=255)
                                                                                              ggAACAGATTTCCACGCTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGACAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGCTCAAGACCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGAAGGGCGtt                                                                                                                                                >  1:245188/1‑251 (MQ=255)
                                                                                                   aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAACAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACAGACTACGAGCTTGCTTAGCGATGGGCGTTAAAcc                                                                                                                                           >  2:269142/1‑251 (MQ=255)
                                                                                                     aTTTCCACGATGCGCGGCGACCCCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGAAATTAAAATGGTGAGTGGGATAAAAACCAGG‑CTGGCCGTTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCGGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccat                                                                                                                                         <  2:151326/251‑1 (MQ=255)
                                                                                                     aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccat                                                                                                                                         >  1:143852/1‑251 (MQ=255)
                                                                                                      tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatc                                                                                                                                        >  1:242048/1‑251 (MQ=255)
                                                                                                       ttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAg                                                                                                                                                           >  1:260314/1‑231 (MQ=255)
                                                                                                                         ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAACTCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCATGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc                                                                                                                     >  2:216312/1‑251 (MQ=255)
                                                                                                                              aTCGCCCCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTgg                                                                                                                >  2:144119/1‑251 (MQ=255)
                                                                                                                                             ccGATGGTGGAAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGAAATTAAAATGGTGAGTGGGATAAAAAGCAGG‑TTTGCAGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCAc                                                                                                 <  2:319897/251‑1 (MQ=255)
                                                                                                                                              cGATGGTGGAAAGATCGCCAAACGCCCGCCCCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGG‑CTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACt                                                                                                <  1:4461/251‑1 (MQ=255)
                                                                                                                                                   gtgtcAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgc                                                                                           >  2:226483/5‑251 (MQ=255)
                                                                                                                                                    gggaaaGATCGCCAAACGGAAGTCCCTGATTACCCGGCACAAAAACATCAATAAAAAGGGGGAGGGGGAAAAAAACCAGG‑CTGGCCAATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCGTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                          <  2:278796/247‑1 (MQ=255)
                                                                                                                                                      gCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGt                                                                                        <  1:247553/251‑1 (MQ=255)
                                                                                                                                                                          caCCTGATTACCCGGCACACAAACAGCAATTAAAAGGGGGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATa                                                                    <  1:99731/251‑1 (MQ=255)
                                                                                                                                                                               gatcccccGACACAAAAACAACAATAAAAAGGGTGAGGGGGAAAAAAACCTGG‑CTTGCGGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATacac                                                                 <  2:242048/244‑1 (MQ=255)
                                                                                                                                                                                       ggcacaaaaaaacaaataaaaaTGGTGAGGGGGAAAAAAAACAGG‑TGGGCCGAT‑CCACCC‑GCCTTATGTCCCAGAAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCGTGCTTAGCGATGGGCGTTAAACCACCCATGATGGGTTTAATCGCCCGTGTCATTAGCAGCATCACTGCCGTCGTTTGTACCGGCAGCTGCAATACTCCTGTGACATc                                                       <  2:16936/233‑1 (MQ=255)
                                                                                                                                                                                                           aaTGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTCTGTCCCAGCAGCAACGCCGGATCCAGGCCCATCAGGCCCTCCTGACCGACGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTCAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTCACCCGCGAGGATGCCg                                    >  1:218690/1‑250 (MQ=255)
                                                                                                                                                                                                                                      gctgcccgaaccccccGCCGTATGTCCCAGAAGCAACGCCGAATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCAGCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCa           <  2:231181/239‑1 (MQ=255)
                                                                                                                                                                                                                                                ccccccGCCGTATGGCCCAGCAGCAACGCCGAATCAAGTCCAAGCAGGAATTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGTCGTTAAACCATCCATGATGGGTTTAATCACCCGTTGCATTAGCAGCAGCACGGCCGCCTTTTTTACCTCCAGCTCCAATACACCTTGGACGTCGTAACGCTCGAGGTTGCCGATTATCAGCCCCATCACAAAGCCTGCGGTGACAg  <  2:218690/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                                    
CAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGA‑GATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAG  >  NC_000913/2173121‑2173608

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