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.96.3% 110.4 / ‑2.6 27intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (17/9);  total (17/10)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01
*NC_0009132,173,3640C.96.3% 110.3 / ‑1.5 27intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (17/9);  total (17/10)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

GCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC  >  NC_000913/2173117‑2173600
                                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                            
gCCCCCGCATTAGCCGCGCGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGCTACACAGGGCGATGCTCATAATGATGACACCCGAGATTAAGATGCGCAACAGATTTCCACGACGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGACACACGGCAGCACCTGATTACCCGGCACACACACAGCAATTCAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACcccccccgt                                                                                                                                                                                                                                           >  2:386804/1‑248 (MQ=255)
              cgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGAGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCagcagc                                                                                                                                                                                                                            >  1:402038/1‑250 (MQ=255)
                ccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCTATCCACAGGGTGATGCTCATAATGATGACCCCCGAGATTAAGGTGTGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACAAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGAAGCACCTGATTACCCGGCACACAATCAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAAc                                                                                                                                                                                                                         <  1:497021/251‑1 (MQ=255)
                  aGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGCTCCCCCGGGTGATCCTCATAATCATGACACCCGAGCTTAAGGTGCGAAACAGATTACCCTGATGCACGGCCACCGCCATCGCCACGACGAAGCCGATGGAGGCAAGATCGCAACACGGCAGCACCTGATTACCCGGCAAACACACAGCAATTAAAATGGAGAGTGGGATAAAACTCAGGCTTGCCGATACCACA‑‑GCCGTATGTACCAGCAGCAACGc                                                                                                                                                                                                                       >  1:286795/1‑251 (MQ=255)
                                               cgcgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                                                                                                                          >  2:131472/1‑251 (MQ=255)
                                                       acaGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                                                                                                                                  <  1:131472/251‑1 (MQ=255)
                                                        caGGGTGATGCTCATAATGATGAACCCCGAGATTAAGGTGCAGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGAGGCCAGATCGCCAAACGGCAGCACCTGATTACCCGGCCCACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGACGATACCACC‑‑GCCGTATGTCCCAGCAGCCACGCCGGATCAAg                                                                                                                                                                                                              >  2:508632/1‑222 (MQ=255)
                                                                      aattgatgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                                                                                                                                   <  1:143545/248‑1 (MQ=255)
                                                                              acatccGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGCAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATAAGGCGTGCCGATACCACC‑‑GCCGTATGTACCAGCAGAAACGCCGGATAACGGCAAATCAGGCACTCCTGACCCCAGAACTTCCCCTGAAAACGACTACGacc                                                                                                                                                           >  2:113965/5‑249 (MQ=255)
                                                                               cacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGAGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGcc                                                                                                                                                          >  1:45986/1‑251 (MQ=255)
                                                                                ccccGAGATTAAGGGGCGGAAAAGATTTCCCCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTCCCCGGCACCCAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                                                                                                                         <  1:281911/250‑1 (MQ=255)
                                                                                     agtttAAGGGGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                 <  1:434158/235‑1 (MQ=255)
                                                                                     agaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                 >  2:434158/1‑238 (MQ=255)
                                                                                      gaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTCCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGACCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                                     >  2:420336/1‑249 (MQ=255)
                                                                                                       gATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCCGATCAAGGCCAATCAGGAACTCCTGGCCGCCCACCTTCGCCTGGAAACCACTACGAGCCTGCTTAGctatcgccgtacaaac                                                                                                                                   >  2:137300/1‑235 (MQ=255)
                                                                                                                         cccccccacccGCCACGAAGAAGCGGAGGGTGGAAAGACCGCAAAACGGACGCCCCGAATTACCCGGCACACAACAAGAAAGAAAAAGGGGGAGTGGGAAAAAAACCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc                                                                                                                <  2:580087/241‑1 (MQ=255)
                                                                                                                                  tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGAAAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAACCCATCCATGATgg                                                                                                                        >  2:270614/1‑234 (MQ=255)
                                                                                                                                  ccGCCACGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCGGCACCTGATTACCCGGCACAAAAACAGCAATTAAAATGGTGAGTGGGAAAAAAATCAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                                                                                                                        <  1:270614/233‑1 (MQ=255)
                                                                                                                                                              aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                                                                                                                                      >  1:5479/1‑176 (MQ=255)
                                                                                                                                                                        ggCAGCCCCTGATTACCCGGCACACAAAAAGCAATAAAAATGGTGAGGGGGATAAAAATCAGCCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACTATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                                 <  1:211362/251‑1 (MQ=255)
                                                                                                                                                                                       cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtttt                                                                             >  2:114113/1‑224 (MQ=255)
                                                                                                                                                                                       cccGGCACACAAACAGCAATAAAAATGGTGAGTGGGATAAAAATCAGGCTTCCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtttt                                                                             <  1:114113/224‑1 (MQ=255)
                                                                                                                                                                                            cacccaAACAGAAATTAAAAGGGTGAGGGGAAAAAAAATCAGGCGTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                      <  1:103817/222‑1 (MQ=255)
                                                                                                                                                                                            cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                      >  2:103817/1‑226 (MQ=255)
                                                                                                                                                                                               acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAc                                          >  2:62466/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGCGGATGCCGATAATCAGCCCCATCACAAAAcccc   >  2:449578/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTaa                                           >  1:276741/1‑210 (MQ=255)
                                                                                                                                                                                                                                       cTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAc  >  2:394691/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                            
GCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC  >  NC_000913/2173117‑2173600

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: ‑