Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,862,518 Δ1 bp pseudogene (1086/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,5160C.100.0% 186.2 / NA 46pseudogene (1088/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 C (0/0);  new base . (23/23);  total (23/23)

GAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATG  >  NC_000913/3862379‑3862650
                                                                                                                                         |                                                                                                                                      
gAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                       >  2:236008/1‑139 (MQ=255)
   tCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                                                                    <  1:224714/139‑4 (MQ=255)
        gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                                                                              <  1:830226‑M1/139‑4 (MQ=255)
        gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                                                                              <  1:236008‑M1/139‑4 (MQ=255)
         tgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCAGGATCCAGTTTTGCGGTAAAAATTAGTCAATCAGAACTCCGCCATGTTcacc                                                                                                                             >  2:1287771‑M1/1‑135 (MQ=255)
             gAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacac                                                                                                                         >  2:837566‑M1/1‑131 (MQ=255)
             gAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacac                                                                                                                         <  2:1392806‑M1/139‑9 (MQ=255)
                  ccACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaag                                                                                                                    <  1:964768‑M1/139‑14 (MQ=255)
                        aGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCTGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                                                                    <  2:1242169/118‑4 (MQ=255)
                        aGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCGGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                                                                       >  1:829349‑M1/1‑120 (MQ=255)
                        aGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCGGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                                                                       <  2:829349‑M1/130‑11 (MQ=255)
                        aGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                                                                    >  1:1242169/1‑115 (MQ=255)
                         gTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagataca                                                                                                             <  2:393357‑M1/139‑21 (MQ=255)
                                  ggTAAAGCACAACCCAATGGTGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcg                                                                                                    >  2:334445‑M1/1‑110 (MQ=255)
                                    tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCTATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                       <  1:1016985/102‑1 (MQ=255)
                                    tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCTATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                       >  2:1016985/1‑102 (MQ=255)
                                            aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggc                                                                                          >  1:1536828‑M1/1‑100 (MQ=255)
                                              cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccg                                                                                        <  2:521457‑M1/139‑42 (MQ=255)
                                                      cGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                                                <  1:837566‑M1/139‑50 (MQ=255)
                                                        aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                                              <  1:334445‑M1/139‑52 (MQ=255)
                                                        aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                                              >  2:125820‑M1/1‑88 (MQ=255)
                                                        aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacac                                                                                                                         <  2:179684‑M1/96‑9 (MQ=255)
                                                        aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacac                                                                                                                         >  1:179684‑M1/1‑88 (MQ=255)
                                                          ttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgt                                                                            <  1:343233‑M1/139‑54 (MQ=255)
                                                               tCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtacc                                                                       >  1:1069979‑M1/1‑81 (MQ=255)
                                                               tCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcacca                                                                                                                            >  2:840553‑M1/1‑81 (MQ=255)
                                                               tCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcacca                                                                                                                            <  1:840553‑M1/86‑6 (MQ=255)
                                                                             gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                                                <  2:407660‑M1/116‑50 (MQ=255)
                                                                             gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                                                >  1:407660‑M1/1‑67 (MQ=255)
                                                                             gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaag                                                                                                                    <  1:1207570‑M1/80‑14 (MQ=255)
                                                                             gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaag                                                                                                                    >  2:1207570‑M1/1‑67 (MQ=255)
                                                                             gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                                                                       <  1:229077‑M1/77‑11 (MQ=255)
                                                                             gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                                                                       >  2:229077‑M1/1‑67 (MQ=255)
                                                                              cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggt                                                        >  2:92927‑M1/1‑66 (MQ=255)
                                                                              cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                                                   <  1:608713‑M1/112‑47 (MQ=255)
                                                                              cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                                                   >  2:608713‑M1/1‑66 (MQ=255)
                                                                                   ttGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaat                                                   >  2:1289635‑M1/1‑61 (MQ=255)
                                                                                             tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaag                                         >  1:803794‑M1/1‑51 (MQ=255)
                                                                                              cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                                     <  2:1468671‑M1/126‑77 (MQ=255)
                                                                                              cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                                     >  1:1468671‑M1/1‑50 (MQ=255)
                                                                                              cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCAAAGAGATACAtt                                                                                                        >  2:759005/1‑72 (MQ=255)
                                                                                              cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCAAAGAGATACAtt                                                                                                        <  1:759005/72‑1 (MQ=255)
                                                                                                  aTCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCTGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaa                                    <  1:428819‑M1/139‑94 (MQ=255)
                                                                                                   tCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaat                                   >  1:1360723‑M1/1‑45 (MQ=255)
                                                                                                   tCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                                              >  2:501591‑M1/1‑45 (MQ=255)
                                                                                                   tCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                                              <  1:501591‑M1/96‑52 (MQ=255)
                                                                                                      aGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattcc                                <  1:1289635‑M1/139‑98 (MQ=255)
                                                                                                                     cttagTCAATCAGACCTCCGCCATATTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaatt                 <  1:125820‑M1/137‑113 (MQ=255)
                                                                                                                        aGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaattccc              >  1:1405855‑M1/1‑24 (MQ=255)
                                                                                                                               cAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaattcccaccagca       >  2:243429‑M1/1‑17 (MQ=255)
                                                                                                                                    cTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaattcccaccagcatggcg  >  1:902205‑M1/1‑12 (MQ=255)
                                                                                                                                    cTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcgg                      <  1:465836‑M1/119‑108 (MQ=255)
                                                                                                                                    cTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcgg                      >  2:465836‑M1/1‑12 (MQ=255)
                                                                                                                                         |                                                                                                                                      
GAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATG  >  NC_000913/3862379‑3862650

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