Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,619,913 A→G E104G (GAA→GGA)  deoB → phosphopentomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,619,9130AG100.0% 169.4 / NA 48E104G (GAA→GGA) deoBphosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/24);  total (24/24)

GCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGG  >  NC_000913/4619782‑4620048
                                                                                                                                    |                                                                                                                                       
gCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCCGTCTGGTCACTGGGGAATTGc                                                                                                                                   <  1:1091074/139‑1 (MQ=255)
  tGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCgg                                                                                                                                <  2:1317834/139‑1 (MQ=255)
    gCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGgtg                                                                                                                              >  1:431455/1‑139 (MQ=255)
          gCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCgg                                                                                                                        <  2:1232032/139‑1 (MQ=255)
            acacGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGtt                                                                                                                      <  2:801812/139‑1 (MQ=255)
                  aGGTTCTACCGGTTTCATTCCGGTGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGttt                                                                                                                >  2:928628/1‑139 (MQ=255)
                          ccGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGcc                                                                                                                                  <  1:649897/113‑1 (MQ=255)
                          ccGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGcc                                                                                                                                  >  2:649897/1‑113 (MQ=255)
                               ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCATTGGGGAATTGc                                                                                                                                   <  1:852613/107‑1 (MQ=255)
                               ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGc                                                                                                                                   >  2:852613/1‑107 (MQ=255)
                                    tCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTcc                                                                                              >  2:252911/1‑139 (MQ=255)
                                        gcggGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATc                                                                                          >  2:256112/1‑139 (MQ=255)
                                            gAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGa                                                                                      <  2:431455/139‑1 (MQ=255)
                                                      cAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTc                                                                            >  1:637305/1‑139 (MQ=255)
                                                       aaCGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAg                                                                                                             >  2:886435/1‑105 (MQ=255)
                                                       aaCGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAg                                                                                                             <  1:886435/105‑1 (MQ=255)
                                                           cTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTcc                                                                                                                           >  2:285860/1‑87 (MQ=255)
                                                           cTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTcc                                                                                                                           <  1:285860/87‑1 (MQ=255)
                                                                 ttATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgct                                                                 <  1:252911/139‑1 (MQ=255)
                                                                  tATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctg                                                                <  1:1364824/139‑1 (MQ=255)
                                                                   aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctgc                                                               <  2:1050151/139‑1 (MQ=255)
                                                                   aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctgc                                                               >  1:1248418/1‑139 (MQ=255)
                                                                        cgcgTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCgg                                                                                                                                >  1:173671/1‑69 (MQ=255)
                                                                        cgcgTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCgg                                                                                                                                <  2:173671/69‑1 (MQ=255)
                                                                          cgTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATaa                                                        >  1:176607/1‑139 (MQ=255)
                                                                            tACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAAc                                                      >  1:708175/1‑139 (MQ=255)
                                                                                    ggCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGaa                                              >  1:1258345/1‑139 (MQ=255)
                                                                                         aCGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACgcgc                                         >  1:304822/1‑139 (MQ=255)
                                                                                                gTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTg                                  <  1:928628/139‑1 (MQ=255)
                                                                                                      cGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGt                            <  1:836357/139‑1 (MQ=255)
                                                                                                      cGGTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGt                            <  2:637305/139‑1 (MQ=255)
                                                                                                       ggTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCgg                                                                                                                                <  1:1389863/38‑1 (MQ=255)
                                                                                                       ggTAAAGATA‑CCCCGTCTGGTCACTGGGGAATTGCCgg                                                                                                                                >  2:1389863/1‑38 (MQ=255)
                                                                                                                a‑CCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTa                  >  2:362261/1‑139 (MQ=255)
                                                                                                                   cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGt                                                                                                                  >  1:752742/1‑41 (MQ=255)
                                                                                                                   cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGt                                                                                                                  <  2:752742/41‑1 (MQ=255)
                                                                                                                   cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc                >  2:747153/1‑139 (MQ=255)
                                                                                                                   cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc                >  2:214579/1‑139 (MQ=255)
                                                                                                                   cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc                >  2:1289840/1‑139 (MQ=255)
                                                                                                                    ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAg                                                                                                             >  2:1233769/1‑45 (MQ=255)
                                                                                                                    ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAg                                                                                                             <  1:1233769/45‑1 (MQ=255)
                                                                                                                      gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGc             >  1:1268263/1‑139 (MQ=255)
                                                                                                                      gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGc             <  1:362261/139‑1 (MQ=255)
                                                                                                                        cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa           <  2:1268263/139‑1 (MQ=255)
                                                                                                                        cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa           <  2:1248418/139‑1 (MQ=255)
                                                                                                                              aCTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTc     <  1:256112/139‑1 (MQ=255)
                                                                                                                                tGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCg                     <  1:249175/121‑1 (MQ=255)
                                                                                                                                tGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCg                     >  2:249175/1‑121 (MQ=255)
                                                                                                                                 ggggAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa           >  2:1311189/1‑130 (MQ=255)
                                                                                                                                 ggggAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa           <  1:1311189/130‑1 (MQ=255)
                                                                                                                                 ggggAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCgg  <  2:176607/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
GCTGGCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATA‑CCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGG  >  NC_000913/4619782‑4620048

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