Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp pseudogene (915‑916/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.92.9% 101.1 / 4.2 29pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (13/13);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.92.9% 100.5 / 1.6 29pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base . (13/13);  minor base T (0/2);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCA  >  NC_000913/2173220‑2173491
                                                                                                                                             ||                                                                                                                                 
gATTTCGACACTGCTCGTGGACCGCCATCGCCACGAAGAAGCCGTTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGGTTGCGGATaccacc                                                                                                                                   <  2:66611/143‑1 (MQ=255)
       aCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                            <  1:229733/143‑8 (MQ=255)
           tgccgggcggcCGCGATGGCTACGAAGAAGCCGGTGGTGGCAAGATCGGCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcc                                                                                                                      <  1:818648/138‑1 (MQ=255)
           tGCGCGGCGACCGCGATCGCCAGGAAGAGGCGGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcc                                                                                                                      <  2:706276/143‑1 (MQ=255)
                ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                 >  1:103525/1‑143 (MQ=255)
                 gCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                >  2:662112/1‑143 (MQ=255)
                 gCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGGCCcagcag                                                                                                                >  2:536096/1‑143 (MQ=255)
                       gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                          <  2:391428/143‑1 (MQ=255)
                        ccATCGTCACGAAGAAGCTGATGGTGGCAAGATCGTCAAATGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                         <  1:545382/143‑1 (MQ=255)
                          gtCGCGACGAAGAAGCCGGTGGTGGCAAGAGCAACATAATGGTTTACGTTATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                       <  2:369026/142‑1 (MQ=255)
                              ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAACATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCTcc                                                                                                   >  1:712570/1‑142 (MQ=255)
                                                   cAAGATCGCCAAGTGGCAGCACCTGATTACCCGGCACACATACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                              <  2:615616/143‑1 (MQ=255)
                                                    aaGATCGCGAAACGGCAGCACCTGATTACCCGGCCCACAAACAGCAATTAAAATGGTGCGTGGGATAAAACTCAGGCTTGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                             <  2:281407/143‑1 (MQ=255)
                                                            cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAATTTGTGATTGGGATAAAAATCATGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                     <  1:41026/143‑1 (MQ=255)
                                                                     gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGAGATAAAAATCAGGCTTGCCGATACCA‑‑CGGCCGTATGTCCCAGCAGCAACTCCGCATCAAGGCGAAACAGGAACTCCTAACCGCCGAACTTCGCCtca                                                            >  2:13442/1‑141 (MQ=255)
                                                                     gaaCCTGATTATCCGGCACACAAACAGCAATGAAAATGGTGAAAACCACTGGTAACGGGTTTGCCGATTTCATCTCCCTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                            <  2:600499/141‑1 (MQ=255)
                                                                           gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                      <  2:103525/143‑1 (MQ=255)
                                                                                 ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGCACTCCTGCCCGCCGAACTTCGCCTGTAAACAACTAcaa                                                >  1:250227/1‑141 (MQ=255)
                                                                                                cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTAACCGCCCAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGaag                                 >  2:16286/1‑141 (MQ=255)
                                                                                                 aaTTGAAAGGGGGAGTGGGATAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCAGACATCCGATCTTCTCTTGTAAACGACTACGAGCCTGCTTAGCGATgg                                <  1:158893/143‑1 (MQ=255)
                                                                                                      aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTATAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                           >  1:457940/1‑143 (MQ=255)
                                                                                                       aaTGGTGAGTGGGATAAAAATCAGGCTTGACGATACCA‑‑CCGTCGTATGTCCCAGAAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                          <  1:690900/143‑1 (MQ=255)
                                                                                                              agggggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                   <  2:599355/140‑1 (MQ=255)
                                                                                                              aGTGGGATAACAATCAGGCTTGCTGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                   <  1:782316/143‑1 (MQ=255)
                                                                                                              aGTGGGATAAAAATCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                   <  1:730442/143‑1 (MQ=255)
                                                                                                                tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcca                 >  2:691455/1‑143 (MQ=255)
                                                                                                                    aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCCCCGCCGGATCACGGCCCATCAGCCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTAAGCGATGGGCGTTAAACCATACatgaa             >  1:481597/1‑142 (MQ=255)
                                                                                                                        ccagcAGGCTTGTCGATCCCA‑‑CCGCCGGACGTGCCAGCAGCCTCGCCGGATCATGGCCGAGCAGGAACTCCTTGCCGCCGATCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt         <  2:464978/139‑1 (MQ=255)
                                                                                                                        aaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCCGGAACTCCTGCCCGCCAAACTTCGCCTGTAAACGACTACGAGCCTGCTAAGCGATGAGCGTTAAACCATCCATGATGGGt         >  1:315504/1‑143 (MQ=255)
                                                                                                                              ggCTTGCCGATACCA‑‑CCGCCATATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTCAACGACTACGAGCCTGCTTAGCGATGCGCGTTAAACCATCCATGATGGGCTTAATc   >  2:314668/1‑143 (MQ=255)
                                                                                                                               gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGCTGGGCGTTAAACCATCCATGATGGGatcaatca  >  2:156299/1‑143 (MQ=255)
                                                                                                                                             ||                                                                                                                                 
GATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCA  >  NC_000913/2173220‑2173491

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: