Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,784,184 T→C I223T (ATA→ACA)  yfjB → NAD kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,784,1840TC100.0% 307.2 / NA 85I223T (ATA→ACA) yfjBNAD kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (50/35);  total (50/35)

CCCCTCTCTGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGA  >  minE/1784116‑1784250
                                                                    |                                                                  
ccccTCTCTGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACaa                                                                  <  1:2426802/71‑1 (MQ=255)
  ccTCTCTGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACaaa                                                                 >  1:707083/1‑70 (MQ=255)
  ccTCTCTGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACaaa                                                                 >  1:1425334/1‑70 (MQ=255)
   ctctctGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACaaa                                                                 >  1:1653247/1‑69 (MQ=255)
   ctctctGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACaaa                                                                 >  1:1836045/1‑69 (MQ=255)
   ctctctGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACaaa                                                                 >  1:519774/1‑69 (MQ=255)
     ctctGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcag                                                              >  1:2178150/1‑70 (MQ=255)
      tctGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcag                                                              >  1:940879/1‑69 (MQ=255)
      tctGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcag                                                              >  1:2174114/1‑69 (MQ=255)
      tctGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAca                                                               >  1:1288398/1‑68 (MQ=255)
       ctGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcagcag                                                           >  1:1141572/1‑71 (MQ=255)
         ggATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcagcag                                                           >  1:1504691/1‑69 (MQ=255)
          gATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACTACCACTGGTCACAAAcagcagc                                                          <  1:1669384/69‑1 (MQ=255)
          gATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcagcagcag                                                        <  1:1606759/71‑1 (MQ=255)
          gATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcagcagcag                                                        >  1:2814414/1‑71 (MQ=255)
          gATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcagcag                                                           <  1:388511/68‑1 (MQ=255)
           aTGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcagcagcag                                                        >  1:209796/1‑70 (MQ=255)
            tGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcagca                                                            <  1:241168/65‑1 (MQ=255)
              cGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACg                                                    <  1:1006848/71‑1 (MQ=255)
               gATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcagcag                                                           <  1:920388/63‑1 (MQ=255)
               gATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGa                                                   <  1:1387348/71‑1 (MQ=255)
                aTTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACg                                                    <  1:2236655/69‑1 (MQ=255)
                aTTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGAt                                                  <  1:695476/71‑1 (MQ=255)
                  tACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATcc                                                >  1:2267239/1‑71 (MQ=255)
                   aCCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCg                                               >  1:3008803/1‑71 (MQ=255)
                     ccTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcagcag                                                           <  1:207788/57‑1 (MQ=255)
                       tGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAAcagcagca                                                         >  1:655401/1‑57 (MQ=255)
                       tGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTg                                           >  1:2672220/1‑71 (MQ=255)
                       tGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTg                                           >  1:1152325/1‑71 (MQ=255)
                         gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGc                                          <  1:2428010/70‑1 (MQ=255)
                         gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCg                                         <  1:2429539/71‑1 (MQ=255)
                         gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCg                                         <  1:1306804/71‑1 (MQ=255)
                          tGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGt                                        >  1:2943750/1‑71 (MQ=255)
                           gCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGtt                                       >  1:1764038/1‑71 (MQ=255)
                             ccATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGtttt                                     >  1:1521639/1‑71 (MQ=255)
                               aTGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATcc                                                <  1:2295863/58‑1 (MQ=255)
                               aTGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATcc                                                <  1:348023/58‑1 (MQ=255)
                               aTGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATcc                                                <  1:203622/58‑1 (MQ=255)
                               aTGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTc                                   >  1:2671272/1‑71 (MQ=255)
                               aTGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTc                                   >  1:1008307/1‑71 (MQ=255)
                               aTGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTc                                   <  1:1377227/71‑1 (MQ=255)
                                tGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCt                                            >  1:1129190/1‑61 (MQ=255)
                                tGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCg                                  >  1:1237937/1‑71 (MQ=255)
                                tGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCAAGATCCGTCTGCGTTTTTCg                                  >  1:1439273/1‑71 (MQ=255)
                                 gTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGt                                        <  1:2357918/64‑1 (MQ=255)
                                 gTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCg                                  >  1:1444403/1‑70 (MQ=255)
                                  ttCCCGCATACGTTGTCAGCACGCCCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCg                                  >  1:2268103/1‑69 (MQ=255)
                                     ccGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCg                             <  1:2805240/71‑1 (MQ=255)
                                     ccGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCg                             >  1:491125/1‑71 (MQ=255)
                                     ccGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCg                             >  1:239107/1‑71 (MQ=255)
                                      cGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTc                                   <  1:23901/64‑1 (MQ=255)
                                      cGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGc                            >  1:1584454/1‑71 (MQ=255)
                                         aTACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCg                          >  1:58333/1‑70 (MQ=255)
                                         aTACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCg                          >  1:1667382/1‑70 (MQ=255)
                                         aTACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGt                         <  1:423991/71‑1 (MQ=255)
                                         aTACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGt                         <  1:2814136/71‑1 (MQ=255)
                                             gatGTCAGCACGACCACTGGTCACAAACACCAGCAGCACGATCCGTCTGCGt                                        >  1:761874/3‑52 (MQ=38)
                                             gTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTaa                       >  1:2733170/1‑69 (MQ=255)
                                             gTTGACAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGt                                              >  1:1522992/1‑46 (MQ=39)
                                               tGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGAc                   >  1:353419/1‑71 (MQ=255)
                                               tGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATAGCCGTAACGAc                   >  1:992045/1‑71 (MQ=255)
                                                gTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGtt                                       >  1:2111875/1‑50 (MQ=255)
                                                gTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGtt                                       >  1:2403766/1‑50 (MQ=255)
                                                gTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGAcc                  >  1:1527539/1‑71 (MQ=255)
                                                 tCAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCt                 <  1:2417879/71‑1 (MQ=255)
                                                  cAGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTg                >  1:117794/1‑71 (MQ=255)
                                                   aGCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCg                                  >  1:2022433/1‑52 (MQ=255)
                                                    gCACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGt                         <  1:2268046/60‑1 (MQ=255)
                                                     cACGACCACTGGTCACAAACAGCAGCAGCACGATCCGt                                              >  1:1772052/1‑38 (MQ=255)
                                                      aCGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTaa                       >  1:622721/1‑60 (MQ=255)
                                                        gACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGtttt                                     >  1:1996538/1‑44 (MQ=255)
                                                        gACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGc                                 <  1:364223/48‑1 (MQ=255)
                                                        gACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGc                                 <  1:1688735/48‑1 (MQ=255)
                                                          ccACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAg        <  1:1355928/71‑1 (MQ=255)
                                                           cACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt       <  1:1301044/71‑1 (MQ=255)
                                                           cACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt       <  1:193756/71‑1 (MQ=255)
                                                            aCTGGTCACAAACAGCAGCGGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGtt      >  1:1050327/1‑71 (MQ=255)
                                                            aCTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGtt      >  1:317570/1‑71 (MQ=255)
                                                             cTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt       <  1:1116149/69‑1 (MQ=255)
                                                             cTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTg     >  1:1731242/1‑71 (MQ=255)
                                                              tGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGa              <  1:1477784/61‑1 (MQ=255)
                                                              tGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATc          <  1:2661059/65‑1 (MQ=255)
                                                              tGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGc    >  1:184301/1‑71 (MQ=255)
                                                                gTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGa  <  1:1629276/71‑1 (MQ=255)
                                                                gTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGa  <  1:903718/71‑1 (MQ=255)
                                                                    |                                                                  
CCCCTCTCTGGATGCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGA  >  minE/1784116‑1784250

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: