Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 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_0009132,173,3610C.96.4% 107.2 / ‑0.3 29pseudogene (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);  major base . (14/13);  minor base T (1/0);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01
*NC_0009132,173,3620C.96.3% 100.8 / 0.8 29pseudogene (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/1);  major base . (14/12);  minor base A (1/0);  total (15/14)
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.89e-01

TTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173223‑2173496
                                                                                                                                          ||                                                                                                                                      
aaccACGAGGGACGGGGACCGCCAGCGCCACGAAGAAGCGGATGGGGGCAAGAGCGCCAAGCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgcc                                                                                                                                     <  1:42893/141‑4 (MQ=255)
  ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt                                                                                                                                   <  1:527301/143‑6 (MQ=255)
          cccGGCGACCGCCATCGCCACGAAGACGCCGTCGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                         <  2:94055/141‑1 (MQ=255)
             ggagACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                      <  1:212150/140‑1 (MQ=255)
                        tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACCAACAGCAATTAAAACGGTGAGTGGGCTAAAAATCAGGCTTGCCGGTACCA‑‑CCGCCGTACGCCCCAGCAGCAACGcccaa                                                                                                           >  2:285156/1‑140 (MQ=255)
                        tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                           >  2:618072/1‑143 (MQ=255)
                             aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGACGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                      >  1:384462/1‑143 (MQ=255)
                             aCGAAGAACCCGCTGGTGCCAAGATCGCCAAACGGCAGCACCTCATTACCCGGCACACAAACAGCAATTCAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                      <  1:134932/143‑1 (MQ=255)
                                     gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                              <  1:80476/143‑1 (MQ=255)
                                         agggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCACTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGACCAAGGCCAATCAgcac                                                                                          >  2:146224/3‑140 (MQ=255)
                                              ggCAAGATCGCCAAACGGCAGCACCTGACTACCCCGCAAAAAGACAGCCATTTAAGTGGCGTGTGGGATCAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTACCTGCAGTAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                     <  2:171192/143‑1 (MQ=255)
                                                       gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAGAAATCAGGCTTCGCGATACCA‑‑CCGCCGTAGGTCCCAGCAGCAACGCCGGAGCAAGGCCAATCAGGGACGCGTGAccgccgc                                                                            >  2:468979/1‑142 (MQ=255)
                                                        ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAGCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAGGGACTCCTGAcagcacga                                                                           >  1:269491/1‑139 (MQ=255)
                                                         cAAACGGCAGCACCTGATTACCCGGCGCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                          <  1:470213/143‑1 (MQ=255)
                                                              ggCAGCACCTGATTACCCGGCACACAAACATAAATGAATATGGTGAGCGGGAGCTACAACAGGAGGGCACGTCCATACTGGTTATGTGCCAGCAGGAGCGCCGGATCAAGGCCCATCGGGAACTCCTGACCGCCGAACGTCGc                                                                     >  2:657387/1‑143 (MQ=255)
                                                                cagcaCCTGATTGCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                   <  1:511765/143‑1 (MQ=255)
                                                                             cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGCCTaag                                                      >  1:485234/1‑141 (MQ=255)
                                                                               cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                    >  1:355790/1‑143 (MQ=255)
                                                                                ggCACACAAACAGCAATTAAAATGGTGAGTGAGATAAAAATCAGGCTTGCCGATACCA‑‑CAGACGTATGTCACAGCAGCAACGCCGGATCAAGGACCATCCGGAACTCCTGACCGCCGAACTTCGCCCGTAAACGACtctaagc                                                   >  2:468835/1‑137 (MQ=255)
                                                                                 gCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCATGCTTGCCGATACCA‑‑CCGCCGTATGTCCCCCCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGcc                                                  <  1:164846/143‑1 (MQ=255)
                                                                                                ctacagtGGGGAGTGGGATAGAAATCGGGCTGGCCGGTACCA‑‑CCGTCGTATGTCGCAGCAGCAGCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                   <  2:533516/137‑1 (MQ=255)
                                                                                                   aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGACTGTAAACGACTACGAGCCCGCTTAGCGATGGGCgct                                >  2:421324/1‑141 (MQ=255)
                                                                                                           aGTGGGATAAAAATCAGGCTTGCCGATACAGGCGCGGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                        <  2:295452/143‑1 (MQ=255)
                                                                                                                ggtAAACATCAGGCTGCCCGATCCCA‑‑CCGCCGTGTGTCCCAGCAGCAGCGCCGGATCAAGGCCAATCACCAGCGGCTGACCGCCGAAGTTGGCCTGTAAACGACGACGAGCCTGCGTAGCGATGGGCGTTAAACCATCCatga                   <  2:463662/141‑1 (MQ=255)
                                                                                                                gacaagaaTCATTCTTCCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCCACGGCGGATCAAGGCCAATCAGGCACTCCTGTCCGCCGAACTTCGCCTGTAAACGACTACTAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                   <  2:205090/137‑1 (MQ=255)
                                                                                                                 aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCCCCGACCTTCGCCTGTCAACCACTACGAGCCTGCTTACCGATGCGCCTCAAACCCTCcatcat                  >  2:87542/1‑143 (MQ=255)
                                                                                                                   aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAGCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCCTGCTTAGCGATGGGCGTTACAGCATCCACGATgg                >  2:204349/1‑143 (MQ=255)
                                                                                                                    agaagcAGGCTCGCTGATACCA‑‑CCGCCGGATGTCCCAGCACTAACGCCGGATCAAGGCGAATCAGGAACTCCTGACCGCCGAATTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg               <  1:239267/138‑1 (MQ=255)
                                                                                                                     aaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt              >  2:447214/1‑143 (MQ=255)
                                                                                                                        tCAGGCTTGCCGGTACCA‑‑CCGCCGGATGTCCCAGCAGCGACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTa           <  2:313035/143‑1 (MQ=255)
                                                                                                                                 ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGGCTGCTTAGCGATGGGCGTTCAACCATCCATGATGGGTTTCATCACCCGt  >  1:658640/1‑143 (MQ=255)
                                                                                                                                          ||                                                                                                                                      
TTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173223‑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: