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% 114.4 / ‑2.7 28pseudogene (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 . (15/12);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
*NC_0009132,173,3620C.96.4% 113.7 / ‑3.1 28pseudogene (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 . (15/12);  minor base G (0/1);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173221‑2173496
                                                                                                                                             ||                                                                                                                                      
tctgccACGATGCACCGCGGCCGCCATCGCTACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                       <  2:62776/139‑2 (MQ=255)
     cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTGGCCGATaccaccgccgta                                                                                                                                  >  1:256954/1‑137 (MQ=255)
            cACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                         >  1:123474/1‑143 (MQ=255)
                               tcggagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATTGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                      <  2:271509/139‑1 (MQ=255)
                                       gCCGATGGGGGCAAGATCGCCAAACGGCAGCAACTGATTACCCGGCACACAAACAGCAATAAAAAAGAT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCCGCAGCAACGGCGGATCACGGCCagaca                                                                                              >  2:180044/1‑139 (MQ=255)
                                         cGATGGTGGCAAGATCGCCAAACGGCAGCAACTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATAAgg                                                                                            >  2:264180/1‑143 (MQ=255)
                                          gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGG‑GAGTGGGATAAAAATCAGGCTTGCCGAGACCA‑‑CCGCCGTATGTCCCAGCAGCAAAGCCGGATCAAGGCCAAACAGGa                                                                                           >  1:359220/1‑143 (MQ=255)
                                                      aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATCAGGAACTCCTGAccgc                                                                               >  1:623464/1‑143 (MQ=255)
                                                       tCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                              <  1:351839/143‑1 (MQ=255)
                                                           cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                          <  2:28989/143‑1 (MQ=255)
                                                            acaCGGCAGCCCCTGATTACCCGGCACACAAACAGCAATGAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                         <  2:5691/141‑1 (MQ=255)
                                                            aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCGCGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCCCGCCGGACt                                                                         >  1:504800/1‑143 (MQ=255)
                                                            aaaCGGCAGCACCTGACTACCCGGCACACAAACAGCAATTAACATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                         <  1:583759/143‑1 (MQ=255)
                                                                ggCAGCACCTGATTCCCCGGCACACAAACAGCAAAT‑AGATGGTGGAGGTGGATAAAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                     <  2:165805/143‑1 (MQ=255)
                                                                ggCAGCACCTGATTACCCGGCACACAAACAGCAATCGTAATTGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                     <  2:165839/143‑1 (MQ=255)
                                                                    gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                 <  2:184022/143‑1 (MQ=255)
                                                                        cTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGCGCCAATGCCAATCAGGAACTCCTGACCGCCGGACTTCACCTGTaaaa                                                             >  1:335510/1‑142 (MQ=255)
                                                                            gtACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                         <  1:777599/142‑1 (MQ=255)
                                                                                ccGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCTGACCCCCGAACTTCGCATGTAAACCCCTACGa                                                     >  2:503641/1‑143 (MQ=255)
                                                                                ccGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTACCAGCAGCAACGCCGGATCAAGGCCTATCAGGACCTCCTGACCGCCGAACTTGCCCGGTAAACGACAACGa                                                     >  2:584761/1‑143 (MQ=255)
                                                                                ccGGCACACAAACAGCAATTAAAAAGGA‑GAGTGGGATAAAAATCAGGCTTGCCGATACAA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                     >  1:473586/1‑143 (MQ=255)
                                                                                        caAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCACAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTGCACCCGTAAACGACTACCAGcccggtg                                             >  2:591953/1‑138 (MQ=255)
                                                                                         aaaCAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                            >  1:135228/1‑143 (MQ=255)
                                                                                          aaCAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                           <  2:256954/143‑1 (MQ=255)
                                                                                                     aaaTGGT‑GAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGCCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                <  2:621927/143‑1 (MQ=255)
                                                                                                              attGGGAGAAAGTTCAGTCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCACAACCGCCGGAGCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                        <  2:92395/141‑1 (MQ=255)
                                                                                                                 gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTGCGCCTGTAAACGACTACGAGCCTGCTTCGCGATGGGCGTTAAAccatccat                     >  2:255164/1‑143 (MQ=255)
                                                                                                                          aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGGTAAACCCCCCATGATGGGttt            >  2:143615/1‑143 (MQ=255)
                                                                                                                             aGGCGTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCAGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACGTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt         <  2:69123/143‑1 (MQ=255)
                                                                                                                                    ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCCTAGCGATGGGCGTTAAACCATCCATGATAGGCTTAATCACCCGt  >  2:500718/1‑143 (MQ=255)
                                                                                                                                             ||                                                                                                                                      
ATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGT‑GAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173221‑2173496

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