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.94.3% 144.2 / 0.7 35pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (15/18);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.3% 140.7 / 0.3 35pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (15/18);  minor base T (0/2);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173223‑2173497
                                                                                                                                           ||                                                                                                                                       
ttgcACGAGGCACGGCGACCGCCAGCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                                                                                                      <  2:284443/140‑4 (MQ=255)
ttCCACGATGCACGGCGACCGCCATCGCCACGAAGACGCCGATGGTGGCAAGATCGCCAAACGGCCGCACCGGATTACCCGGCACAAAAACAGCAATCAAAATGGTGAGCGGCAGAAAAACCAGCCTCACCGATaccgccgcc                                                                                                                                      >  2:808366/1‑140 (MQ=255)
ttCCACGATGCACGGCGACCGCCATCGACACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                                                                                                      <  2:657491/143‑4 (MQ=255)
   cACGATGCACTGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACAAACAGCAATTACAATGGTGAGAGGGAGAAAAATCAGGCGTGCCGATaccaccgccgga                                                                                                                                   >  2:754746/1‑137 (MQ=255)
          cccGGCGACGGCCATCGCCACGAAGAAGCGGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAACATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                          <  2:786111/141‑1 (MQ=255)
           aCGGCGACCGCCATCGCGACGAAGAGGCCGATGGGGGCAAGAGCGCCAAACGGCAGCACCTGATTACCCGGCACACCAACAGCAATTAAAACGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                         <  2:575422/143‑1 (MQ=255)
             ggCGACCGCCATCGCCACGAAGACGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                       <  1:763292/143‑1 (MQ=255)
             ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCCTCCCAACAGCAATTAAAATGGTGAGTGGGATAAAACTCAGGCTTGCCCATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                       >  1:565927/1‑143 (MQ=255)
             ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                       >  1:71925/1‑143 (MQ=255)
                     ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCCACGcc                                                                                                               >  1:158871/1‑143 (MQ=255)
                             accgagaagCCGATGGTGGCAACATCGCCAAACGGCAGCACGTGCGTCCCCGGCACACCAACGGCAATTAAAATGGTGGGTGGGATACACAGCCGGAAAGCCGTGAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                       <  1:937395/139‑1 (MQ=255)
                             aCGAAGAAGCCGCTGGTGGCAAGATCGCCGAGCGGCAGCACCGGATTACCCGGCACACAAACAGCAATTAACATGGTGAGTGGGAGAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                       <  1:487227/143‑1 (MQ=255)
                               gaagaaGCCGATGGCGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAACTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                     <  1:333235/143‑1 (MQ=255)
                                      ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                              <  2:807366/143‑1 (MQ=255)
                                      ccGAGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                              <  1:168156/143‑1 (MQ=255)
                                       cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                             <  1:827779/143‑1 (MQ=255)
                                         ctcgcggCACGCTCGCCAGTCGGCACCACCTGATTACCCGGCCCACAAACAGCAATTAAAATGGTCAGTGGGATGAAACTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGTATCAAGGCCAATCAGGaa                                                                                           <  1:404959/138‑1 (MQ=255)
                                         aTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                           <  1:203558/143‑1 (MQ=255)
                                         aTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAACGGCGAGCGCGTTACGTATCATGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                           >  1:458803/1‑143 (MQ=255)
                                               gCAAGATCGCCAAACGGAAGCACCTGATTACCCGGCACACGAACAACAAACCACAATGATATTTGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                     <  1:979187/143‑1 (MQ=255)
                                                    aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCTCACGGCCAATCCGGAACTCCTGAccgc                                                                                >  2:107792/1‑143 (MQ=255)
                                                           aaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                         <  1:655190/143‑1 (MQ=255)
                                                           aaCGCCCGCACCTCATCACCCGGCACACCACCAGCAATTAAAATGCTCACTGGGATAAAAATCAGGCTCGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                         <  1:174218/143‑1 (MQ=255)
                                                            aCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTAGCCGATTC‑‑CTGTGCAGTACAGCCCAACAGCAACGCCGGATCAAGGCCAATCAGGATCTCCTGACCGCCGAACTTc                                                                        <  2:700832/143‑1 (MQ=255)
                                                                      cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                              >  2:577861/1‑143 (MQ=255)
                                                                         aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACGTCGCCCGTAAACGAc                                                           >  2:1954/1‑143 (MQ=255)
                                                                              ccGGCACACAGAAAGCAGTTTAAATGGTGAGTGGGATAAAAATCAGGCTTGCCTATACAA‑‑CCGCCGTGTGTCCCAGCAGCAACACCGGATCAAGGCAAATCGGGAACTACTGACCACCGAAATACGCCTATAAAAGAATAcaa                                                      >  2:355323/1‑141 (MQ=255)
                                                                                   acacaAACAGCAATTAAAATGGTGAGTGGTATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCGGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTg                                                 <  2:972858/143‑1 (MQ=255)
                                                                                    cacgAACAGCAATTAAAATGGTGAGTGGGAGAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGc                                                <  2:278801/143‑1 (MQ=255)
                                                                                       aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCCACGCCGGATCAAGGCCCATCCGGAACTCCTGACCGCCGCACTTCGCCCGTAAACGGCCACCACCCTGccca                                             >  1:431663/1‑140 (MQ=255)
                                                                                                 ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                    >  1:809807/1‑143 (MQ=255)
                                                                                                  tAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCg                                   >  2:45384/1‑143 (MQ=255)
                                                                                                    aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                 >  1:503385/1‑143 (MQ=255)
                                                                                                       tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGCAATGACATCCTCCTTGCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGGTaaa                              >  1:455529/1‑143 (MQ=255)
                                                                                                            ggTGGGATAAAAATCAGGCTTGCCGGCACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                         <  1:447980/142‑1 (MQ=255)
                                                                                                                  aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCACTGCCGGACCAAGACCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                   <  1:454480/143‑1 (MQ=255)
                                                                                                                       aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt              <  1:731651/143‑1 (MQ=255)
                                                                                                                                  ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTCAATCACCCGt   >  2:478555/1‑143 (MQ=255)
                                                                                                                                   cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGACTGTAAACGACTACGAGCCTGCTTCGCGATGGGCGTTAAACCATCCATGATGGGGTTAACCATCCGTg  >  2:780573/1‑143 (MQ=255)
                                                                                                                                           ||                                                                                                                                       
TTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173223‑2173497

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