Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 2,173,363 Δ2 bp intergenic (‑1/+1) gatC ← / ← gatC pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA2,173,3610C.97.4% 162.8 / ‑0.4 39pseudogene (2/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (21/16);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01
*NC_000913_3_pae_tpiA2,173,3620C.97.3% 154.4 / ‑0.0 39pseudogene (1/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base . (20/16);  minor base A (1/0);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAC  >  NC_000913_3_pae_tpiA/2173219‑2173492
                                                                                                                                              ||                                                                                                                                  
aGATTGCCACGATGCACGGCGACCGCCATCGCCACGAAGAACCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt                                                                                                                               <  2:592991/149‑6 (MQ=255)
        aCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                     <  1:5720/149‑1 (MQ=255)
        aCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTAATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                     >  1:573796/1‑149 (MQ=255)
              cACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCCgccgca                                                                                                               >  2:72244/1‑149 (MQ=255)
                         ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCAACGCCGGATCa                                                                                                    >  1:366618/1‑149 (MQ=255)
                          cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                   <  2:543101/149‑1 (MQ=255)
                           aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                  <  1:280655/149‑1 (MQ=255)
                                 aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTAAGTGGGATAAACATCCGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAACAGCAACGCCCGGTCAATGCCAAt                                                                                            >  2:97455/1‑149 (MQ=255)
                                       aagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                      <  1:417718/149‑1 (MQ=255)
                                       aagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGATTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                      <  1:183165/149‑1 (MQ=255)
                                       aacgctaTGGTGGCAAGATCGCCAAAGGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                      <  2:471053/143‑1 (MQ=255)
                                         gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCCGCAACGCCGGATCAAAGCCAACCAGGCACt                                                                                    >  2:116132/1‑149 (MQ=255)
                                         gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                    >  1:880074/1‑149 (MQ=255)
                                          ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                   <  1:655674/149‑1 (MQ=255)
                                           cGATGGTGGCAAGAACGCCAAACGACAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAAGCCAATAAGGAACTcc                                                                                  >  1:809829/1‑149 (MQ=255)
                                              tggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                               >  1:96510/1‑149 (MQ=255)
                                                  ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                           <  2:880074/149‑1 (MQ=255)
                                                   gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACCCAAACAAATAGCAAAATGGTGCGTGGTAAAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTACTGAccgcc                                                                          >  1:460024/1‑149 (MQ=255)
                                                   gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAACATCAGGCTTGCCGATACACCCGCCGTCTGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAGCTCCTGAccgcc                                                                          >  1:313745/1‑149 (MQ=255)
                                                          ctgtaaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                   <  2:140348/145‑1 (MQ=255)
                                                           gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGATGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCCCCGAACTTCg                                                                  >  1:736419/1‑149 (MQ=255)
                                                             cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATGAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                >  1:813170/1‑149 (MQ=255)
                                                              aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTAGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAATTCGCCt                                                               >  1:240138/1‑149 (MQ=255)
                                                                    cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGGAAAc                                                         >  2:274364/1‑149 (MQ=255)
                                                                         cactgttACACGGCACACAAACAGAAATTAGATGGGTGAGTTGTATAAAAAACAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCGGTAGACGTCTa                                                    <  2:398788/144‑1 (MQ=255)
                                                                                aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCACGAACTCCTGACCGCCGACATTCGCCTGTAAACGACTACGaacac                                             >  2:21268/1‑145 (MQ=255)
                                                                                 cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCtc                                            >  1:675648/1‑148 (MQ=255)
                                                                                               aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTGCGAGCCTGCTTAGCGATGGGCg                              >  1:391865/1‑149 (MQ=255)
                                                                                                       aaaTGGTGAGTGGGATAAAAATCAGGCTTGCAAGTA‑‑AGACCGCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTAACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTGAACGAGGGGCGTTAAAcac                      >  2:202728/1‑147 (MQ=255)
                                                                                                       aaaTGGTGAGTGGGATAAAAACCAGGCTTGGCGTTACCA‑‑CCGCCGTATGTCCCAGGAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcca                      <  1:61816/149‑1 (MQ=255)
                                                                                                            gTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                 >  1:344242/1‑149 (MQ=255)
                                                                                                              gAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAGCGACGAACTTCGCCTGTAACCGACTACGAACCTGCTTGGCGATGGGGGTTAAAccaccactga               >  2:842600/1‑142 (MQ=255)
                                                                                                                gTGGTATAACAACCCGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg             <  2:313745/149‑1 (MQ=255)
                                                                                                                gTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg             <  1:353145/149‑1 (MQ=255)
                                                                                                                gTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg             <  1:449362/149‑1 (MQ=255)
                                                                                                                gTGGGATAAAAACCAGGCTTGCCGATAC‑‑CGCCGGCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg             <  2:61811/149‑1 (MQ=255)
                                                                                                                 tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                            <  2:401712/101‑1 (MQ=255)
                                                                                                                 tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                            >  1:401712/1‑101 (MQ=255)
                                                                                                                       aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGTCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa      <  2:247095/149‑1 (MQ=255)
                                                                                                                           aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGAAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCAATGATGGGTTTAATCAc  >  2:778163/1‑149 (MQ=255)
                                                                                                                                              ||                                                                                                                                  
AGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAC  >  NC_000913_3_pae_tpiA/2173219‑2173492

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