Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,361 Δ2 bp 100% pseudogene (917‑918/1358 nt) gatC ← galactitol‑specific PTS enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.97.2% 152.5 / ‑2.5 36pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (16/19);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
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.100.0% 155.2 / NA 35pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (16/19);  total (16/19)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173219‑2173497
                                                                                                                                              ||                                                                                                                                       
aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCACTTCAAATGGTGAGTGGGATAAAAAAAAGGCTTGACAATaccac                                                                                                                                          >  2:457816/1‑143 (MQ=255)
    ttCCACGATGCACGGCGACCGCCATGGCCCCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACAACCAGCAATTAAAATGGAGAGAGGGATAAAAATCAGGCTTGCCGGTAGCACcgcg                                                                                                                                      >  1:845294/1‑143 (MQ=255)
      ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt                                                                                                                                    >  1:643611/1‑138 (MQ=255)
       cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCAAACAAACAGCAATGAAAATGGTGAGAGGGATAAAAATCAGGCTTGCCAATgccacggccgta                                                                                                                                   >  2:545838/1‑136 (MQ=255)
       cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGAAGCACCGGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgta                                                                                                                                   >  1:493155/1‑137 (MQ=255)
          gATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                              <  1:370950/143‑1 (MQ=255)
          gATCCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGAGTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                              <  2:563143/143‑1 (MQ=255)
            tGCACGGCGACCGCCAACGCCACGAAGAAGGCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACAGCAATTAAAATGGTGCGTGGGATAACAATCAGGCTTGCCGATACCA‑‑CCGCCGCATGTcc                                                                                                                            <  2:153062/143‑1 (MQ=255)
                  tcgcccGCCATCCCCTCCAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                      <  2:201612/139‑1 (MQ=255)
                       cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCGGTATGTCCCAGCAGCAACg                                                                                                                 >  2:490420/1‑143 (MQ=255)
                         ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTAGCCGGCACACAACCAGCAATTAAACACGTGCGCGGGCTAAAAATCAGACCTGCCGATACCA‑‑CCGCCGTCCGCCCCAGCAACAACAcc                                                                                                               >  2:577421/1‑143 (MQ=255)
                          cacggCCACGAAAAAGCCGATGCTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                                                                                              <  2:120792/139‑1 (MQ=255)
                              gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATGAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACGA‑‑CCGCCGTATATCCCAGCAGCAGCGCGGGATc                                                                                                          >  1:650926/1‑143 (MQ=255)
                                cagcgataagCCGAGGGTGGCAAGATCGCCAAACTGCAGCTCCTGATTACCCGGCACACAAACAGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                        <  2:541266/136‑1 (MQ=255)
                                   gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTCCCGCTACCA‑‑CCGCCGTATGTCTCGGCAGCCTCCCCGGACCACGGc                                                                                                     >  1:323371/1‑143 (MQ=255)
                                          ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCACCGCCGGACCAAGGCAAATCAg                                                                                              >  2:752112/1‑143 (MQ=255)
                                                 tggCAAGATCGCCAAACGGCCGCACCTGATTCCCCGGCACACAGATAGCAATTAACATGGTGAGTGGGACAAAAATTAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc                                                                                       <  2:220465/143‑1 (MQ=255)
                                                 tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCACCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc                                                                                       <  1:841463/143‑1 (MQ=255)
                                                    cAAGACCGCCAAACGGCAGCACCTCATTACCCGGCACACCAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGCCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                                    <  1:8964/143‑1 (MQ=255)
                                                        cttgCCAAACGCGAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTTCTGAccgc                                                                                <  2:697666/140‑1 (MQ=255)
                                                        aTCGCCAAACGGCAGCACCTGATTACCCGGCTCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCCAGCAGGAACCCccgcccgc                                                                                >  2:824669/1‑143 (MQ=255)
                                                            ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                            <  1:828211/143‑1 (MQ=255)
                                                            ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                            <  1:828210/143‑1 (MQ=255)
                                                             cAAACGGCAGCACCTGATGACCCGGCACACAAACAGCAGTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                           <  2:241690/143‑1 (MQ=255)
                                                                  ggCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGCCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      <  1:442748/143‑1 (MQ=255)
                                                                  gcctgcaGCTGTTGACCCTGCCCCCACACAGCAATTAGAATGGTGAGTGGGATAAAGATCAGGCTTGCCAATACCA‑‑TCGCCGGATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      <  2:632393/139‑1 (MQ=255)
                                                                          cTGATTACCCGGCACACAAACAGCAATTAAAAAGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCTCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCAAATAAGGAACGCCTGACCGCCGCACTTCGCCTGTAAAc                                                              >  1:660481/1‑143 (MQ=255)
                                                                                aCCCGGCATACAAACAGAAATTAAAATGGTGAGTGGGATAAAAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                        <  1:520251/143‑1 (MQ=255)
                                                                                   cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGACGATCAGAAACTCCTGACCGCAGAACTTCGCCTGTAAACGACTACGAg                                                     >  2:758880/1‑143 (MQ=255)
                                                                                      cacccaAACAGCAATGAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                  <  1:357865/143‑1 (MQ=255)
                                                                                         acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGGCAATCAGGAACTCCTAACCGCCGAACTTCGCCTGTAAACAACTACGAGCCTGCt                                               >  1:307381/1‑143 (MQ=255)
                                                                                             ccAGCAAATAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                           <  1:214483/142‑1 (MQ=255)
                                                                                                   aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGAGCCTGCCTAGCGATggg                                     >  2:88771/1‑143 (MQ=255)
                                                                                                       aaaTGGTGAGTGGGAAAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCCTGCTCAGCGATGGGCGtt                                 >  2:744059/1‑143 (MQ=255)
                                                                                                          tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGAGAGTGCTTAGCGATGGGCGTTaaa                              >  2:139360/1‑143 (MQ=255)
                                                                                                               aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                         >  1:865107/1‑143 (MQ=255)
                                                                                                                     aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTACTGACCGCCGAACTTCGCCGGTAAAGGACTACGAGCCTGCTTAGCGATGGGGGTTAAAACATGCatgat                   >  2:642581/1‑143 (MQ=255)
                                                                                                                           aTCAGGCGTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt             <  1:337277/143‑1 (MQ=255)
                                                                                                                                gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCCTGCTTCGCGGTGGGCGTTAGACCCTCCAGTCTGGGTTTAATCa        >  1:396212/1‑143 (MQ=255)
                                                                                                                                      cGATACCA‑‑CCGCCGTAGGGCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAACTCCTGACCGCCGAACCGCGCCTGTAAACGACTACGAGCCTGCTCAGCGAGGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg  >  1:768468/1‑143 (MQ=255)
                                                                                                                                              ||                                                                                                                                       
AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173219‑2173497

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: