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.95.0% 70.4 / 0.8 21pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (7/12);  minor base A (1/0);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01
*NC_0009132,173,3620C.95.0% 70.3 / 1.0 21pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (7/12);  minor base A (1/0);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

GGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATG  >  NC_000913/2173236‑2173478
                                                                                                                             ||                                                                                                                    
ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGTAATTAAAATGTTGATTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                                                                                    <  1:604859/129‑3 (MQ=255)
       gCCATCGCAACGGAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATGAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                           <  2:290521/129‑1 (MQ=255)
         cATCGCCACGAAGGAGCAGATGGTTGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcc                                                                                                         <  2:429209/129‑1 (MQ=255)
                aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAATATGGTGTGTGGGATAAAAATCAGGCTTGCAGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                  <  1:662127/129‑1 (MQ=255)
                   aagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCAAc                                                                                               >  2:418190/1‑129 (MQ=255)
                      tagCCCATTGTAGCACGATCGCCAAACGTCACCATCTGCTTGCCCGGCACACAACCAGCAATTAAAAGGGTGAGGGAGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                            <  2:8631/128‑1 (MQ=255)
                       agCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAATATTTTTTTTTTTCTATATCA‑‑CCGACGTATGTCCCAGCAGCAACGCCg                                                                                           <  1:154990/129‑1 (MQ=255)
                        gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTTGGATAAATTTTATTTTTGCTGTTAC‑‑TTTCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                          <  1:130951/129‑1 (MQ=255)
                            aTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTCTGTCCCAGCAGCAACGCCGGACCa                                                                                      >  2:524577/1‑129 (MQ=255)
                             tggtggCAAGATCGCCAATCTGCAGCCCCTGATTAACCGAAAAAAAAACAACAATAAAAAAAGAAAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CAGCCATATGTCCCAGCAGCAACGCCGGATaaa                                                                                     >  1:201940/1‑129 (MQ=255)
                               gtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                   >  2:483550/1‑129 (MQ=255)
                                               aCGGCCGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGTCGATATCA‑‑CTGCCGTATGTCCCAGTAGTATCGCGGGATCAAGGCGGATCAGGAACTCCt                                                                   <  1:465669/129‑1 (MQ=255)
                                                  gcagcaCCTGATTACCGGGCACACAAACAGCAATTAAAATGGCGAGTGGAAAAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                <  2:252538/129‑1 (MQ=255)
                                                         cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑TTGCCGTATGTTTCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                         <  1:368920/129‑1 (MQ=255)
                                                                                 aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGC‑‑AAAAAAAACAAACATATAACAAAGCAGCAACGCCGAATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCACCTGTAAAAGACTACGAGc                                 >  1:483619/1‑129 (MQ=255)
                                                                                      aaaTGGTGAGTGGGATAATAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                            <  2:293799/129‑1 (MQ=255)
                                                                                       agtGGTGAGTGTGAGAACGATCAGGGTCGGCGGTACCA‑‑CCGGCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGATCTCCTGACCGCCGAACTGCGCCTGTAAACGACTCCGAGCCTGCtt                           <  2:146878/127‑1 (MQ=255)
                                                                                            tgatttggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAATTCCTGATCAATTCCAATCATGATCTCCTGACCTCCGAACTTCGCCAGTAAACGACTACGAGCCTGCTTAGCGa                      <  1:626468/123‑1 (MQ=255)
                                                                                             gAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCACCGTATGTCCAAGCAACAACGCCAGAACAAGGCCAATCAGGAACTGGTGACCGCCAAACATCCCCACTAGACGACTACGAGCCTGCTTAGCGAt                     >  1:122869/1‑129 (MQ=255)
                                                                                               gTGGGATAAAAATCATGCTTGCCGATACCA‑‑CCACCGTAAGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                   >  1:491626/1‑129 (MQ=255)
                                                                                                  ggATAAAAATCAGGCTTGCCGCTACCA‑‑CCGCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACCTACCCTGTAAACGACTACGAGCCTGCTGAGCGAtacgcg                >  2:628511/1‑124 (MQ=255)
                                                                                                                cTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCTGCCACGCCAGATCAATGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAATCGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg  <  1:418190/129‑1 (MQ=255)
                                                                                                                             ||                                                                                                                    
GGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATG  >  NC_000913/2173236‑2173478

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: