Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,732,027 A→G 6.8% M716T (ATG→ACG)  clpB ← protein disaggregation chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,732,0270AG 6.8% 181.8 / 3.9 59M716T (ATG→ACG) clpBprotein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base G (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

TGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGAC  >  NC_000913/2731883‑2732174
                                                                                                                                                |                                                                                                                                                   
tGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAAt                                                                                                                                                 >  4:159306/1‑149 (MQ=255)
    aGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATgacg                                                                                                                                             >  2:552260/1‑149 (MQ=255)
     gACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATgacga                                                                                                                                            <  3:469557/149‑1 (MQ=255)
        caccacTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCg                                                                                                                                         >  1:166499/1‑149 (MQ=255)
                  tCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATACTCCACTTCACCGAAGCGTTCCTGAATCAGATCTGAACCAAGGTCAGAGGTCATAATCACGACCGTATTACGGaa                                                                                                                               <  4:350504/149‑1 (MQ=255)
                  tCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGaa                                                                                                                               <  2:189933/149‑1 (MQ=255)
                     aTACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTc                                                                                                                            >  3:338133/1‑149 (MQ=255)
                     aTACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTc                                                                                                                            <  4:225095/149‑1 (MQ=255)
                               tGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTtctc                                                                                                                  >  4:378888/1‑149 (MQ=255)
                                gAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTcc                                                                                                                 <  2:220153/149‑1 (MQ=255)
                                gAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTcc                                                                                                                 >  2:260073/1‑149 (MQ=255)
                                 aaTTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATTTGCGCATAATCTAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGATAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTccc                                                                                                                <  4:348543/149‑1 (MQ=255)
                                  aTTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCGTAATGACGACCGTa                                                                                                                                       >  3:411805/1‑125 (MQ=255)
                                  aTTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCGTAATGACGACCGTa                                                                                                                                       <  4:522439/125‑1 (MQ=255)
                                  aTTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCGTAATGACGACCGTa                                                                                                                                       >  3:522439/1‑125 (MQ=255)
                                  aTTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCGTAATGACGACCGTa                                                                                                                                       <  4:411805/125‑1 (MQ=255)
                                    tCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTg                                                                                                             <  1:262060/149‑1 (MQ=255)
                                           ggaAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCgtca                                                                                                      >  2:309568/1‑149 (MQ=255)
                                            gcaGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCgtcag                                                                                                     <  3:461363/147‑1 (MQ=255)
                                                ttATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGtcagtcag                                                                                                 >  3:548145/1‑149 (MQ=255)
                                                 tATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGa                                                                                                <  3:55984/149‑1 (MQ=255)
                                                       ttACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCca                                                                                          >  4:498217/1‑149 (MQ=255)
                                                          ccACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCcatca                                                                                       >  3:189778/1‑149 (MQ=255)
                                                           cacaCCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCcatcat                                                                                      >  3:503440/1‑149 (MQ=255)
                                                             cacCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATaa                                                                                                                                                  <  4:293915/87‑1 (MQ=255)
                                                             cacCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATaa                                                                                                                                                  >  3:293915/1‑87 (MQ=255)
                                                              acCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCa                                                                                   >  4:145859/1‑149 (MQ=255)
                                                                 gAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTa                                                                                >  2:292963/1‑149 (MQ=255)
                                                                    cACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCt                                                                             <  2:249814/149‑1 (MQ=255)
                                                                         cctctTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCTGATCGGAACAGGGGTTAGATGTCATACTCACGACCGTATTACGGAAGGCGACCGCTCTCCCTTGCCCATCAGTCAGACCCCCCTCGTCCAGTACCAGCAAc                                                                        <  4:553297/148‑1 (MQ=255)
                                                                           tctTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAg                                                                      >  3:550056/1‑149 (MQ=255)
                                                                             tttCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGaa                                                                    <  3:159306/149‑1 (MQ=255)
                                                                               tCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATg                                                                  <  3:160144/149‑1 (MQ=255)
                                                                                 aTGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGtt                                                                >  1:126423/1‑149 (MQ=255)
                                                                                    tgCGCATATTCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGaa                                                             <  4:513610/149‑1 (MQ=255)
                                                                                    tgCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGGTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGaa                                                             <  2:341980/149‑1 (MQ=255)
                                                                                        cATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACa                                                         >  3:148749/1‑149 (MQ=255)
                                                                                         aTAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACAt                                                        >  1:450757/1‑149 (MQ=255)
                                                                                           aaTCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATcc                                                      <  4:338133/149‑1 (MQ=255)
                                                                                            aTCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCg                                                     <  3:140525/149‑1 (MQ=255)
                                                                                               cAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGAt                                                  >  1:378353/1‑149 (MQ=255)
                                                                                                 gTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGccc                                                                                           <  4:516398/106‑1 (MQ=255)
                                                                                                 gTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGccc                                                                                           >  3:516398/1‑106 (MQ=255)
                                                                                                   tCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATgcgc                                              >  3:214631/1‑149 (MQ=255)
                                                                                                    cACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCt                                             <  3:166281/149‑1 (MQ=255)
                                                                                                          aGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTcc                                       <  1:180930/149‑1 (MQ=255)
                                                                                                            cGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCAc                                     <  1:552260/149‑1 (MQ=255)
                                                                                                             gTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACt                                    >  4:411003/1‑149 (MQ=255)
                                                                                                                      tCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCcagc                           >  4:228466/1‑149 (MQ=255)
                                                                                                                      tCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCcagc                           >  4:319352/1‑149 (MQ=255)
                                                                                                                       cAGATCGGAACCGAGGTTAGAGGTCATAATGACGACTGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCAGATGCGCTTTTTCCACTTCATCcagca                          <  3:378888/149‑1 (MQ=255)
                                                                                                                           tCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGAt                      <  1:334231/149‑1 (MQ=255)
                                                                                                                           tCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGAt                      <  3:145859/149‑1 (MQ=255)
                                                                                                                                   gAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAAt              <  1:553948/149‑1 (MQ=255)
                                                                                                                                     ggTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATAc            <  2:450757/149‑1 (MQ=255)
                                                                                                                                      gTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACg           <  2:126423/149‑1 (MQ=255)
                                                                                                                                         agagGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGgac        <  1:424208/149‑1 (MQ=255)
                                                                                                                                             gTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGgacgacg    >  4:237390/1‑149 (MQ=255)
                                                                                                                                               cATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGgacgacgac  <  3:86313/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCAGCTCTTTCATGTGCGCATAATCCAGTTCACCGAAGCGTTCCTGAATCAGATCGGAACCGAGGTTAGAGGTCATAATGACGACCGTATTACGGAAGTCGACCGTTCTCCCTTGCCCGTCAGTCAGACGCCCATCATCCAGTACCTGCAACAGAATGTTGAAGACATCCGGATGCGCTTTTTCCACTTCATCCAGCAGGATGACGGAATACGGACGACGAC  >  NC_000913/2731883‑2732174

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