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.91.3% 77.5 / 3.5 23intergenic (‑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 . (13/8);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01
*NC_0009132,173,3640C.91.3% 78.4 / 2.6 23intergenic (‑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 . (13/8);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGT  >  NC_000913/2173123‑2173603
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                               
gCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAACC‑‑GCCGTATGTcc                                                                                                                                                                                                                                      >  2:150462/1‑251 (MQ=255)
                   aatgCAAACAAATAGTCTGCGTCGCGATCAAAATGGGGATGCTCATAATGATGACACTCGAGATTAAGATGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGGCAAAC‑GGCAGCACCTGGTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGTGATAAAAATCAGGCTTGCCGATACCATC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                    >  2:145916/3‑250 (MQ=255)
                       aaaaccaatagtttggtgagagagccACAGGGTGATGCTCCTAAGGTTGACACCCGAGATTAAGGTGGGGAACAAATTTCCCCGATGCACGCCG‑ACACCCAGCGCCACGAAGAAGCCGATGGTGGCAAGTTCGCCAAAA‑GGCAGCACCTGTTTACCCGGCAAACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCGTGCCGATACCACC‑‑GCCGTATGTACCAGCAGCAACGCCGGATCAGGGc                                                                                                                                                                                                               <  2:225034/227‑1 (MQ=255)
                                       gTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                                               >  2:33455/1‑251 (MQ=255)
                                            gATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                                                                                                                          >  2:130027/1‑251 (MQ=255)
                                                     gtgtaagcacATAATGAAGACCCCCGAGATAAAGGAGAGAAACAGATTTCCACGATGAACGGCG‑ACCGCCATCGCCACGAAAAAGCCGATGGGGGCAAAAGCGCAAAAG‑GGAAGCCCCTGATGACCCGGCACAAAAACAGCAATGAAAAGGGGGAGGGGGATAAAAATCAGGCTTGCCGTTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCGACTGACCGCCTAACTTCg                                                                                                                                                                                 <  2:39381/242‑1 (MQ=255)
                                                       ggttctcagaaaTGATGACCCCCGAGATTAGGGTGCGGAACGGATTTCCCGGAGCCACCGCG‑CCCGCCACCGCCAGGAAGAGGCCGTTGGTGGAAAGATCGCCAAAC‑GGAAGCCCCTAATTCCCCGGCACACAACCAACAATAAAAATGGTGAGTGGAAAAAAAACCGGGCTGGCCATTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                                                                                                                               <  2:223555/241‑1 (MQ=255)
                                                              cATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGAACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                                        >  1:246364/1‑251 (MQ=255)
                                                                            ccGAGATAAGGGTGCGAACCAGATTCCCCCGATGC‑CGGCGGACCGCCACCACCACGAAGAAGCCGTGGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAAAAGCAATTAAAAGGGTGAGTGGGAAAAAAACCAGGCTGGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGc                                                                                                                                                          <  2:289449/251‑1 (MQ=255)
                                                                              gagaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                                                                        >  1:270790/1‑251 (MQ=255)
                                                                                       gTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATAGCCACGAAGAAGCCGATGGGGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATggg                                                                                                                                               >  1:32997/1‑251 (MQ=255)
                                                                                                              cACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt                                                                                                                        >  2:260008/1‑251 (MQ=255)
                                                                                                                           cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg                                                                                                            >  1:280492/1‑251 (MQ=255)
                                                                                                                               gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAt                                                                                                        >  2:56815/1‑251 (MQ=255)
                                                                                                                                         gccagatgggggcaatacacctaaaC‑AGCAGCACCTGTTTACACGCCACACAAAGAGTAATTAAAACGGGGAAAGGGAAAAAAATACGGCTTGCCAATACCACC‑‑GCCGTATGTCCCAGCACCAACGCCGGATCAAGGCCAATCAGGAACTCCCGACCGCCGAACTTCTCCTGTTAACGACTACGAGCTGGCTTAGCGATGGCCGTTAAACCAGCCATGATGAGTGTAACCACCCGTGGCATTAGCAGCATc                                                                                              <  2:123698/229‑1 (MQ=255)
                                                                                                                                                     aaaGATC‑CCAAACCGGCAGCACCTGATTCCCCGGCACACAAACAGAAATAAAAAGGGTGAGGGGGATAAAAACCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtt                                                                                  <  1:260008/250‑1 (MQ=255)
                                                                                                                                                            cccaaaccgccagcaCCTAATTACCCGGCAAACAAACAGCAATTAAAATGGGGAGGGGGATAAAAACCAGGCTGGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                            <  2:250071/241‑1 (MQ=255)
                                                                                                                                                                                    ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GACGTATGTCCCAGCAGCAAAGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAcc                                                    >  1:191653/1‑250 (MQ=255)
                                                                                                                                                                                             aaaCAGCAATTAACCTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAAAC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAccc                                           >  2:127307/1‑251 (MQ=255)
                                                                                                                                                                                                        aaaaGGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                <  1:237345/251‑1 (MQ=255)
                                                                                                                                                                                                           aTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                             <  1:130027/251‑1 (MQ=255)
                                                                                                                                                                                                            gggggAGTGGGATAAAAATCAGGCTGGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATa                            <  2:204647/247‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGCTTAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTACACGCCTTCGAGCTTGCTTAGCGATGGGCGAAAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGATTTTACCGCCAGCTGCCATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAACCCCACGGt  >  1:289577/1‑251 (MQ=255)
                                                                                                                                                                                                                                       tgccgatccccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACg    <  1:33455/241‑1 (MQ=255)
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                               
GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCG‑ACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGT  >  NC_000913/2173123‑2173603

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