Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,639,808 A→G T245A (ACG→GCG) 
*147Q (TAA→CAA) 
ugpQ →
yhhA ←
glycerophosphodiester phosphodiesterase, cytosolic
conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,639,8080AG100.0% 268.1 / NA 75T245A (ACG→GCG) 
*147Q (TAA→CAA) 
ugpQ
yhhA
glycerophosphodiester phosphodiesterase, cytosolic
conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/36);  total (39/36)

CAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGT  >  minE/2639740‑2639873
                                                                    |                                                                 
cagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCg                                                                 >  1:2659390/1‑71 (MQ=255)
cagAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCg                                                                 >  1:3008712/1‑71 (MQ=255)
 agagTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCgg                                                                <  1:78075/71‑1 (MQ=255)
  gagTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGc                                                               >  1:1919425/1‑71 (MQ=255)
    gTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGcc                                                              >  1:936476/1‑70 (MQ=255)
       gctgcGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCaa                                                           >  1:424424/1‑70 (MQ=255)
       gctgcGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCaa                                                           <  1:1900269/70‑1 (MQ=255)
        ctgcGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATa                                                         <  1:657370/71‑1 (MQ=255)
           cgcTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGtt                                                      >  1:259850/1‑71 (MQ=255)
           cgcTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGtt                                                      >  1:505414/1‑71 (MQ=255)
             cTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGtttt                                                    >  1:2536317/1‑71 (MQ=255)
                gggTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCg                                                                 >  1:2114974/1‑55 (MQ=255)
                gggTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTc                                                   <  1:2925719/69‑1 (MQ=255)
                gggTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCa                                                  <  1:104149/70‑1 (MQ=255)
                gggTGTGGATTGCATCTGTACCGACGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCaa                                                 <  1:1885567/71‑1 (MQ=255)
                    gtgGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTCCAACgg                                              <  1:395348/70‑1 (MQ=255)
                      ggATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAAt                                           >  1:524122/1‑71 (MQ=255)
                      ggATTGCATCTGTACCGATGCGATTGACGTGATTGGTCAGAACTTTGCGGCCCAATAGTTTTCAACGGAAt                                           >  1:1336645/1‑71 (MQ=255)
                       gATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATg                                          >  1:1753928/1‑71 (MQ=255)
                       gATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATg                                          <  1:2787598/71‑1 (MQ=255)
                        aTTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACTGAATGt                                         <  1:391434/71‑1 (MQ=255)
                        aTTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGt                                         <  1:2806839/71‑1 (MQ=255)
                        aTTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGt                                         <  1:2092011/71‑1 (MQ=255)
                         ttGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGtttt                                                    >  1:1753297/1‑59 (MQ=255)
                           gCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAg                                      >  1:722914/1‑71 (MQ=255)
                             aTCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGc                                    <  1:50785/71‑1 (MQ=255)
                             aTCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGc                                    <  1:2155001/71‑1 (MQ=255)
                              tCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAg                                                        <  1:2860589/50‑1 (MQ=255)
                              tCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAg                                                        <  1:2203590/50‑1 (MQ=255)
                               cTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAg                                                        <  1:93531/49‑1 (MQ=255)
                               cTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCt                                   >  1:115774/1‑70 (MQ=255)
                                tGTACCGATGCGATTGACGTGATTGGTCGGAACTTTGCGGCCCAATAGtt                                                      <  1:1410328/50‑1 (MQ=255)
                                 gTACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGaa                                            >  1:2840847/1‑59 (MQ=255)
                                  tACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTc                                                   >  1:2460527/1‑51 (MQ=255)
                                  tACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTgcg                                <  1:457015/70‑1 (MQ=255)
                                  tACCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGGTTTCAACGGAATGTCAGGCTgc                                 <  1:1486216/69‑1 (MQ=255)
                                   aCCGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTgcggc                              >  1:1170940/1‑71 (MQ=255)
                                     cGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGt                            >  1:2387003/1‑71 (MQ=255)
                                     cGATGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGt                            <  1:1444159/71‑1 (MQ=255)
                                       aTGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACgg                                              <  1:2429918/51‑1 (MQ=255)
                                       aTGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTc                                        <  1:2967247/57‑1 (MQ=255)
                                        tGCGATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCt                         >  1:1653681/1‑71 (MQ=255)
                                           gATTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTg                        <  1:1665846/69‑1 (MQ=255)
                                            aTTGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGc                     >  1:948735/1‑71 (MQ=255)
                                             ttGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                    >  1:1070928/1‑71 (MQ=255)
                                              tGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTg                           <  1:2403407/63‑1 (MQ=255)
                                              tGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTgg                       >  1:2154289/1‑67 (MQ=255)
                                              tGACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCt                    >  1:2067352/1‑70 (MQ=255)
                                               gACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCaa                                                 >  1:2039510/1‑40 (MQ=255)
                                               gACGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGa                  >  1:551427/1‑71 (MQ=255)
                                                aCGTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGCCTGa                  <  1:1878728/70‑1 (MQ=255)
                                                  gTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTTCTGGGCTGAtt                >  1:3019586/1‑70 (MQ=255)
                                                  gTGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTc               <  1:1413170/71‑1 (MQ=255)
                                                   tGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTc                                        <  1:407681/45‑1 (MQ=255)
                                                   tGATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCa              <  1:2724096/71‑1 (MQ=255)
                                                    gATTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAg             >  1:929849/1‑71 (MQ=255)
                                                     aTTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCa              >  1:2793536/1‑69 (MQ=255)
                                                     aTTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAg             <  1:579988/70‑1 (MQ=255)
                                                     aTTGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGc            <  1:2422553/71‑1 (MQ=255)
                                                       tGGTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCAt          >  1:320498/1‑71 (MQ=255)
                                                        ggTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCAt          <  1:1056191/70‑1 (MQ=255)
                                                        ggTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCAt          <  1:184746/70‑1 (MQ=255)
                                                         gTCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATg         <  1:514641/70‑1 (MQ=255)
                                                          tCCGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTgcg                                >  1:324306/1‑46 (MQ=255)
                                                           ccGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCg      >  1:805732/1‑71 (MQ=255)
                                                            cGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGc     >  1:516019/1‑71 (MQ=255)
                                                            cGAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGc     >  1:1746941/1‑71 (MQ=255)
                                                             gAACTTTGCGGCCCAATAGTTTTCAACGGAATGTCa                                       <  1:852640/36‑1 (MQ=255)
                                                              aaCTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTgaggc                              >  1:2058063/1‑44 (MQ=38)
                                                              aaCTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCg   >  1:2867108/1‑71 (MQ=255)
                                                               aCTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGcc    >  1:1205541/1‑69 (MQ=255)
                                                               aCTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGt  >  1:620028/1‑71 (MQ=255)
                                                               aCTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGt  >  1:768711/1‑71 (MQ=255)
                                                                cTTTGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTg                        <  1:1447858/48‑1 (MQ=255)
                                                                 tttGCGGCCCAATAGTTTTCAACGGAATGTCAGGCTgc                                 <  1:798498/38‑1 (MQ=255)
                                                                    |                                                                 
CAGAGTTGCTGCGCTGGGGTGTGGATTGCATCTGTACCGATGCGATTGACGTGATTGGTCCGAACTTTACGGCCCAATAGTTTTCAACGGAATGTCAGGCTGCGGCGTGCTGGGCTGATTCAGCATGTCGCCGT  >  minE/2639740‑2639873

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