Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,731,822	T→C	K784K (AAA→AAG)	clpB ←	ClpB chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,731,822	0	T	C	100.0%	175.2 / NA	49	K784K (AAA→AAG)	clpB	ClpB chaperone
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/26);  total (23/26)

GCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGTTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCA  >  NC_000913/2731685‑2731955                                                                                                                                          |                                                                                                                                      gCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGctt                                                                                                                                     >  1:828940/1‑140 (MQ=255) gCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGAATTGCATAACCACGTTCTTCCAGACGctt                                                                                                                                     >  1:406420/1‑140 (MQ=255) gCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCACAAATGTGGATTGCAGAACCACGTTCTTCCAGACGctt                                                                                                                                     >  1:308033/1‑140 (MQ=255)   gacggggTTTTCGATATGCTGCTGAGTTGGACGTTTCAGAGGGCGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGt                                                                                                                                   <  1:370444/136‑1 (MQ=255)                aTCTGCTGCTGAATTGCACGTTTCAGAGGACGTGCACCATAGACCGAATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCaa                                                                                                                      <  1:631719/140‑1 (MQ=255)                   tgctgctgAATTGCACGTTTCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTTGATTTCATCACCACGTTCTTCCAGAAGCTTGTACAGACGTTTcctgtc                                                                                                                   >  2:370444/1‑135 (MQ=255)                     ctgctgAATTGCACGTTTCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGaa                                                                                                                 <  1:623941/140‑1 (MQ=255)                           aaTTGCACGTTTCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCACAACCACGTTCTTACAGCCGCTTGTACAGACATTTCAACTGAATCTgcg                                                                                                           >  2:302103/1‑140 (MQ=255)                              tGCACGTTTCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCGGACGCTTGTACAGACGTTTCAACTGAATCTGCGCaa                                                                                                        <  2:100437/140‑1 (MQ=255)                                 aCGTTTCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGCTTCAACTGAATCTGCGCAATCg                                                                                                     <  2:534867/140‑1 (MQ=255)                                     ttCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGc                                                                                                 >  2:66759/1‑140 (MQ=255)                                      tCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCa                                                                                                <  1:66759/140‑1 (MQ=255)                                        agagGACGTGCACCATAGACCGTATCGTAGCCGTTCTGGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAAt                                                                                              <  1:302103/140‑1 (MQ=255)                                          agGACGTGCACCATAGATCGGAACGTAACCGTTCTCGCTTAGCAGTTTAAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAAtgt                                                                                            <  2:308033/140‑1 (MQ=255)                                          agGACGTGAACCATAGTCCGGAACGTAACCATTCTAGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAAtgt                                                                                            <  2:406420/140‑1 (MQ=255)                                               gTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGt                                                                                       >  1:715983/1‑140 (MQ=255)                                                    ccATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCAcc                                                                                  <  1:574401/140‑1 (MQ=255)                                                      aTAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCAcc                                                                                  >  1:509807/1‑138 (MQ=255)                                                      aTAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCAcc                                                                                  <  2:509807/138‑1 (MQ=255)                                                      aTAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCa                                                                                >  2:853655/1‑140 (MQ=255)                                                       tAGACCGGGTCGTAACCTTTCTCGCTCAGCTGTTTCAGCCCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAg                                                                               <  2:405283/140‑1 (MQ=255)                                                          aCCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAg                                                                               >  2:789530/1‑137 (MQ=255)                                                          aCCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAg                                                                               <  1:789530/137‑1 (MQ=255)                                                          aCCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCgg                                                                            >  1:492954/1‑140 (MQ=255)                                                            cGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATGTGCGCAATCGAGGCTATGTGCAGTTCACCCAGCGGAt                                                                          >  1:250439/1‑140 (MQ=255)                                                            cGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGAt                                                                          <  2:594645/140‑1 (MQ=255)                                                               aTCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGa                                                                       >  2:769444/1‑140 (MQ=255)                                                                   tAACCGTTCTCGCTCAGCAGTTCCAGCGCCTCGTCTGAACTGTGGATTCCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGac                                                                   <  2:221085/140‑1 (MQ=255)                                                                    aaCCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGacc                                                                  <  2:208910/140‑1 (MQ=255)                                                                      ccGTTCTCGCTCAGCTGTTTCACCGCCTCCTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGaccac                                                                <  1:258701/140‑1 (MQ=255)                                                                         ttCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACcaccac                                                             >  1:140620/1‑140 (MQ=255)                                                                         ttCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGATGTCCAGCCAGCGTATGGAAGACcaccac                                                             >  1:248066/1‑140 (MQ=255)                                                                          tctcGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACt                                                            >  1:470523/1‑140 (MQ=255)                                                                            tcGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTc                                                          >  2:437760/1‑140 (MQ=255)                                                                              gCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCAt                                                        >  2:837226/1‑140 (MQ=255)                                                                                    cagTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGACAATGTGCAGTTAACCAAGCGGATGAAAGACCACCACTTCATCGATAc                                                  >  2:344916/1‑140 (MQ=255)                                                                                           aGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAAt                                           <  2:828940/140‑1 (MQ=255)                                                                                              gcCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGaa                                        >  2:434940/1‑140 (MQ=255)                                                                                                       gAAATGTGGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCcggacg                               <  2:492954/140‑1 (MQ=255)                                                                                                            gtgGATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCCGACGGAAGt                          <  2:140620/140‑1 (MQ=255)                                                                                                              ggATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTa                        <  1:837226/140‑1 (MQ=255)                                                                                                               gATTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTAt                       >  2:786198/1‑140 (MQ=255)                                                                                                                aTTTCATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATg                      <  1:769444/140‑1 (MQ=255)                                                                                                                    cATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATg                      <  2:54322/136‑1 (MQ=255)                                                                                                                    cATAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATg                      >  1:54322/1‑136 (MQ=255)                                                                                                                      tAACCACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTAc                >  2:802410/1‑140 (MQ=255)                                                                                                                         ccACGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCAATTCATCGACACGTTTAGTGAATTCCGGACGAAAGTTATGGCTTACcac             >  2:389970/1‑140 (MQ=255)                                                                                                                            cGTTCTTCCAGACGCTTGTACAGGCGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACAcc          >  2:142978/1‑140 (MQ=255)                                                                                                                            cGTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACAcc          >  2:656425/1‑140 (MQ=255)                                                                                                                             gTTCTTCCAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCg         >  2:669424/1‑140 (MQ=255)                                                                                                                             gTTCTTCCAGACGCTTGCACAGACGTTTCAACTGACTCTTCGCAATCGCGGCAATGTGCTGTTAACCCAGCGGATGGAAGACCACCACTTCATCAATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCg         <  2:248066/140‑1 (MQ=255)                                                                                                                                    cAGACGCTTGTACAGGCGTTTCGACTGAAGCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCa  <  1:389970/140‑1 (MQ=255)                                                                                                                                    cAGACGCTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCa  <  2:715983/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      GCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGACGTGCACCATAGACCGGATCGTAACCGTTCTCGCTCAGCAGTTTCAGCGCCTCGTCAGAAATGTGGATTTCATAACCACGTTCTTCCAGACGTTTGTACAGACGTTTCAACTGAATCTGCGCAATCGAGGCAATGTGCTGTTCACCCAGCGGATGGAAGACCACCACTTCATCGATACGGTTAATGAATTCCGGACGGAAGTTATGGCTTACCACACCGAGCACCA  >  NC_000913/2731685‑2731955

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