Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,196,508 C→A intergenic (+2563/‑1229) ftsQ → / → ylmE cell division protein FtsQ/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,196,5080CA96.7% 205.1 / 0.5 61intergenic (+2563/‑1229)ftsQ/ylmEcell division protein FtsQ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (24/35);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

ATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAG  >  CP000730/1196372‑1196643
                                                                                                                                        |                                                                                                                                       
aTTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGc                                                                                                                                     >  2:1003738/1‑141 (MQ=255)
  tttCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCaa                                                                                                                                   <  1:211134/141‑1 (MQ=255)
  tttCGGGACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCaa                                                                                                                                   <  2:580157/141‑1 (MQ=255)
    tCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAAAATCACATGGACGTAATAGAGGTAGCACTGGATTCGAAACAAGCGTAAATACTTATAGCAATg                                                                                                                                 >  1:208415/1‑141 (MQ=255)
     cGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGc                                                                                                                                >  2:15767/1‑141 (MQ=255)
      ggTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCa                                                                                                                               <  1:638421/141‑1 (MQ=255)
      ggTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCa                                                                                                                               >  2:314511/1‑141 (MQ=255)
        tACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAAc                                                                                                                             <  1:211602/141‑1 (MQ=255)
            gTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCt                                                                                                                         <  1:653952/141‑1 (MQ=255)
               attaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTaaa                                                                                                                      <  1:31132/141‑1 (MQ=255)
                      cTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAAt                                                                                                               >  1:339469/1‑141 (MQ=255)
                        gAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACCTCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCa                                                                                                             >  2:615303/1‑141 (MQ=255)
                        gAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCa                                                                                                             >  1:471229/1‑141 (MQ=255)
                        gAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCa                                                                                                             >  2:130525/1‑141 (MQ=255)
                             aCAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                                                                        <  1:600873/141‑1 (MQ=255)
                             aCAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                                                                        <  2:215375/141‑1 (MQ=255)
                             aCAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTAATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                                                                        <  2:373177/141‑1 (MQ=255)
                               aaGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACtt                                                                                                      >  1:767031/1‑141 (MQ=255)
                                aGATGAGATTGTTGTAACAGTTATTGCAACTGGTGTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTc                                                                                                     <  1:954603/141‑1 (MQ=255)
                                   tGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCaaa                                                                                                  <  1:381966/141‑1 (MQ=255)
                                    gagaTTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAt                                                                                                 >  1:519304/1‑141 (MQ=255)
                                    gagaTTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAt                                                                                                 <  1:320483/141‑1 (MQ=255)
                                     agaTTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAtt                                                                                                <  2:705328/141‑1 (MQ=255)
                                        ttgttgTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtcat                                                                                             >  1:410104/1‑141 (MQ=255)
                                         tgttgtAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtcatc                                                                                            <  2:223454/141‑1 (MQ=255)
                                           ttgtAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCaa                                                                                          <  2:666252/141‑1 (MQ=255)
                                              tAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATg                                                                                       >  2:339046/1‑141 (MQ=255)
                                                aCAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGca                                                                                     >  2:228216/1‑141 (MQ=255)
                                                    ttATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAg                                                                                 <  1:313179/141‑1 (MQ=255)
                                                       ttGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtc                                                                                               >  1:838110/1‑124 (MQ=255)
                                                       ttGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtc                                                                                               <  2:838110/124‑1 (MQ=255)
                                                        tGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                                                                             >  1:185438/1‑141 (MQ=255)
                                                            aCTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGAt                                                                         <  1:246939/141‑1 (MQ=255)
                                                                gTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgtg                                                                     <  1:350056/141‑1 (MQ=255)
                                                                gTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgtg                                                                     <  1:121599/141‑1 (MQ=255)
                                                                   ttGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTaa                                                                  <  2:301223/141‑1 (MQ=255)
                                                                    tgatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAg                                                                 >  1:86029/1‑141 (MQ=255)
                                                                     gatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGt                                                                <  2:689059/141‑1 (MQ=255)
                                                                     gatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGt                                                                <  1:570408/141‑1 (MQ=255)
                                                                          cAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAg                                                           >  1:792469/1‑141 (MQ=255)
                                                                              ccAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAaca                                                       <  1:130525/141‑1 (MQ=255)
                                                                                    tCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACa                                                 >  2:652430/1‑141 (MQ=255)
                                                                                    tCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACa                                                 >  2:492/1‑141 (MQ=255)
                                                                                    tCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACa                                                 >  2:579111/1‑141 (MQ=255)
                                                                                      acaTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAAc                                               >  1:141479/1‑141 (MQ=255)
                                                                                       aaTGGTCGTAAATCTGGGAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACt                                              <  2:682803/140‑1 (MQ=255)
                                                                                          ggTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTaaa                                           <  1:41056/141‑1 (MQ=255)
                                                                                              gTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAaagaag                                       >  2:976424/1‑141 (MQ=255)
                                                                                               tAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAagaaga                                      <  2:519304/141‑1 (MQ=255)
                                                                                                  aTCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAgatga                                   <  2:290772/141‑1 (MQ=255)
                                                                                                   tCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAgatgat                                  <  1:1003879/141‑1 (MQ=255)
                                                                                                           cACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAGCTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAg                          <  1:495630/141‑1 (MQ=255)
                                                                                                           cACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAg                          <  2:21846/141‑1 (MQ=255)
                                                                                                             cTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCt                        <  2:767031/141‑1 (MQ=255)
                                                                                                                      gAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGaa               >  1:108829/1‑141 (MQ=255)
                                                                                                                       aaCAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGaaa              <  1:623517/141‑1 (MQ=255)
                                                                                                                         cAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATa            <  1:832453/141‑1 (MQ=255)
                                                                                                                         cAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATa            <  1:466737/141‑1 (MQ=255)
                                                                                                                          aaGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATag           >  1:383445/1‑141 (MQ=255)
                                                                                                                           aGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATaga          >  1:421496/1‑141 (MQ=255)
                                                                                                                                   tACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACGCATACAACTAAAGAAGATGATATTCCAAGCTTCATTAGAAATAGAGaagaaag  <  2:637613/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
ATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAG  >  CP000730/1196372‑1196643

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