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.90.9% 80.0 / 1.5 22intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (17/3);  total (18/4)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
*NC_0009132,173,3640C.90.9% 81.0 / 0.4 22intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (17/3);  total (18/4)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173122‑2173612
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                         
aGCATTAGCCGCCAGCTGGGAGCGCAAACCACTAGCTTGCGTCGCGATCCCCAGGGTGATGCTCATAATGATGACACCCGAGCTTACGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTc                                                                                                                                                                                                                                                 >  2:457683/1‑251 (MQ=255)
         cgccAGTTGGGTGTGCAAACCAATAGGTTGGGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAgca‑gca                                                                                                                                                                                                                                       <  1:277749/251‑1 (MQ=255)
                 gggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGAGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGTTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCgg                                                                                                                                                                                                                                >  1:489328/1‑250 (MQ=255)
                   gtgtGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATc                                                                                                                                                                                                                             >  1:149953/1‑251 (MQ=255)
                         aaaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGcc                                                                                                                                                                                                                       >  1:340364/1‑251 (MQ=255)
                                                               cATAATGATGCCACCGGAGATTAGGGTGCGGAACAGATGTCCGCGTAGCCCGGCGCCGGCCATCGCCCCAAAGAAGCCGATGGTGGCAAGAGCGCCAACCGGCAGCACATGATTGCCCGCCACACAAACAAAAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGAGAC‑‑CTCCGCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                                                 <  1:131240/251‑1 (MQ=255)
                                                                aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTCGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACg                                                                                                                                                                                >  2:117581/1‑251 (MQ=255)
                                                                       tgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                          >  2:169650/1‑250 (MQ=255)
                                                                                     aaGGTGCGAAACAGATTCCCACGAGGCACGGCGACCGCCACCGCCCGGAAGAAGCGGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCGGTAAGCGACTACGAGCCTGCTTAGCGAt                                                                                                                                                           <  1:457683/251‑1 (MQ=255)
                                                                                                                     gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc                                                                                                                            >  2:423046/1‑250 (MQ=255)
                                                                                                                               gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAt                                                                                                                 >  1:114570/1‑251 (MQ=255)
                                                                                                                                                            cccAAACGGCACCCCCTGATTCCCCGGCACACAAACAGCAATAAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCAATCCCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGTGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCg                                                                                    <  1:463213/250‑1 (MQ=255)
                                                                                                                                                               aaaCGGCAGCACCTGGTTACCCGGCCCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                  >  1:136428/1‑250 (MQ=255)
                                                                                                                                                                              aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCtt                                                                  >  2:141348/1‑251 (MQ=255)
                                                                                                                                                                                    cGGCACACAAACAGCAGTTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGCCCAATCCGGAACTCCTGACCGCGGACCTTCGCCTGTAAACGACTACGACCCTGCTTCGCGCTGGGCTTTAAACAATCCATGATGGGTTTAATCACCCGTGGCATTCGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATCCACCTTTGACac                                                            >  2:40695/1‑250 (MQ=255)
                                                                                                                                                                                            aaaCAGCAATTAAAATGGTGGGTGGGATAAAAATCAGGCTTGCCGATA‑‑ACACCGCCGTATGTCCCAGCA‑GCAACGCAGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCAGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAccc                                                    >  2:37838/1‑251 (MQ=255)
                                                                                                                                                                                                aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCGATACACCTTTGACATCGTAACCcgcg                                                 >  1:492880/1‑250 (MQ=255)
                                                                                                                                                                                                       aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                         >  1:298337/1‑251 (MQ=255)
                                                                                                                                                                                                           tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                      >  1:276176/1‑250 (MQ=255)
                                                                                                                                                                                                            ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATa                                     >  2:161305/1‑250 (MQ=255)
                                                                                                                                                                                                              tGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATAGTCACCCGCGAGGATGCCGATAATc                                  >  1:264478/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaatccagcatttccgatcccgccgccgTATTTCCCAGAACGCA‑CCCCGGAACAAGGCCAACCGGGACCGCCTGACCGCCAAACTCCGCCTGTAAACGACTCCGAGCCTGCTGAGCGAGGGGCGTCAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCTGCAGCTCTGCCGCCGTTTTTTCCGCCAGCTGCAATACTCCTTTGACACCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac                       <  1:40695/229‑1 (MQ=255)
                                                                                                                                                                                                                                ggCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTCAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAccc                >  2:466440/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gataccaccgccgTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCa                                                                                                      >  1:304975/8‑157 (MQ=255)
                                                                                                                                                                                                                                             accaccgccgTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAgg     >  2:197190/5‑251 (MQ=255)
                                                                                                                                                                                                                                                cccGCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc  <  1:197190/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                         
AGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCA‑GCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173122‑2173612

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