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.96.8% 240.0 / 1.9 62pseudogene (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);  major base . (30/30);  minor base A (1/1);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

GTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATG  >  NC_000913/3862385‑3862650
                                                                                                                                   |                                                                                                                                      
gtgtGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTc                                                                                                                                >  2:1118209‑M1/1‑138 (MQ=255)
 tgtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTca                                                                                                                               <  2:577529‑M1/139‑3 (MQ=255)
  gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCC‑ACCATGTTcac                                                                                                                              >  1:825187‑M1/1‑136 (MQ=255)
     cGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccac                                                                                                                           >  1:1120119‑M1/1‑133 (MQ=255)
       gAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTca                                                                                                                               >  2:228314‑M1/1‑131 (MQ=255)
       gAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTca                                                                                                                               <  1:228314‑M1/133‑3 (MQ=255)
         aaaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacacca                                                                                                                       <  1:1021180‑M1/139‑11 (MQ=255)
           aCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaa                                                                                                                     <  1:1118209‑M1/139‑13 (MQ=255)
           aCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaa                                                                                                                     >  2:331862‑M1/1‑127 (MQ=255)
           aCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaa                                                                                                                     <  1:331862‑M1/139‑13 (MQ=255)
             cACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaaga                                                                                                                   <  1:545518‑M1/139‑15 (MQ=255)
                      cagcagGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacatta                                                                                                          >  1:824431‑M1/1‑116 (MQ=255)
                              tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgac                                                                                                  >  2:1294682‑M1/1‑108 (MQ=255)
                                 aGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacatt                                                                                               >  1:409099‑M1/1‑105 (MQ=255)
                                   cacaACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcacca                                                                                                                            <  1:596659‑M1/108‑6 (MQ=255)
                                   cacaACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcacca                                                                                                                            >  2:596659‑M1/1‑103 (MQ=255)
                                      aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggc                                                                                          >  1:925675‑M1/1‑100 (MQ=255)
                                      aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggc                                                                                          >  1:1187209‑M1/1‑100 (MQ=255)
                                        cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGtt                                                                                                                                 >  1:1313394/1‑98 (MQ=255)
                                        cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGtt                                                                                                                                 <  2:1313394/98‑1 (MQ=255)
                                        cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccg                                                                                        <  2:925675‑M1/139‑42 (MQ=255)
                                             tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                                                   <  1:1022153‑M1/139‑47 (MQ=255)
                                               gCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcac                                                                                 <  1:1311756‑M1/139‑49 (MQ=255)
                                                  aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                                              >  2:479397‑M1/1‑88 (MQ=255)
                                                    ttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgt                                                                            <  2:694571‑M1/139‑54 (MQ=255)
                                                                 atatCCGCGTGGTTGCTGCACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtagcgcaaacag                                                               <  1:882034‑M1/139‑67 (MQ=255)
                                                                       gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacacca                                                                                                                       >  1:40344‑M1/1‑67 (MQ=255)
                                                                       gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacacca                                                                                                                       <  2:40344‑M1/77‑11 (MQ=255)
                                                                        cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggt                                                        >  1:1269330‑M1/1‑66 (MQ=255)
                                                                        cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                                            <  1:772444‑M1/103‑38 (MQ=255)
                                                                        cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                                            >  2:772444‑M1/1‑66 (MQ=255)
                                                                         gTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagaga                                                                                                                 <  2:75985‑M1/81‑17 (MQ=255)
                                                                         gTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagaga                                                                                                                 >  1:75985‑M1/1‑65 (MQ=255)
                                                                          tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtga                                                      >  2:37992‑M1/1‑64 (MQ=255)
                                                                           ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                                     <  2:1269330‑M1/139‑77 (MQ=255)
                                                                            gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                                     >  1:1106329‑M1/1‑62 (MQ=255)
                                                                            gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                                     <  2:1106329‑M1/138‑77 (MQ=255)
                                                                            gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                                                   >  2:165856‑M1/1‑62 (MQ=255)
                                                                            gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                                                   <  1:165856‑M1/108‑47 (MQ=255)
                                                                             ttGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaat                                                   <  1:479397‑M1/139‑79 (MQ=255)
                                                                             ttGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcacacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatg                                                   >  1:800138‑M1/1‑61 (MQ=255)
                                                                                    aCATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCC‑ACCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacagg                                            <  2:825187‑M1/139‑86 (MQ=255)
                                                                                       tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgta                                                                         >  1:374486‑M1/1‑51 (MQ=255)
                                                                                       tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgta                                                                         <  2:374486‑M1/107‑57 (MQ=255)
                                                                                        cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaagg                                        >  2:1133398‑M1/1‑50 (MQ=255)
                                                                                                   tttGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagc                             <  1:32861‑M1/139‑101 (MQ=255)
                                                                                                       cGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacg                                                          >  2:735176‑M1/1‑35 (MQ=255)
                                                                                                       cGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacg                                                          <  1:735176‑M1/106‑72 (MQ=255)
                                                                                                       cGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgt                                                                            <  2:1260868‑M1/88‑54 (MQ=255)
                                                                                                       cGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgt                                                                            >  1:1260868‑M1/1‑35 (MQ=255)
                                                                                                            aaaCTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgatgccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggta                    >  2:1206227‑M1/1‑30 (MQ=255)
                                                                                                            aaaCTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacg                                                          <  2:169408‑M1/101‑72 (MQ=255)
                                                                                                            aaaCTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacg                                                          >  1:169408‑M1/1‑30 (MQ=255)
                                                                                                             aaCTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaa                   <  2:1120119‑M1/139‑111 (MQ=255)
                                                                                                               cttagTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaatt                 <  1:1206227‑M1/137‑113 (MQ=255)
                                                                                                                  aGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaattccc              <  1:1133398‑M1/139‑116 (MQ=255)
                                                                                                                    tCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtg                                                       >  1:1003528‑M1/1‑22 (MQ=255)
                                                                                                                    tCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtg                                                       <  2:1003528‑M1/96‑75 (MQ=255)
                                                                                                                        tCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaattcccaccagc        >  2:403080‑M1/1‑18 (MQ=255)
                                                                                                                          aGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattcaag                              >  1:830953‑M1/1‑16 (MQ=255)
                                                                                                                          aGACCTCCG‑CCATGTTcaccacaccaaagagatacatgacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccag                              <  2:830953‑M1/115‑100 (MQ=255)
                                                                                                                              cTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacatggaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacagtaaggtaaattccagcggttcggtaattcccaccagcatggcg  >  2:760575‑M1/1‑12 (MQ=255)
                                                                                                                                   |                                                                                                                                      
GTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATG  >  NC_000913/3862385‑3862650

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