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.5% 63.8 / 3.5 19intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (16/1);  total (16/3)
Fisher's exact test for biased strand distribution p-value = 1.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01
*NC_0009132,173,3640C.89.5% 64.1 / 3.2 19intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (16/1);  total (16/3)
Fisher's exact test for biased strand distribution p-value = 1.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

AGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGC‑AAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGA  >  NC_000913/2173174‑2173605
                                                                                                                                                                                              ||                                                                                                                                                                                                                                                  
aGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGC‑AAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                                                                                                                                     >  1:304113/1‑251 (MQ=255)
      aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGC‑AAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg                                                                                                                                                                               >  1:103941/1‑251 (MQ=255)
                 gtcgacaccgtagtttaggttGCGGAACATATGTCCACGATGCACGGCGACCGCTACCGCCACGAAGAAGCCGATG‑TGGCCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCCCACAACCAGCATTTAAAAGGGTGAGTGGGATAAAAATCAGCCTTCCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                                                                                                                                     <  1:6832/230‑1 (MQ=255)
                           gagaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGACATCGCCACGAAGAAGCCGATGGTGGC‑AAGATCGCCAAACGGCAGCACCTGATTAACCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATAAGGCTTGCAGATACCACC‑‑ACCGTATGTCACAGCAGCAACGCCGG‑CTCAAGGCCAAACAGGACCACCTGACCGCCGAAATTCGCCttgcacagaatacgagcttgcct                                                                                                                                                          >  2:184733/1‑229 (MQ=255)
                                        ggAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGC‑AAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGtt                                                                                                                                             >  2:283816/1‑251 (MQ=255)
                                               aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGC‑AAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                                                                                                                                        >  1:337479/1‑201 (MQ=255)
                                                  tCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGC‑AAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatcc                                                                                                                                    >  1:343360/1‑250 (MQ=255)
                                                          cacggcgccccaccacccccaagaagaaGCCGAGGGGGGA‑AAGATCCCAAAACGGCAGCACCTGATTACCCGGCACAAAAACAGCAATAAAAAGGGTGAGTGGGATAAAAACCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt                                                                                                                           <  2:103941/229‑1 (MQ=255)
                                                                       cATCGCCACGAAGAAGCCGATGGTGGC‑AAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg                                                                                                              >  1:207113/1‑251 (MQ=255)
                                                                                  agaagCCGAGGGTGGA‑AAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACATCAATTAAAATGGTGAGGGGGAAAAAAACCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcagc                                                                                                   <  1:211870/251‑1 (MQ=255)
                                                                                                   aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GACGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGCTGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTCGCATTAGCAGCATCACTGCCGCCGttt                                                                                   >  2:240843/1‑251 (MQ=255)
                                                                                                      atctccAAACGGCAGCACCTGATTAACCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATAACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                                 >  1:293458/5‑250 (MQ=255)
                                                                                                      aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCACTAGCAGCATCACTGCCGCCGttttt                                                                                 >  1:293448/1‑250 (MQ=255)
                                                                                                                               cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc                                                       >  1:205659/1‑251 (MQ=255)
                                                                                                                                                   tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGc                                   >  1:335495/1‑251 (MQ=255)
                                                                                                                                                            gAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCTAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATactc                           >  1:197249/1‑247 (MQ=255)
                                                                                                                                                                        aaTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcaca               >  1:31706/1‑250 (MQ=255)
                                                                                                                                                                            aGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATACCGATAATCAGCCCCATCACAAAAcc          >  1:131572/1‑251 (MQ=255)
                                                                                                                                                                                    cGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTg   >  1:277979/1‑250 (MQ=255)
                                                                                                                                                                                     cccacccccccccCGTTTGCCCCAGCAG‑AACCCCGGGATCAAGGCAAACCAGGACCCCCTACCCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAACCCACCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTg   <  2:31706/244‑1 (MQ=255)
                                                                                                                                                                                      aataccccccGCCGTATCCCCAAGCCGCACCGCCGG‑AACAGGGCCAATCAGGACTTCCGGACCGCGGAACGTCGCCGGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCAGGATGGGTTTAAGCACCCGTGGCATTAGCAGCATCACGGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGa  <  2:207113/245‑1 (MQ=255)
                                                                                                                                                                                              ||                                                                                                                                                                                                                                                  
AGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGC‑AAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGG‑ATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGA  >  NC_000913/2173174‑2173605

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