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.90.0% 70.6 / 0.8 20pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (9/9);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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.90.0% 68.9 / 0.4 20pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (9/9);  minor base T (1/1);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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.

TCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATC  >  NC_000913/2173224‑2173490
                                                                                                                                         ||                                                                                                                                
tCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCATGATCGCCAAACGGCGGCACATGATTACACGGCACACAAACACCAATTAAAATGGTGAGGGCGATAAAAACCAGCCTTGCAGCCaccaccgccg                                                                                                                              >  2:256402/1‑139 (MQ=255)
 ccATGATGAACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGTCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt                                                                                                                             <  2:73154/143‑6 (MQ=255)
   acgacgCACGGCGCCCGCCATCGCCACGAACAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                           <  1:100009/143‑8 (MQ=255)
     gATGCACGGCGACCGCCATCGCCACGAAGCAGCGGCGCGAGTCGAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                       <  2:593857/143‑1 (MQ=255)
                 cctgcctCGCCACGTAGATGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                           <  2:607382/137‑1 (MQ=255)
                        ggCCACGAAGCAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                    <  2:82141/142‑1 (MQ=255)
                        cGCCACGAGGAAGCCGATGGTGGCAAGATCGCCCAGCGGCAGCACCTGATTACCCGGCACACAAACAGCAATGAACATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                    <  2:775365/143‑1 (MQ=255)
                         gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTCATTACCCGGCACACAAACCGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                   <  1:289764/143‑1 (MQ=255)
                                agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAAGCCa                                                                                            >  2:800800/1‑143 (MQ=255)
                                    gCCGATGATGGCAAGATCGCCAAACGGCAGCACCTGATTACACGGCACACAAACAGCAATCAAAATAGTGAGCAGGATAAAAATCAGGCTAGCCGCTACCA‑‑CCGACGAATGGCCCAGAAGCCCCGCCGGATACACGCCAACCa                                                                                        >  2:328432/1‑143 (MQ=255)
                                      cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGAGAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAgg                                                                                      >  2:527371/1‑143 (MQ=255)
                                           gtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACCGCACTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                 >  1:461497/1‑143 (MQ=255)
                                                   accGCCAAACGGCAGCACCTGATGACCCGGCACACAAACAGCAATGAAAATGGTGAGTGGGATAAAAATCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                         <  2:744409/141‑1 (MQ=255)
                                                             ggCAGCACCTGATTACCCGGCACACAAGCAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCACCGCCGGATCAAGGCCAATCAACACCTCCTGACCGCCGAACTTCGc                                                               <  2:157283/143‑1 (MQ=255)
                                                                 gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAACTCAGGCTTGCACATGC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCCAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTgc                                                           >  1:117858/1‑142 (MQ=255)
                                                                        aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAGAATCATGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc                                                    >  1:434451/1‑143 (MQ=255)
                                                                          tACCCGGCACACAAACGGCAATGAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                  <  2:795108/143‑1 (MQ=255)
                                                                                          aGCAATTAAAATGGTGAGTAGGATAAAAACCACGCTTGCCGACAC‑‑CTCAGCCGTGTTCCCCAACATCAACGCCAGATCAACGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                  <  1:387336/143‑1 (MQ=255)
                                                                                                            tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAGTCAGGGACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGACTGCTGAGCGATGGGCGTTCAAccatcca                >  1:45990/1‑143 (MQ=255)
                                                                                                                 tAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAACAAGGCCAATCAGGAACTCCTGACCGCCGTACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCAATGGGCGTTAAACCATCCatgatg           >  2:587129/1‑143 (MQ=255)
                                                                                                                     aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGAGCCTGCTTAGCGATGGGCGTAAAACCATCGATGATGGGtg       >  1:591778/1‑142 (MQ=255)
                                                                                                                     aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATAAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGAAAAACCATCCATGATGGGtt       >  1:810592/1‑143 (MQ=255)
                                                                                                                          ggCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc  <  1:840989/143‑1 (MQ=255)
                                                                                                                                         ||                                                                                                                                
TCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATC  >  NC_000913/2173224‑2173490

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: