Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,862,497 C→T pseudogene (1107/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,4970CT100.0% 194.8 / NA 61pseudogene (1107/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 T (30/31);  total (30/31)

GCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTC  >  NC_000913/3862362‑3862622
                                                                                                                                       |                                                                                                                             
gCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                            <  1:501250/139‑1 (MQ=255)
      gTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATc                                                                                                                      >  1:16079/1‑139 (MQ=255)
         gAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                            <  1:1305467/130‑1 (MQ=255)
         gAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                            >  2:1305467/1‑130 (MQ=255)
         gAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGa                                                                                                                   >  1:451233/1‑139 (MQ=255)
            cTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCt                                                                                                                <  2:451233/139‑1 (MQ=255)
              gCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAAt                                                                                                                       <  1:1336349/130‑1 (MQ=255)
              gCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAAt                                                                                                                       >  2:1336349/1‑130 (MQ=255)
               cAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCa                                                                                                                     <  1:594605/131‑1 (MQ=255)
               cAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCa                                                                                                                     >  2:594605/1‑131 (MQ=255)
                aGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGc                                                                                                            >  2:655809/1‑139 (MQ=255)
                 gAATCAGTGTACGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                           <  2:628768/139‑1 (MQ=255)
                  aaTCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                            >  1:1372767/1‑121 (MQ=255)
                  aaTCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                            <  2:1372767/121‑1 (MQ=255)
                    tCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                                             >  2:1188889/1‑134 (MQ=255)
                    tCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                                             <  1:1188889/134‑1 (MQ=255)
                       gtgtGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                           <  2:1067440/133‑1 (MQ=255)
                       gtgtGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                           >  1:1067440/1‑133 (MQ=255)
                       gtgtGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTc                                                                                                    >  1:495239‑M1/1‑138 (MQ=255)
                       gtgtGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTc                                                                                                    >  1:78729‑M1/1‑138 (MQ=255)
                         gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                                                  <  2:78729‑M1/139‑4 (MQ=255)
                         gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                                                  <  2:495239‑M1/139‑4 (MQ=255)
                                  aCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaa                                                                                         <  2:16079‑M1/139‑13 (MQ=255)
                                              agcagGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattac                                                                             <  2:250439‑M1/139‑25 (MQ=255)
                                                   ggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGTCATGTTcaccacaccaaagagatacattacggtcg                                                                        >  2:299350‑M1/1‑110 (MQ=255)
                                                     tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAATATCGGGATCCAGTTTCGCGGTAAAAATTAGTCAAGCAGAGCTCCGCCATGTTcaccacaccgagaagataaattacgttcgac                                                                      >  2:1197041‑M1/1‑108 (MQ=255)
                                                     tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccac                                                                                               <  2:899467‑M1/114‑7 (MQ=255)
                                                     tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccac                                                                                               >  1:899467‑M1/1‑108 (MQ=255)
                                                             aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                <  2:1315377‑M1/137‑38 (MQ=255)
                                                             aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                >  1:1315377‑M1/1‑100 (MQ=255)
                                                             aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                                           <  2:429296‑M1/110‑11 (MQ=255)
                                                             aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                                           >  1:429296‑M1/1‑100 (MQ=255)
                                                               cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccg                                                            >  1:1154569‑M1/1‑98 (MQ=255)
                                                                    tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                       <  1:299350‑M1/139‑47 (MQ=255)
                                                                         aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCTATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                >  1:321862‑M1/1‑88 (MQ=255)
                                                                         aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                <  2:321862‑M1/125‑38 (MQ=255)
                                                                         aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaa                                                                                         >  2:495035‑M1/1‑88 (MQ=255)
                                                                         aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaa                                                                                         <  1:495035‑M1/100‑13 (MQ=255)
                                                                            tGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggt                                                                          <  2:1230775‑M1/112‑28 (MQ=255)
                                                                            tGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggt                                                                          >  1:1230775‑M1/1‑85 (MQ=255)
                                                                                          atCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                           >  2:973241/1‑66 (MQ=255)
                                                                                          atCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                           <  1:973241/66‑1 (MQ=255)
                                                                                             cgcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacg                              <  2:982693‑M1/139‑72 (MQ=255)
                                                                                              gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                            >  1:88782/1‑45 (MQ=255)
                                                                                              gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                            <  2:88782/45‑1 (MQ=255)
                                                                                              gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacgg                             >  1:943807‑M1/1‑67 (MQ=255)
                                                                                              gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcacc                                                                                                 >  2:414388‑M1/1‑67 (MQ=255)
                                                                                              gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcacc                                                                                                 <  1:414388‑M1/71‑5 (MQ=255)
                                                                                               cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                            <  1:1399568/44‑1 (MQ=255)
                                                                                               cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                            >  2:1399568/1‑44 (MQ=255)
                                                                                                 tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaa                                                                                         <  1:1263786‑M1/76‑13 (MQ=255)
                                                                                                 tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaa                                                                                         >  2:1263786‑M1/1‑64 (MQ=255)
                                                                                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCt                                                                                                                <  2:1056763/52‑1 (MQ=255)
                                                                                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCt                                                                                                                >  1:1056763/1‑52 (MQ=255)
                                                                                                       cTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaa                    >  2:913833‑M1/1‑58 (MQ=255)
                                                                                                               cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGttcacc                                                                                                 >  1:752221/1‑50 (MQ=255)
                                                                                                               cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGttcacc                                                                                                 <  2:752221/54‑5 (MQ=255)
                                                                                                                         ttttGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccag  <  2:1237011‑M1/139‑100 (MQ=255)
                                                                                                                         ttttGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccag  <  1:655809‑M1/139‑100 (MQ=255)
                                                                                                                            tGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacg                              >  2:730009‑M1/1‑37 (MQ=255)
                                                                                                                            tGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacg                              <  1:730009‑M1/108‑72 (MQ=255)
                                                                                                                                       |                                                                                                                             
GCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTC  >  NC_000913/3862362‑3862622

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