Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,533,126 A→G F20L (TTC→CTC)  yjhP ← KpLE2 phage‑like element; putative methyltransferase YjhP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,533,1260AG100.0% 190.0 / NA 54F20L (TTC→CTC) yjhPKpLE2 phage‑like element; putative methyltransferase YjhP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/29);  total (25/29)

ATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATC  >  NC_000913/4532990‑4533262
                                                                                                                                        |                                                                                                                                        
aTGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCgg                                                                                                                                       >  2:566296/1‑140 (MQ=255)
  gCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGt                                                                                                                                     >  2:57365/1‑140 (MQ=255)
       tCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATgg                                                                                                                                >  1:1115321/1‑140 (MQ=255)
         cccGTAATGCCATGATCCCTGGTCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGAt                                                                                                                              <  1:30164/140‑1 (MQ=255)
          ccGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATg                                                                                                                             <  1:178962/140‑1 (MQ=255)
            gTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATg                                                                                                                             <  2:224896/138‑1 (MQ=255)
            gTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATg                                                                                                                             >  1:224896/1‑138 (MQ=255)
                  ccATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTc                                                                                                                     >  1:905069/1‑140 (MQ=255)
                   cATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTTCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGtt                                                                                                                                    >  2:1084136/1‑124 (MQ=255)
                   cATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGtt                                                                                                                                    <  1:1084136/124‑1 (MQ=255)
                           cTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCt                                                                                                            <  1:335613/140‑1 (MQ=255)
                             ggCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGa                                                                                                          <  1:566296/140‑1 (MQ=255)
                              gCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGAt                                                                                                         <  2:381063/140‑1 (MQ=255)
                                   ggTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTaa                                                                                                    >  1:1126108/1‑140 (MQ=255)
                                   ggTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTaa                                                                                                    >  2:269987/1‑140 (MQ=255)
                                   ggTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTaa                                                                                                    >  2:978680/1‑140 (MQ=255)
                                   ggTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTaa                                                                                                    <  1:234334/140‑1 (MQ=255)
                                   ggTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTaa                                                                                                    >  2:882059/1‑140 (MQ=255)
                                   ggTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTaa                                                                                                    >  2:1114644/1‑140 (MQ=255)
                                       cAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAAt                                                                                                >  2:633900/1‑140 (MQ=255)
                                         gagCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATAc                                                                                              >  1:943725/1‑140 (MQ=255)
                                          agCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACg                                                                                             <  1:1008033/140‑1 (MQ=255)
                                            cATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTg                                                                                           >  2:1193648/1‑140 (MQ=255)
                                            cATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTg                                                                                           <  1:269987/140‑1 (MQ=255)
                                              tctcGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTggg                                                                                         <  2:349488/140‑1 (MQ=255)
                                               ctcGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGTTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGa                                                                                        >  2:586363/1‑140 (MQ=255)
                                               ctcGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGa                                                                                        >  2:103276/1‑140 (MQ=255)
                                                tcGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGTTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGat                                                                                       <  1:586363/140‑1 (MQ=255)
                                                tcGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGat                                                                                       <  1:1005829/140‑1 (MQ=255)
                                                tcGCCCGAGCAGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGat                                                                                       <  1:1133469/140‑1 (MQ=255)
                                                     cGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCa                                                                                  <  1:633900/140‑1 (MQ=255)
                                                              gCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAAt                                                                         <  2:1115321/140‑1 (MQ=255)
                                                                cGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTc                                                                       <  2:905069/140‑1 (MQ=255)
                                                                 gAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGAGTTCACTTTCGCTGATGGTAAAAATACGTGGGATACCCATTATTTAATTcc                                                                      <  2:900711/140‑1 (MQ=255)
                                                                  aGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCt                                                                     >  2:251231/1‑140 (MQ=255)
                                                                      cAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTa                                                                 <  1:1114644/140‑1 (MQ=255)
                                                                             gCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGtt                                                          <  1:781634/140‑1 (MQ=255)
                                                                             gCGGGGGCTCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTGCTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGtt                                                          <  2:224784/140‑1 (MQ=255)
                                                                              cGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTa                                                         >  2:1190484/1‑140 (MQ=255)
                                                                              cGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTa                                                         >  1:191833/1‑140 (MQ=255)
                                                                                  tGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGcc                                                     <  1:1191217/140‑1 (MQ=255)
                                                                                    cccGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCtt                                                   >  2:681056/1‑140 (MQ=255)
                                                                                            cATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGa                                           <  1:251231/140‑1 (MQ=255)
                                                                                            cATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGa                                           <  1:1193648/140‑1 (MQ=255)
                                                                                                     caCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAg                                  >  1:310509/1‑140 (MQ=255)
                                                                                                        gcgACCCAGCGTGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGAtt                               <  1:681056/140‑1 (MQ=255)
                                                                                                                   tGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGt                    <  1:882059/140‑1 (MQ=255)
                                                                                                                   tGGCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGt                    <  1:978680/140‑1 (MQ=255)
                                                                                                                    ggCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTAtt                         >  1:1007560/1‑134 (MQ=255)
                                                                                                                    ggCGTACTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTAtt                         <  2:1007560/134‑1 (MQ=255)
                                                                                                                          cTTCTCTTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCt             >  1:398490/1‑140 (MQ=255)
                                                                                                                             ctctTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCg          >  1:1160931/1‑140 (MQ=255)
                                                                                                                               ctTCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGt        >  2:828728/1‑140 (MQ=255)
                                                                                                                                ttCGGTGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGtt       <  1:828728/140‑1 (MQ=255)
                                                                                                                                     tGAGCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtcatc  <  1:103276/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
ATGCCAGTCCCCGTAATGCCATGATCCCTGGCCCAGGTGCAGAGCATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATC  >  NC_000913/4532990‑4533262

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