Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,862,676 G→C 32.7% pseudogene (927/1617 nt) glvC ← pseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6760GC32.7% 15.8 / 4.9 32pseudogene (927/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base G (6/14);  new base C (3/6);  total (12/20)
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.72e-01

CGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAAC  >  NC_000913/3862456‑3862914
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCTTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGCTTAGGTTTCGGGGATACGCCAAACCCCCACtttttaccca                                                                                                                                                                                                                                >  1:51704/1‑233 (MQ=255)
   ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTATGTTGCGGGGAACACAAAGACCCCCCACTTTTTcccactt                                                                                                                                                                                                                              >  1:471571/1‑235 (MQ=255)
    gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTTTGGGGAACAGACAACCCCGCACCTTTTAccctgc                                                                                                                                                                                                                               >  1:875311/1‑234 (MQ=255)
                     tCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGTTTCGGGAAACCGCAAGCCCCGCCCCTTTTACCCCGTTTTCCgggcgcgcggggaa                                                                                                                                                                                                            >  1:7231/1‑227 (MQ=255)
                             tGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGGATCAGCCAAACCGCCCCCTTTACCCCAATTTTCGTTGCGGGCGGGGAGTACAACg                                                                                                                                                                                                    >  1:64891/1‑238 (MQ=255)
                                                tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGATTAGGGGAACAGCAAAGCCACCCCATTTTACCCGGTTTTCAGGCCGGCAGGGAAGAACATACCCAAAAAAAAGCCCACGGc                                                                                                                                                                                >  1:592437/1‑238 (MQ=255)
                                                             tgtGTTCCCCCCCACAACAAAGAAATACTTTACGGGAAACAGTGGGGCCCCAAGCTCGGTGTGTCCCACAAACAGAAGCGGTAAAAGAAACGGGAGTGTAAATCCAAGCGGTTCGGTAATCCCCACCAGCTGGGGGGTTAGGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTg                                                                                                                                                                  <  2:556607/237‑1 (MQ=255)
                                                                                     aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTCGGGGATCAGCAAACCCGCCCCCTTTTCCCCGTTTTTCGGGGCGGCAGTGGAGGACATAGCTAAAAAAATAGCCACGGCCCCACAGAATTTTTAGGTTTCGTGGAGGGGAAAAACg                                                                                                                                          >  1:877167/1‑238 (MQ=255)
                                                                                      tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGGAATTCCCACCAGCATGGCGGGTTAAGGTTCGGGGGACGAGCAACCCCCCCCCTTTTCCCCGATTTTCCGTGGGGGCGAGAAGGAAAACGCCTAAAAAAAGCCCCCGAGCCCCAAAAGTTTTTAGTGAACCTGGAGGGCaaaaacac                                                                                                                                         >  2:158676/1‑234 (MQ=255)
                                                                                                                       gggtgtCCCGAAAACAGAAACTGGAAAATGAACAGGAAGGTAAATCCCAGCGGTTCGGAAATGCCCCCCAGCAGGGGGGTTAAGGTTGGGGGAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGgcgc                                                                                                        <  1:436192/236‑1 (MQ=255)
                                                                                                                                  aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTCAAGCTGCGGGGAGCCGGCAGCCCCCCCCCTTTTAACCAGTTTTCGGGGCGGCAGTGAAGTTCACTCCTAAAAAAAAGCCCCCCGGCCCCAAGAACTTTTAGTTGACCGTCCCGGCACAAACGCCCTCCCGGGACAAACAATTTCacggacccgcgcacaaaacaaac                                                                                              >  2:358178/1‑216 (MQ=255)
                                                                                                                                             aggaaagaaacgagAAGGTAAATCCCAGGGTTTCTATAATCCCCACCAGCATGGCGGTTAAGGTTGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                                  <  1:158676/225‑1 (MQ=255)
                                                                                                                                               gaaagaaagaGGAAGGTAATTTCCAGGGGTCCGGAAACCCCCCCCAGCAGGGCGGTTAAGGTTGCGGAAACACGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                                <  1:491366/229‑1 (MQ=255)
                                                                                                                                                 aaggaaagggAAGTAAAATCCCATGGGTTCGGTAATCCCCCCCAGCAGGGCGTTTAAGGTGGGGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                             <  1:668578/231‑1 (MQ=255)
                                                                                                                                                    gaaaggaaaGATAAATCCCGGGGGTCGGGTAATCCCCCCCAGGTGGGGGTTTAGGGTGGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                           <  1:189262/231‑1 (MQ=255)
                                                                                                                                                      aagggaatttAAATTCCGGGGTTTGGTAAATCCCCCAACGCAGGGGTGTTAGGGTTGGGGAAATAGGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                                        <  1:821282/230‑1 (MQ=255)
                                                                                                                                                          gAAGGAAAATTCCGGGGTTCCGGTAATCCCCACAAGCTGGGCGTTTAAGGTGGCGGAAATCAGCAGCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTa                                                                      <  1:84242/237‑1 (MQ=255)
                                                                                                                                                                   ctcacagggtttgtgaaTACCCACCCAGCTGGGCGTTTAAGTTTGCGGGAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTCCCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATg                                                            <  1:579891/224‑1 (MQ=255)
                                                                                                                                                                               cgTAATCCCCCACAGGATGGGCGGTTAGGTTTGGGGAAATAAGCAGGCCCGCCCCTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAAc                                                 <  1:498370/237‑1 (MQ=255)
                                                                                                                                                                                        cacgaGCATGGCGTTGAGGTGTGCGGAAATCAGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                                      <  1:358178/239‑1 (MQ=255)
                                                                                                                                                                                         ccccaGCTTGGCGTTTAGGGTTGGGGGAATCGGAATGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTTCATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                                      <  1:672513/237‑1 (MQ=255)
                                                                                                                                                                                          cccccgcgtGGCGGTTAGGGTGGGGAAAAGCAGCACGCCCGCCACGTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAGGCCCACGGCCCCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                                      <  2:247855/231‑1 (MQ=255)
                                                                                                                                                                                              ccaacgtggggattAAGGTGGCGGGAACGAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAATTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                                    <  1:526147/227‑1 (MQ=255)
                                                                                                                                                                                                  gcggggggTTAAAGTTGGGGGAAATCCGCAACCCCGCCACTTTGACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAg                                 <  1:430526/233‑1 (MQ=255)
                                                                                                                                                                                                                ttGCGGGAAACAGCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt                    <  1:471910/237‑1 (MQ=255)
                                                                                                                                                                                                                    cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGAta               >  1:624129/1‑238 (MQ=255)
                                                                                                                                                                                                                    cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCCCGGCGCCCAAAACGAACTCTTCCAGCACCGGCCCCCCGATCATCTGAATTCCCCCTTCAACCACTGCCCGGCCAAAAATaaaactgccgtacatg               >  2:861609/1‑226 (MQ=255)
                                                                                                                                                                                                                         gaaataGCATGCCCGCCAGTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACAGCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa             <  1:861609/233‑1 (MQ=255)
                                                                                                                                                                                                                         aaaacAGGAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa             <  1:182513/234‑1 (MQ=255)
                                                                                                                                                                                                                           gaaccaagaaGCCCGACACTTTTGCCCGATTTTCCGGTCGGGCAGTGAAGTACATCGCTAAAGAAAGGCCCTCGGCACAAAAGATTCTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa             <  1:835374/230‑1 (MQ=255)
                                                                                                                                                                                                                           aaacgcagacgCCCGCCACGTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGaa             <  1:134641/228‑1 (MQ=255)
                                                                                                                                                                                                                           aTCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCCCAAACTGAACTCTTGCCGATGCTGGGCCCAGTACACCTGACTGCCGCCTTCAAACGCTGCCGGACCAAAGAGAAACTTGCCCTAGATGAAggg          >  1:286019/1‑237 (MQ=255)
                                                                                                                                                                                                                            attaacctaGCCCCCCCCTTTTACCCGATTTTCCGGTGCGGGAGTGAAGTACAGCGCTAAAGAAAGGCCCCCGGCACCAAAGATTTTTGAGTTACCGGGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAg            <  2:624129/230‑1 (MQ=255)
                                                                                                                                                                                                                                     gCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCACCAGCTGCCGGACCAAAGATAAAATGTCCGCAGATGAAGTGGTGTAAAc  >  1:451896/1‑239 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
CGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAAC  >  NC_000913/3862456‑3862914

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