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.1% 148.7 / ‑2.4 35pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (22/12);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-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.97.1% 147.8 / ‑3.2 35pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (22/12);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 3.71e-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: 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.

ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173221‑2173496
                                                                                                                                            ||                                                                                                                                      
aTTTCCACGATGCACCGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                       <  1:435181/143‑2 (MQ=255)
    ccACGCTGCACGGCCACCGCCATCGCCACGAGGAAGCCGATGGTGCCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt                                                                                                                                   <  1:95398/143‑6 (MQ=255)
      aCGATGCACGGCGACCGCCACAGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                                 <  1:859733/143‑8 (MQ=255)
                  gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                                                                                                   >  1:819745/1‑143 (MQ=255)
                   cccGCCATCGCCACGAAGACGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAGATGGTGAGCGGGATGCAAATGAGGCTTGCCGATACAA‑‑CCGCCGTATGTCCCCGCCGCaa                                                                                                                  >  1:366236/2‑143 (MQ=255)
                     cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCCGCAGCAACg                                                                                                                >  1:732684/1‑143 (MQ=255)
                        cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCAGACAAACAGCAATTAAAATGGTGAGTGGGATACAAATCAGGCATGCCGATACAA‑‑CCGCCGGATGTACCAGCAGCAACGCCg                                                                                                             >  2:601134/1‑143 (MQ=255)
                          tCGCAACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                           <  2:651809/143‑1 (MQ=255)
                           cGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGCGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTAAGTCCCAGCAGCCACGCCGGAt                                                                                                          >  1:370955/1‑143 (MQ=255)
                               gcggtacagCCGATGGGCGAAAGGTCGCCAAACGTCAGCACCTGATTACCCGGCACCCGATCTGCAATTGAAATGGTGAGTGGGAGAATAATCTGGCTTGCCGATGCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                      <  2:789307/136‑1 (MQ=255)
                                cgcagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGCGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                     <  1:447815/140‑1 (MQ=255)
                                cGAACAACCCCAGGGTGGCAAGAGCTCAAAACAGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTTGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGTATCAAgg                                                                                                     <  2:660142/143‑1 (MQ=255)
                                             ggtggCAAGATCGCCACACGGCAGCGCTTGATTACCCGGCACAGAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                        <  2:761370/143‑1 (MQ=255)
                                                     gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGAccg                                                                                >  1:325087/1‑143 (MQ=255)
                                                      aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGAccgc                                                                               >  1:547315/1‑143 (MQ=255)
                                                       tCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCCACGCCCGATAAAGGCCAATCAGGAACTCCTGAccgcc                                                                              >  1:335202/1‑143 (MQ=255)
                                                          ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTACCACCCAGCAGCAACACCGGATCAAGGCCAATCAGGAACTCCAGACCGCCGaa                                                                           >  2:684789/1‑143 (MQ=255)
                                                              aCGGCAGCACCTTATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTTCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                       <  1:174708/143‑1 (MQ=255)
                                                                 gcagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGAGACAA‑‑CCGCCGTATGTCCCAGCAGCAACGCCAGATCAAGGACAATCAGGCACTCCCGACCGACGCACTGCGcc                                                                    >  2:2340/1‑143 (MQ=255)
                                                                      aCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                               >  1:27086/1‑143 (MQ=255)
                                                                        cTGATGACCCGGCACACGAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCGTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                             <  2:699580/143‑1 (MQ=255)
                                                                             tACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                        >  2:284937/1‑143 (MQ=255)
                                                                                ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACAGCCGAACTGCGCCTGTAAACGACTCCGa                                                     >  2:391508/1‑143 (MQ=255)
                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTg                                                <  1:672621/143‑1 (MQ=255)
                                                                                             aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTCAACGACTACGAGCCTGCTCAGCGa                                        >  1:99104/1‑143 (MQ=255)
                                                                                                  ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                   <  1:188781/143‑1 (MQ=255)
                                                                                                  ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTAGGTCGCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGCTGGGc                                   >  2:822824/1‑143 (MQ=255)
                                                                                                             aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGTCGTTCAAccatc                        >  2:809095/1‑143 (MQ=255)
                                                                                                                    tAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCGAAGCAGCAACGCCGGATCACGGCAGATCAGGAGCTCCTGCCCGCCGAACTTCGCGTGTAAACGACTACAAGCCTGCTTAGCGCTGGGCGCTACAccatccaacatg                 >  1:819743/1‑138 (MQ=255)
                                                                                                                    tAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                 <  1:91235/143‑1 (MQ=255)
                                                                                                                     aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACAGCAGAACTTCGCCTGTAAACGACTACGAGCCTACTTAGCCATGGGCAGTAAAGAATCAatgatga                >  2:588102/1‑142 (MQ=255)
                                                                                                                         aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTAGCCTGTAAACGACTACGAGCCTGCTTCGCGATGGGCGGTAACCCATCCAAGACGGGttt            >  2:447526/1‑143 (MQ=255)
                                                                                                                           cAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa          <  1:308130/143‑1 (MQ=255)
                                                                                                                            aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAGCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt         <  2:819745/143‑1 (MQ=255)
                                                                                                                              gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATCGGCGTTAAACCATCCATGATGGCTTCAATCa       >  2:903123/1‑143 (MQ=255)
                                                                                                                                 tGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc    <  1:284937/143‑1 (MQ=255)
                                                                                                                                   ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt  >  2:247037/1‑143 (MQ=255)
                                                                                                                                   ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAAGCCAATCAGGAACTCCTGACCGCCCAACTTCGCCTGTAAACGACTACAAACCTGCTTACCGATGAGCGTTAGACCATCCATAATGGAGTTAATCACCCGt  >  2:755203/1‑143 (MQ=255)
                                                                                                                                            ||                                                                                                                                      
ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173221‑2173496

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: