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.96.1% 208.1 / ‑1.1 51pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base . (20/29);  minor base A (1/0);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-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.
*NC_0009132,173,3620C.96.1% 207.4 / ‑0.3 51pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base . (20/29);  minor base A (1/0);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-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.

GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173220‑2173497
                                                                                                                                             ||                                                                                                                                       
ggtctCCACGACGCACGGCGACCGCCCTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATCAGAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccacc                                                                                                                                         <  1:943924/139‑1 (MQ=255)
 aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATCGTGAGTGGGATAAAAATCAGGCTCGCCGATaccaccg                                                                                                                                        >  2:255377/1‑142 (MQ=255)
  tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAGCAGGCTTGCCGATcccaccgc                                                                                                                                       >  1:620226/1‑141 (MQ=255)
     ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAGATCAGGCTTGCCGATaccaccgccgt                                                                                                                                    <  2:724432/143‑6 (MQ=255)
      cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgta                                                                                                                                   >  1:1083796/1‑137 (MQ=255)
                  ccaCCGCCATCGCCACGAACAAGACGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCGCACAAACAGCAATTAACATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                     <  2:213697/141‑1 (MQ=255)
                  cGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                     <  2:517967/143‑1 (MQ=255)
                     ccGCCATCGCCACGAAGAAGCCGATGGTAGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                                  <  1:652441/143‑1 (MQ=255)
                      cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                                                                                                 >  1:9652/1‑143 (MQ=255)
                      cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                                                                                                 >  1:440624/1‑143 (MQ=255)
                       gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                                <  2:893536/143‑1 (MQ=255)
                        ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAACATCAGGCTTGCCGCTAC‑‑CCCCGCCGTACGTCCCAGCAGCACCGcc                                                                                                               >  2:42875/1‑143 (MQ=255)
                         cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                                                                                              <  2:528011/143‑1 (MQ=255)
                               cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACATCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                        <  2:821067/143‑1 (MQ=255)
                               cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCAACGCCGGATCaa                                                                                                        >  1:85027/1‑143 (MQ=255)
                                aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTGGCCGATACCA‑‑ACGCCGGATGGACCGGCAGCAACGCCGCATCaac                                                                                                       >  2:750375/1‑142 (MQ=255)
                                  gacgcggCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTGCCCTGCATACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                     <  2:108301/137‑1 (MQ=255)
                                  gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                     >  1:1103929/1‑143 (MQ=255)
                                    agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                   <  1:81400/143‑1 (MQ=255)
                                      aagCCGATGGTGGCAAGATCGCCAAACGGCCGCGCCGGATTACACGGCACACAAACAGCAAGTACACGGGGGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTCTGTCCCAGCAGCAACGCCGGGGCAAGGCAAAt                                                                                                 >  2:634411/1‑143 (MQ=255)
                                        cccGATGGCCGCAGGCTCGCCAAGCGGCAGCACCCGCTCCCCCGGCACACAAACAGCAATTAAAACGGTGAGTGGGATCAAAATCAGGCCTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                               <  2:172285/142‑1 (MQ=255)
                                          cgagtgtggCAGGATCTCCAAACGGCATCACCTCAGTACCCGGCACACAAGCAGCAAGTAAAATGGTGATTGCGATAAACAGCAGGCTTGCCGAGACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                             <  2:1022205/138‑1 (MQ=255)
                                          cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                             <  1:1008685/143‑1 (MQ=255)
                                          cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGCGAGTGGGATCAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCCCTCAgg                                                                                             >  2:221410/1‑143 (MQ=255)
                                                 ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTCGCCGATACCA‑‑CCGCCGTATGCCCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCt                                                                                      >  2:47458/1‑143 (MQ=255)
                                                    aaGGTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                   <  1:723864/143‑1 (MQ=255)
                                                    aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                   <  1:354786/143‑1 (MQ=255)
                                                    aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                   >  2:455662/1‑143 (MQ=255)
                                                    aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                   <  1:647321/143‑1 (MQ=255)
                                                    aaGATCGCCAAACGCCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                   <  2:315490/143‑1 (MQ=255)
                                                                cGGCAGCACCTGATTACCCGGCACGCAGCCAGCAATTAAAGTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                       <  2:666737/143‑1 (MQ=255)
                                                                 gtcagccCCTGATCACCCCGCACACGCACAGCCATTACAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      <  2:371490/141‑1 (MQ=255)
                                                                       aCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGGTGAAGATGGGGCATGCCGCTACAAAAAGTTCG‑‑GGTTGGAGCAGCAACGCCGGCTCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTac                                                                >  1:257468/1‑142 (MQ=255)
                                                                         cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                              <  1:272962/143‑1 (MQ=255)
                                                                                cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGCACACGACtccc                                                       >  1:735237/1‑140 (MQ=255)
                                                                                 ccGGCACACAAACAGCAATTAAAATGGTGACTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                      <  1:1017883/143‑1 (MQ=255)
                                                                                 ccGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGGCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAGACGACTACGa                                                      >  1:419537/1‑143 (MQ=255)
                                                                                  cGGCACACAAACAGCAGTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACCTCGCCTGTAATGGACTACGAg                                                     >  1:71737/1‑143 (MQ=255)
                                                                                  cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAGCTTCGCCTGTAAACGACTACGAg                                                     >  1:106395/1‑143 (MQ=255)
                                                                                            aCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACCTCGCCTGTAAACGACTACGAGCCTGCTTatc                                           >  2:572243/1‑141 (MQ=255)
                                                                                              aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACCAGCCTGCTTAGcca                                         >  2:502045/1‑141 (MQ=255)
                                                                                              aGCAATGAACATGGTGAGGGGGATAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAATGCCAATCAGGAACTCCTGTCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                         <  1:179694/143‑1 (MQ=255)
                                                                                               gCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGGACTCCTGACCGCCAAACTTCCCCTGTAGACGACGCAGAACGTGCTTAGGGAt                                        >  2:760674/1‑143 (MQ=255)
                                                                                                 aaTTAAAACGGCGAGGGGGATAAAAATCAGGCGTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCCGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                      <  1:422719/143‑1 (MQ=255)
                                                                                                  aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCACGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAACCTGCTTACCGAtcga                                     >  1:646590/1‑140 (MQ=255)
                                                                                                  aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGGAGGTCCGAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACGACGCGCCTGCTTAGCGATggg                                     >  2:1051984/1‑143 (MQ=255)
                                                                                                     aaaaCGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGTACTGCTGACCTCTTAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                  <  1:160592/143‑1 (MQ=255)
                                                                                                            tGAGTGGGATAAAAATCAGGCGTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcca                           <  2:356240/143‑1 (MQ=255)
                                                                                                                  ggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAGGGCCAGTCCGGGACTCCTGAGCGCCGAACGTCGCCTGTAACCGACTACGCGCCTGCTTAGCGAGGGGCGGTAAACCGTCCatg                     >  1:839163/1‑143 (MQ=255)
                                                                                                                     tAAAAATCAGGCTTGCCGATACCA‑‑CGGCCGTATGCCCCAGCAGCAGCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                  <  1:455662/143‑1 (MQ=255)
                                                                                                                      gaaacgcAGGCTTGGCGGTACCA‑‑CCGCCGCATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                 <  1:828530/137‑1 (MQ=255)
                                                                                                                      aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                 >  2:344482/1‑143 (MQ=255)
                                                                                                                        caaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt               <  2:858396/142‑1 (MQ=255)
                                                                                                                        aaaGCAGGCTTGCCGATACCA‑‑CCGCCGTACGTCCCAGCAGCAACGCCTGATCAAGGCCAATCAGGACCTGCGGACCGCCGAACGCCGCCTGTAAACTACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt               <  1:268488/143‑1 (MQ=255)
                                                                                                                            cAGGCTTGCCGATACCA‑‑CCGCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa           <  2:124106/143‑1 (MQ=255)
                                                                                                                                     cGATACCA‑‑CCGCCGTATGTCCCAGCATCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg  <  1:313143/143‑1 (MQ=255)
                                                                                                                                             ||                                                                                                                                       
GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173220‑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: