Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_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_vch_tpiA2,173,3610C.95.8% 203.8 / 0.4 48pseudogene (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 (1/1);  new base . (24/22);  total (25/23)
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.56e-01
*NC_000913_3_vch_tpiA2,173,3620C.95.8% 203.1 / ‑0.3 48pseudogene (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 (0/0);  major base . (24/22);  minor base G (1/1);  total (25/23)
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.98e-01

GAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGC  >  NC_000913_3_vch_tpiA/2173215‑2173505
                                                                                                                                                   ||                                                                                                                                               
gAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGCCACACAAACGGCGTACAAAATGGTGAATGGGAGA‑AAAAGCAGGCTTACCGATaccaccg                                                                                                                                                >  1:607743/1‑148 (MQ=255)
     gATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATaccaccgccgta                                                                                                                                           <  1:727321/149‑7 (MQ=255)
            aCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                                  <  1:870432/149‑1 (MQ=255)
               aTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTCGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                               <  1:67631/149‑1 (MQ=255)
                     ggCGACCGCCATCGCCACTAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCAAACAAGCTGCAATCAGAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                                                                                                         >  2:899179/1‑149 (MQ=255)
                      gCGACCGGCATAGCCACGAAGAAGCCGATGGTGGCAAGAGCGACAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGCATA‑AAAATCAGGCTTGACGACACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                                        <  1:567201/149‑1 (MQ=255)
                      gCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                                        <  2:333743/149‑1 (MQ=255)
                          ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                                    <  2:310372/149‑1 (MQ=255)
                          ccGCCACAGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                                    <  1:985601/149‑1 (MQ=255)
                           cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                   >  2:69862/1‑149 (MQ=255)
                              cATCGCCACGAAGAAGCCGATGGTGGCAATATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTTGGATA‑AAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                                <  2:764100/149‑1 (MQ=255)
                                tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                              >  1:367599/1‑149 (MQ=255)
                                     aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCACCGCCGGATCACGGCCAAt                                                                                                         >  2:607636/1‑149 (MQ=255)
                                        aagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCCCATTTCCTAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                      <  1:1031912/149‑1 (MQ=255)
                                           aagCCGATGGTGGCAAGAACGCCAAACGGCAGCACCTGACTACCCGGTACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCTCCGTATGTCCCAGCAGCCACGCCGGATCAAGGCCAATCAGGaa                                                                                                   <  1:798584/149‑1 (MQ=255)
                                             gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTAATTACCCGGCACACAAACAGCAATTAAAACGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCCGCAACGCCGGATCACGGCCAGTCAGGAACt                                                                                                 >  1:296005/1‑149 (MQ=255)
                                               cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACGCAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc                                                                                               >  1:942078/1‑149 (MQ=255)
                                               cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc                                                                                               >  1:1011218/1‑149 (MQ=255)
                                                     tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                         <  1:69862/149‑1 (MQ=255)
                                                          aGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                                    <  2:1011218/149‑1 (MQ=255)
                                                                 cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAGATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                             >  1:758331/1‑149 (MQ=255)
                                                                 cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                             >  2:45549/1‑149 (MQ=255)
                                                                 cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                             >  1:452115/1‑149 (MQ=255)
                                                                 cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                             >  1:370560/1‑149 (MQ=255)
                                                                 cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCAcc                                                                             >  1:452116/1‑149 (MQ=255)
                                                                 cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCAAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                             >  1:192497/1‑149 (MQ=255)
                                                                          gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                                    <  2:370560/149‑1 (MQ=255)
                                                                              cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGCTA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCACGCCGCACCTCCCCTGTACACAACTAc                                                                >  2:401442/1‑149 (MQ=255)
                                                                              cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                >  2:531826/1‑149 (MQ=255)
                                                                                 aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                             >  2:394694/1‑149 (MQ=255)
                                                                                     cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATACAAGATGCTG‑ATGGCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTg                                                         <  1:607636/149‑1 (MQ=255)
                                                                                        ggCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACAACTACTAACCTGCtt                                                      >  1:62225/1‑149 (MQ=255)
                                                                                          cacacaAACAGCAATTAAAATGGTTAGTGGGATA‑AAAATCAGGCCTGCCGATACCA‑‑CCGCCGTATGTCCCGGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                                    <  2:551942/149‑1 (MQ=255)
                                                                                           acacaAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                                   <  2:367599/149‑1 (MQ=255)
                                                                                             acaAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                                 >  2:916726/1‑149 (MQ=255)
                                                                                               aaaCAGCAATTAAAATGGTGAGAGGGATA‑AAAATCAGGATTGCCGATACCA‑‑CCGCCGTATGTCAAAGAAGCAAAGCAGGATAAAGGCCAATAAGGAACTCCTGACGGCCGAACTTCGCCGGTAAACGACAAAGAGCATGATCAGCGATg                                               >  2:188965/1‑149 (MQ=255)
                                                                                                        ttAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                      >  2:331415/1‑149 (MQ=255)
                                                                                                         tAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                                     <  1:1647/149‑1 (MQ=255)
                                                                                                         tAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                                     <  2:520843/149‑1 (MQ=255)
                                                                                                            aaTGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                  <  1:331415/149‑1 (MQ=255)
                                                                                                                  gAGTGGGATA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                            <  1:45549/149‑1 (MQ=255)
                                                                                                                         aTA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                     >  1:495587/1‑149 (MQ=255)
                                                                                                                          tA‑AAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCGGCAGCAACGCCGGATCAGGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGCGCCTGCTTAGCGATGCACGTTACACACTCCATGATggttgta                    >  1:270067/1‑144 (MQ=255)
                                                                                                                           a‑aaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                   <  2:648271/149‑1 (MQ=255)
                                                                                                                           a‑aaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                   >  2:862092/1‑149 (MQ=255)
                                                                                                                           a‑aaaaTCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATAAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                   <  2:212215/149‑1 (MQ=255)
                                                                                                                                aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAAGCGCCCGCTTAGCGATGGGCGTTAAACCCTCCACGATGGGTTTAATCAc               >  2:375643/1‑149 (MQ=255)
                                                                                                                                     gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                        <  1:1014961/71‑1 (MQ=255)
                                                                                                                                     gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                        >  2:1014961/1‑71 (MQ=255)
                                                                                                                                     gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGTCTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg          <  1:675044/149‑1 (MQ=255)
                                                                                                                                               atacccccGCCGTATTTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGACCCTGCTTCGCGATGGGCGTTAAACCACCCATGATGGGTTTAATCACCCGTGGCATTagc  <  2:927525/145‑1 (MQ=255)
                                                                                                                                                   ||                                                                                                                                               
GAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATA‑AAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGC  >  NC_000913_3_vch_tpiA/2173215‑2173505

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