Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,039,960 Δ1 bp intergenic (+42/‑36) yihD → / → yihE hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,039,9590A.100.0% 441.4 / NA 96intergenic (+41/‑37)yihD/yihEhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (48/48);  total (48/48)

Only 100 of 104 total aligned reads displayed.

GCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACA  >  NC_000913/4039885‑4040030
                                                                          |                                                                       
gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                        >  1:4298333/1‑75 (MQ=255)
gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                        <  1:3656460/76‑2 (MQ=255)
gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                        >  1:3014260/1‑75 (MQ=255)
gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                        >  1:395855/1‑75 (MQ=255)
 cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAGGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatc                                                                       <  1:3839807/76‑3 (MQ=255)
 cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatc                                                                       >  2:1016063/1‑74 (MQ=255)
 cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatc                                                                       >  1:1685894/1‑74 (MQ=255)
  tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcc                                                                      >  1:3427878/1‑73 (MQ=255)
     tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTc                                                                  <  2:2204755/76‑1 (MQ=255)
     tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTc                                                                  >  2:2962113/1‑76 (MQ=255)
     tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTc                                                                  <  1:1044327/76‑1 (MQ=255)
     tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTc                                                                  <  2:4591041/76‑1 (MQ=255)
      aCGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCg                                                                 <  2:2840475/76‑1 (MQ=255)
       cgcgCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGt                                                                <  1:1537777/76‑1 (MQ=255)
       cgcgCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGt                                                                >  1:1877627/1‑76 (MQ=255)
           cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCAtt                                                            >  2:2514625/1‑76 (MQ=255)
            tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTc                                                           <  2:1981075/76‑1 (MQ=255)
            tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTc                                                           >  1:3662212/1‑76 (MQ=255)
             cgcgGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCg                                                          >  1:185214/1‑76 (MQ=255)
              gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGt                                                         >  2:3915650/1‑76 (MQ=255)
              gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGt                                                         >  1:3940662/1‑76 (MQ=255)
               cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTa                                                        >  2:1072247/1‑76 (MQ=255)
               cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTa                                                        >  2:1044868/1‑76 (MQ=255)
               cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTa                                                        >  2:2883179/1‑76 (MQ=255)
               cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTa                                                        >  1:281429/1‑76 (MQ=255)
               cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTa                                                        >  2:1798641/1‑76 (MQ=255)
               cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATACCCTTCGTCATTCGTa                                                        <  2:4162254/76‑1 (MQ=255)
                ggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTaa                                                       >  2:1540515/1‑76 (MQ=255)
                ggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTaa                                                       <  2:4641930/76‑1 (MQ=255)
                ggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTaa                                                       >  1:2660661/1‑76 (MQ=255)
                  cgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAAtt                                                     >  1:326759/1‑75 (MQ=255)
                  cgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATg                                                     <  2:1776101/76‑1 (MQ=255)
                   gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGt                                                    <  2:1997652/76‑1 (MQ=255)
                    tAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtt                                                   <  2:2799924/76‑1 (MQ=255)
                    tAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATAGTGAAATGATATCCTTCGTCATTCGTAATGtt                                                   >  1:4266557/1‑76 (MQ=255)
                     aaTTAAAGAGTAAAAGCTTGTCAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGttt                                                  <  2:1283032/76‑1 (MQ=255)
                     aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGttt                                                  <  1:1129656/76‑1 (MQ=255)
                      aTTACAGAGTAAAAGCTTGTAAGCGGCGCCACCACACTCATCGTGACATGATATCCTTCGTCATTCGTACTGtttt                                                 >  2:3002332/1‑76 (MQ=255)
                      aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                 >  2:2506194/1‑76 (MQ=255)
                      aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                 <  1:3260377/76‑1 (MQ=255)
                      aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                 >  2:2992231/1‑76 (MQ=255)
                         aaaGCGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGCAATTCTACATGTTTTCCg                                              >  2:3637655/1‑76 (MQ=255)
                         aaaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCg                                              >  2:3633301/1‑76 (MQ=255)
                          aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCgg                                             <  2:1077360/76‑1 (MQ=255)
                          aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCgg                                             >  1:3604308/1‑76 (MQ=255)
                          aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGCTTTCCgg                                             >  1:158452/1‑76 (MQ=255)
                            gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgat                                           <  1:4744278/76‑1 (MQ=255)
                              gTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatga                                         <  1:3915650/76‑1 (MQ=255)
                              gTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatga                                         <  1:1589570/76‑1 (MQ=255)
                               tAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatgat                                        <  2:3239526/76‑1 (MQ=255)
                               gaaaaGCGTGTGAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTGTTCCGgatgat                                        <  2:193588/75‑1 (MQ=255)
                                  aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                     >  2:2803965/1‑76 (MQ=255)
                                  aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                     >  1:111525/1‑76 (MQ=255)
                                   aGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATgggg                                    >  2:134351/1‑75 (MQ=255)
                                   aGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGa                                    <  1:3673281/76‑1 (MQ=255)
                                   aGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGa                                    >  1:4399400/1‑76 (MQ=255)
                                    gCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGAt                                   <  2:809579/76‑1 (MQ=255)
                                     ctggTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                  <  2:4322040/73‑1 (MQ=255)
                                     cTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                  <  2:2400671/76‑1 (MQ=255)
                                     cTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                  <  2:2838886/76‑1 (MQ=255)
                                     cTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                  <  2:1964662/76‑1 (MQ=255)
                                      ttGTAATCGGCGCCACCAAAAGCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                 <  2:2193763/76‑1 (MQ=255)
                                      ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                 <  2:3081399/76‑1 (MQ=255)
                                       tGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                <  1:1248409/76‑1 (MQ=255)
                                         tAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                              <  1:3026223/76‑1 (MQ=255)
                                          aaGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                             <  1:3273012/76‑1 (MQ=255)
                                          aaGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                             >  2:3572745/1‑76 (MQ=255)
                                          aaGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                             <  2:2056384/76‑1 (MQ=255)
                                           aGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAAc                            <  1:737849/76‑1 (MQ=255)
                                              ggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgc                         >  1:61815/1‑76 (MQ=255)
                                              ggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgc                         >  1:2620304/1‑76 (MQ=255)
                                              ggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgc                         <  1:271431/76‑1 (MQ=255)
                                              ggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTAGTGATGTTTTCCGGATGAGGTGATCAATAAcctc                         >  1:95918/1‑73 (MQ=255)
                                                cgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                       >  1:2168166/1‑76 (MQ=255)
                                                cgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                       >  2:3761025/1‑76 (MQ=255)
                                                cgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                       >  1:3138824/1‑76 (MQ=255)
                                                 gcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCt                      >  1:451293/1‑76 (MQ=255)
                                                   caccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                    <  1:1134157/76‑1 (MQ=255)
                                                   caccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                    <  1:2880199/76‑1 (MQ=255)
                                                   caccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                    <  2:3671309/76‑1 (MQ=255)
                                                    accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtttt                   >  2:1487624/1‑76 (MQ=255)
                                                       aaaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                <  2:2794308/76‑1 (MQ=255)
                                                        aaaTCATCGTGAAATGATATCCTGCGGCATTCGGAATGGGTGCCGGATGAGGGGATGAATAACAGAGATTGTACtg               >  1:1429529/1‑75 (MQ=255)
                                                          atcatcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTc             >  1:1275808/1‑76 (MQ=255)
                                                           tcatcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTcc            <  2:4252253/76‑1 (MQ=255)
                                                           tcatcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTcc            <  1:1414595/76‑1 (MQ=255)
                                                            catcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCa           >  2:843767/1‑76 (MQ=255)
                                                             atcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAg          <  1:2605149/76‑1 (MQ=255)
                                                             atcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAg          >  1:2245458/1‑76 (MQ=255)
                                                               cGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGac        >  1:1162334/1‑76 (MQ=255)
                                                               cGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGac        >  1:4817035/1‑76 (MQ=255)
                                                               cGTGAAATGATAGCCTTCGTCATTCGTAATGTTTTCCGGATGATGGTAGTCATACAAGCGCTTTTACGTTCCAgcg        >  1:1167736/1‑74 (MQ=255)
                                                                gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca       <  2:2895969/76‑1 (MQ=255)
                                                                gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca       <  2:964708/76‑1 (MQ=255)
                                                                 tGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGacac      <  2:3172430/76‑1 (MQ=255)
                                                                 tGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGacac      <  1:2956508/76‑1 (MQ=255)
                                                                  gAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACt     >  1:244431/1‑76 (MQ=255)
                                                                    aaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac   <  2:326759/76‑1 (MQ=255)
                                                                    aaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac   <  1:3003822/76‑1 (MQ=255)
                                                                    aaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac   >  1:4348042/1‑76 (MQ=255)
                                                                          |                                                                       
GCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACA  >  NC_000913/4039885‑4040030

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