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.80.0% 68.7 / 11.9 25pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base . (15/5);  total (17/8)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*NC_0009132,173,3620C.80.1% 69.9 / 8.2 25pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (15/5);  total (17/8)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.

ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAG  >  NC_000913/2173221‑2173504
                                                                                                                                            ||                                                                                                                                              
aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                               >  2:440430/1‑142 (MQ=255)
     cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGACGGTGGCAAGATCGCCAGACGCCAGCACCTGAGTACCCGGCACACAAACAGCACTCAAAATGGTGCCGGGCAGAAAACTCAGGCTCACCGACaccaccgccgga                                                                                                                                          >  2:457858/1‑137 (MQ=255)
      aCGATGCACGGCGACCGCCATCCCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                                         <  1:60206/143‑8 (MQ=255)
             aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                                >  2:66357/1‑143 (MQ=255)
                 cGACCGCCATCGTCATGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                            <  2:481587/143‑1 (MQ=255)
                         aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGTGGGCTAAAAATCAGGCTCGCCGATACCCCCCCCGTAGGTCCCAGCCGCGACCCCgg                                                                                                                    >  2:561661/1‑143 (MQ=255)
                                cGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTAGAAATTA‑‑TCTCCTCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                             >  1:339716/1‑143 (MQ=255)
                                                ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                             >  1:237855/1‑143 (MQ=255)
                                                  cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                                           <  2:80395/143‑1 (MQ=255)
                                                   aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                          >  2:202676/1‑143 (MQ=255)
                                                          gcAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                                   <  2:80511/142‑1 (MQ=255)
                                                          ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGCGTGGGATACAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                                   >  1:514265/1‑143 (MQ=255)
                                                           cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCGGGAACTCCTGACCGCCGacc                                                                                  >  1:232795/1‑141 (MQ=255)
                                                             aaCGGCAGCACCTGATGACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGACCAATCAGGAACTCCTGACCGCCGAACtt                                                                                <  1:544547/143‑1 (MQ=255)
                                                              aCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGATTGCCGATTCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                               >  1:305366/1‑143 (MQ=255)
                                                                  cagcaACTGATCACGCGGCGCACAAACAGCGTTCAAAAGGGTGCGTGGGAGAAAAACCAGCCTTGCCGACAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGGCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                           <  2:515339/143‑1 (MQ=255)
                                                                           aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGAACTTCGCCTGCAAACGAc                                                                  >  2:632139/1‑143 (MQ=255)
                                                                                 cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                            >  1:707408/1‑143 (MQ=255)
                                                                                    cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGGATGTCCAAGCAGCATCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                         >  1:597217/1‑143 (MQ=255)
                                                                                         aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGGCCGCCAAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                    <  1:440430/143‑1 (MQ=255)
                                                                                                  ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCCTCGCGATGGGc                                           >  1:522498/1‑143 (MQ=255)
                                                                                                         ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAGCGAAGACGAGCCTGCGCAGCGAAGGGCGCTCAAc                                    >  2:252278/1‑143 (MQ=255)
                                                                                                          gTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGCTAAAcc                                   >  1:107339/1‑143 (MQ=255)
                                                                                                          gTGAGTGGGATAAAAATCAGGCTCGCCGATACCA‑‑CCGCCGGATGTCACAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAATTCGCCTGTAAACGACTACGAGCCTGCTTAGCGCAGGGCGCTAAAcc                                   >  2:297543/1‑143 (MQ=255)
                                                                                                           tGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAACcca                                  >  1:19630/1‑143 (MQ=255)
                                                                                                             aGTGGGATAAAAATGATGCTTGCCG‑‑ACCCCCCCTACGTATGTGCCAGCAGCAACGCCGGATCAAGGCCAATCTGGTTCTCCTGACCGCCGATCGTCGCCTGTAACCTACTTCGAGCCTGCTTAGCGATGGGCGTTAAAccatc                                <  2:18211/143‑1 (MQ=255)
                                                                                                                        aaTCAGGCTTGCCGCTAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCCGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt                     <  2:194641/143‑1 (MQ=255)
                                                                                                                                   ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCCACTACGAGCCTGCATCGCGAAGGTCGTTAAACCATCCATGATGGGGTTAATCACCCGt          >  2:240480/1‑143 (MQ=255)
                                                                                                                                             cccGCCGGATGTCCCCGCAGCCACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTag  <  2:305366/143‑1 (MQ=255)
                                                                                                                                            ||                                                                                                                                              
ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAG  >  NC_000913/2173221‑2173504

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: