Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,361 Δ2 bp 100% intergenic (‑488/+920) 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,3610C.94.7% 73.6 / ‑1.0 19intergenic (‑488/+921)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (10/8);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*NC_0009132,173,3620C.94.7% 71.1 / ‑1.0 19intergenic (‑489/+920)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (10/8);  minor base A (1/0);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCC  >  NC_000913/2173093‑2173585
                                                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                               
gCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGCTGATGCTCATAATGATGACACCCGAGCTTAACGTGCGGAACAGCTTTCCACGATGCACGCCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTaccaccgccga                                                                                                                                                                                                                            >  1:50918/1‑270 (MQ=255)
                                                             tAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATg                                                                                                                                                                                                                       >  1:214216/1‑217 (MQ=255)
                                                                               acaGGGTATAGCTCAAAAGGTGGAAACCCAGGATTAAGGAGCGGAACAGATTTCCACGCTGCACGGCGACGGCCACCGCCACGAAGACGCCGCAGGTGGCAGGATAGCCAAACTGCATCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGCGGGATAAAAATCAGGCGTGCCGATACCA‑‑CCGCCGTATGTCCCATCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                            <  2:50918/274‑1 (MQ=255)
                                                                                        tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGACGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAACAATCAGTCTTGCCGATaccacc                                                                                                                                                                                                                                 >  1:165422/1‑182 (MQ=255)
                                                                                        tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCATGACCGCCGACCTTCGCCTg                                                                                                                                                           >  1:185814/1‑250 (MQ=255)
                                                                                        tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                  >  1:125940/1‑275 (MQ=255)
                                                                                           tCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGCCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                <  1:284594/274‑1 (MQ=255)
                                                                                                   atgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCCGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGCCCACGGACCTGCTGAGCGAggggggggta                                                                                                                       >  2:284593/1‑269 (MQ=255)
                                                                                                         cccGAGATTACGGTGCGGAACAGATTTCCCGGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGAACCTGATTACCCGCCACACAAACAGCAATTAAAATGGTGAGTGAGATAAAAATCAGGCTTGCCGATaccaccgccgta                                                                                                                                                                                                                           <  1:241974/171‑7 (MQ=255)
                                                                                                                       ggggAACAGATTTCCACAATGCACGGCGACCGCCATCGCCCCGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCAAACAAACAGCAATTAAAATTGTGAGTGGGATAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                                                                                                   <  1:348281/273‑1 (MQ=255)
                                                                                                                                    ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                                                                                                                     <  1:346506/244‑1 (MQ=255)
                                                                                                                                                                     aagCCGATGGGGGAAAGACGGCCAAACTGAAGAACCTGATTACCCGGAAAACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACTA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGGTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTGTTTACCGCCAg                                                     <  2:125940/275‑1 (MQ=255)
                                                                                                                                                                                                                  ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                                                                                                >  1:365347/1‑139 (MQ=255)
                                                                                                                                                                                                                  ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGACGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAACCGCGAGGATGCCGATAAt        >  1:20988/1‑275 (MQ=255)
                                                                                                                                                                                                                      cacaAACATCAAATTAAATGGTGAGTGGGATAAAAATCAGGGGTGCCGAGACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAGCAAGGGCGATCAGGAACTCCTGAGCGGCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                            <  1:1534/171‑1 (MQ=255)
                                                                                                                                                                                                                      cacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                            <  2:365346/139‑1 (MQ=255)
                                                                                                                                                                                                                         aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCAGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTCATCACCCGTGGGATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGACGCAATACACATTTCACATCGTAACCCGCCAGCATGCCGATAACCAGccc  >  1:80070/1‑274 (MQ=255)
                                                                                                                                                                                                                                     aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGCCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTACACGACTACGAGCCTGCTTAGCGATGGTGGTTAAACCATCCCTGATGGGTTTAATc                                                                                                 >  1:241918/1‑167 (MQ=255)
                                                                                                                                                                                                                                                  gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc                                                                                                                                                                                  >  2:348846/1‑73 (MQ=255)
                                                                                                                                                                                                                                                      aaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGATCTCCt                                                                                                                                                                              <  1:348847/73‑1 (MQ=255)
                                                                                                                                                                                                                                                            aGGCTTTCCGATA‑‑ACACCGCCGTATGTCCCAGCAGCAATGTCTGATCAATGCCTATCAGGAACTTCTATCTGCCGAACTTCCCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                                                                                                         >  2:321677/1‑120 (MQ=255)
                                                                                                                                                                                                                                                                 tGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTag                                                                                   >  1:129568/1‑153 (MQ=255)
                                                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                               
GCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCC  >  NC_000913/2173093‑2173585

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: