Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,873,099 2 bp→CT coding (238‑239/324 nt) ygbE ← DUF3561 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,873,0990TC100.0% 180.1 / NA 53D80G (GAC→GGC) ygbEDUF3561 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/23);  total (30/23)
*NC_0009132,873,1000CT100.0% 195.5 / NA 53D80N (GAC→AAC) ygbEDUF3561 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/23);  total (30/23)

TAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGTCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGG  >  NC_000913/2872957‑2873237
                                                                                                                                              ||                                                                                                                                         
tAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCaat                                                                                                                                    >  2:528959/1‑151 (MQ=255)
 aaTCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCaatg                                                                                                                                   >  2:734532/1‑151 (MQ=255)
      cATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCaatgaatg                                                                                                                               <  2:12033/150‑1 (MQ=255)
      cATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGc                                                                                                                              <  2:722959/151‑1 (MQ=255)
         tattaGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAg                                                                                                                            <  2:82749/150‑1 (MQ=255)
             aGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCa                                                                                                                        <  1:276148/150‑1 (MQ=255)
                tACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATcc                                                                                                                    >  1:270094/1‑151 (MQ=255)
                 aCGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATccc                                                                                                                   >  2:256783/1‑151 (MQ=255)
                       tggaTTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCcaccac                                                                                                             >  2:516275/4‑151 (MQ=255)
                         aaTTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACt                                                                                                            >  1:168539/1‑150 (MQ=255)
                          aTTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTAc                                                                                                          >  1:422336/1‑151 (MQ=255)
                            tATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACcg                                                                                                        >  1:34265/1‑151 (MQ=255)
                            tATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACcg                                                                                                        >  1:832347/1‑151 (MQ=255)
                               gAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACcgcg                                                                                                      >  2:763614/1‑150 (MQ=255)
                                aaCAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGAc                                                                                                    >  2:35119/1‑151 (MQ=255)
                                 aCAGAGGAGACTGGAAAGTACCAATTAGCCCAGTAGCCACATAAACAGAGCGCCCAACATAATGCCTACTGTCACCAGGGAGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAACGCAGCGCAATCCCCACCACTACCGCGACg                                                                                                   >  2:609213/1‑151 (MQ=255)
                                   agagGAGACGGGAAAGTACACATTAGCCCAGTCGCCACATAAAGAGCGCGCCAAACATAATGAATACTGTCACCAGGGTGAAAACCATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACggg                                                                                                 >  2:254695/1‑151 (MQ=255)
                                        agaCGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCa                                                                                            <  2:168515/151‑1 (MQ=255)
                                            gggAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAgcgc                                                                                        >  2:212436/1‑151 (MQ=255)
                                               aaaGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCag                                                                                     <  1:219076/151‑1 (MQ=255)
                                               aaaGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCag                                                                                     <  1:212462/151‑1 (MQ=255)
                                                       aaTTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCagaaagaa                                                                               >  2:475052/1‑149 (MQ=255)
                                                          tAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGc                                                                           <  1:3566/150‑1 (MQ=255)
                                                               cAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGgtg                                                                     <  2:517587/151‑1 (MQ=255)
                                                                aGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGgtgt                                                                    <  1:807078/151‑1 (MQ=255)
                                                                aGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGgtg                                                                     >  2:327509/1‑150 (MQ=255)
                                                                    gCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGAt                                                                <  1:510495/151‑1 (MQ=255)
                                                                       acgtAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGGAAACCATACTCTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAACCCCCACCACTACCGCGACGGGCATCAGCGCCGGAAAGAAAGGCAAGGTGTAGATaa                                                              <  2:653318/147‑1 (MQ=255)
                                                                       acaTAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATaa                                                              >  2:243158/1‑150 (MQ=255)
                                                                        caTAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATaa                                                              >  1:491533/1‑149 (MQ=255)
                                                                            aaCAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAaaga                                                           >  2:189004/1‑148 (MQ=255)
                                                                                  gcgcCAACCATAATGCCTACCGTCACCAGGGTGAAATCAATACTGTAGCGTAGCTTTCCGCTCATCAATGTATGCAGCGCATTCCCCACCACTACCGCAACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGTACAGCg                                                  <  2:648586/151‑1 (MQ=255)
                                                                                   cgcCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCgt                                                 <  1:580920/151‑1 (MQ=255)
                                                                                     ccAAACATAATGCCTACTGTCACCAGGGTGAAAACAATCCTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCAGCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCg                                                  <  2:547945/148‑1 (MQ=255)
                                                                                       aaaCATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGtt                                              >  1:580432/1‑150 (MQ=255)
                                                                                         aCATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTaga                                           <  1:146365/151‑1 (MQ=255)
                                                                                          cATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTagag                                          <  1:174980/151‑1 (MQ=255)
                                                                                             aaTGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGcc                                        >  2:885380/1‑150 (MQ=255)
                                                                                              aTGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTT‑GAGCCa                                       >  2:745202/1‑149 (MQ=255)
                                                                                                    aCTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATc                                 <  1:435057/150‑1 (MQ=255)
                                                                                                    aCTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATc                                 >  1:302341/1‑150 (MQ=255)
                                                                                                    aCTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCa                                <  2:416901/151‑1 (MQ=255)
                                                                                                         cACCAGGGTGAAAACAATACTGTAGCTTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATc                           <  2:778379/151‑1 (MQ=255)
                                                                                                              gggTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCAt                      >  2:572837/1‑151 (MQ=255)
                                                                                                               ggTGAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATc                     >  1:311492/1‑151 (MQ=255)
                                                                                                                  gAAAACAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGcc                  >  2:525569/1‑151 (MQ=255)
                                                                                                                     aaCAATACTGTAGCGTAGCTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAgc               >  2:645006/1‑151 (MQ=255)
                                                                                                                              gtagcgtagcTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATaa       <  2:233720/150‑1 (MQ=255)
                                                                                                                                 gcgtagcTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGACCCATAAATCAACATCGGCATCGCCATCGCAAATAAc      <  2:291931/148‑1 (MQ=255)
                                                                                                                                    tagcTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGCAAGGCAAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCAATAAATCAACATCGGCATCGCCAGCGCAAATAACCAgg  >  2:466394/1‑149 (MQ=255)
                                                                                                                                      gcTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAAt         <  1:335324/140‑1 (MQ=255)
                                                                                                                                      gcTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAAt         >  2:335298/1‑140 (MQ=255)
                                                                                                                                       cTTTCCGCTCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAgg  >  1:868240/1‑146 (MQ=255)
                                                                                                                                              ||                                                                                                                                         
TAATCTCATTATTAGTTACGGTTTGAATTATGAACAGAGGAGACGGGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCTACTGTCACCAGGGTGAAAACAATACTGTAGCGTAGCTTTCCGTCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGG  >  NC_000913/2872957‑2873237

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