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.96.8% 128.0 / ‑1.9 31pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (14/16);  total (15/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 = 5.77e-01
*NC_0009132,173,3620C.96.8% 127.5 / ‑2.3 31pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (14/16);  total (15/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 = 8.71e-01

AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173219‑2173496
                                                                                                                                              ||                                                                                                                                      
aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGCTGGTGGCCCGATCGCCAAACAGCAACAACAAATTAACAAGAAAAAAAACAGCAAAAAAAATGGTGAGAGGGATAAAAATCAGGCTTGCCGATaccac                                                                                                                                         >  1:609418/1‑143 (MQ=255)
        aCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                                 >  1:481129/1‑136 (MQ=255)
              cACGGCGACCGCCATCGCCACTAAGAAGCCGATGATGGCAATATCGCCAAACGGCAGCATCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                         <  1:25722/143‑1 (MQ=255)
                 tgCGACCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                      <  1:482753/142‑1 (MQ=255)
                 ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATAGAAAGACGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                      >  1:337309/1‑143 (MQ=255)
                     aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGACTTACCAATACCA‑‑CCACCGTATGTCCCAGCAGCaa                                                                                                                  >  1:303520/1‑143 (MQ=255)
                       cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTTCCGATACCACCC‑‑CCGTATGTCCCAGCAGCAACg                                                                                                                >  1:562065/1‑143 (MQ=255)
                                 aCGAAGAAGCCGCTGGTGGCAGGCTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                      <  2:749920/143‑1 (MQ=255)
                                      gaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTAGGTCCCAGCAGCCCCGCCGGATCAAGGCCaa                                                                                                 >  1:487521/1‑143 (MQ=255)
                                      gaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGACTTGCCGATACCA‑‑CCACCGTATGTCCCAGCAGCAACGCCAGATCAAGGCCaa                                                                                                 >  1:260206/1‑143 (MQ=255)
                                          ccGATGGAGGTAAGATAGCAAAACGGCAGCACCAGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                             <  1:444146/143‑1 (MQ=255)
                                                gtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                       <  2:481129/143‑1 (MQ=255)
                                                  ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTACGTACCAGCAGCAACGCCGGATCAAGGCCAATCAGAAACTccc                                                                                     >  2:574066/1‑142 (MQ=255)
                                                   gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGAGAGAGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAACAAGGCCAAGAAGGacaccctg                                                                                    >  2:628526/1‑136 (MQ=255)
                                                     aaGCTCGCCACACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                  <  2:729175/143‑1 (MQ=255)
                                                         tCGCCAAACGCCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAACAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                              <  1:422707/143‑1 (MQ=255)
                                                                  ggCAGCACCTGATTACCCGGCGCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                     <  2:761053/143‑1 (MQ=255)
                                                                   gcagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATCGTAAGGTCGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCACAGCAGCAACGGCGGATCAAGGCCAATAAGGAACTCCTGAGCACCGAACTTCGcc                                                                    >  2:432723/1‑143 (MQ=255)
                                                                        aCCTGATTACCCGGCACACAAACAGCAATTAACATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCCGATCAAAGCCAATCAGGAACTCATGACCGCAGAACTTCGCCTGTac                                                               >  1:185301/1‑142 (MQ=255)
                                                                           tGATTACCTGGCACACTAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACg                                                            <  1:745613/143‑1 (MQ=255)
                                                                            gATTACCCGGCATACAATCAGCAATTACAATGGTGAGTGTGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                           <  1:66719/143‑1 (MQ=255)
                                                                            gATTACCCCGCCGACAGACAGCAATTATTATGTTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                           <  2:399705/143‑1 (MQ=255)
                                                                                    ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCAACGCCGGATCAAGTCCAATCAGGAAATCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                   >  1:589279/1‑143 (MQ=255)
                                                                                                  aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAA‑‑CCGCCGTATGTCCCAGCAGCAACGCAGGATCAAGGCCAATCAGGAACTCCGGACCGCCGGACTTCGCCTGTACACGACAACGAACCTGCTGCGCGATga                                     >  2:298653/1‑142 (MQ=255)
                                                                                                   aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                         <  1:557439/74‑1 (MQ=255)
                                                                                                   aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                         >  2:557439/1‑74 (MQ=255)
                                                                                                    ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTCATCGACTACGAGCCTGCATAGCAATGAGc                                   >  1:27772/1‑143 (MQ=255)
                                                                                                    ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTAACCGCAGAAATTCACCTGTAAACTACTACAAGCCTGCTTAGCGATAGGc                                   >  1:536887/1‑143 (MQ=255)
                                                                                                                   ggATAACAATCAGGCTTGCCGACACCA‑‑CCGCCTTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                    <  1:264355/143‑1 (MQ=255)
                                                                                                                       aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                <  1:742900/143‑1 (MQ=255)
                                                                                                                          cgtCAGGCCTGCCGATATCA‑‑CCGCCGTACGGCCCAGCAACAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGGAAACGCCTACGAGCCGGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt             <  1:544846/141‑1 (MQ=255)
                                                                                                                                gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa       <  1:583194/143‑1 (MQ=255)
                                                                                                                                     ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCGGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACGAGCCTGCTTCGCGCTGGGCGTTAAACCACTCATCATGTGTTATATCACCAGt  >  1:380854/1‑143 (MQ=255)
                                                                                                                                              ||                                                                                                                                      
AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173219‑2173496

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: