Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,530,484	T→C	K250K (AAA→AAG)	zipA ←	cell division protein ZipA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,530,484	0	T	C	100.0%	187.5 / NA	54	K250K (AAA→AAG)	zipA	cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/27);  total (27/27)

TTCATCGGCAATATGCTGCGCAGATTGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGTTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTC  >  NC_000913/2530349‑2530618                                                                                                                                        |                                                                                                                                       ttCATCGGCAATATGCTGCGCAGATTGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCa                                                                                                                                    <  2:453298/140‑1 (MQ=255) ttCATCGGCAATATGCTGCGCAGATTGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCa                                                                                                                                    <  2:323594/140‑1 (MQ=255)    aTCGGCAATATGCTGCGCAGATTGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCa                                                                                                                                    >  1:915752/1‑137 (MQ=255)    aTCGGCAATATGCTGCGCAGATTGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCa                                                                                                                                    <  2:915752/137‑1 (MQ=255)      cGGCAATATGCTGCGCAGATTGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCata                                                                                                                               >  2:986298/1‑140 (MQ=255)             atGCTGCGCAGATTGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCa                                                                                                                        >  1:602890/1‑140 (MQ=255)               gctgcGCAGATTGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAgg                                                                                                                      <  2:164382/140‑1 (MQ=255)               gctgcGCAGATTGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAgg                                                                                                                      <  1:560473/140‑1 (MQ=255)                        aTTGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAAc                                                                                                             >  2:262347/1‑140 (MQ=255)                          tGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGc                                                                                                           >  1:467326/1‑140 (MQ=255)                                 tCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCg                                                                                                    >  1:454389/1‑140 (MQ=255)                                    gCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTg                                                                                                 <  1:484443/140‑1 (MQ=255)                                        gAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCAt                                                                                             <  2:454389/140‑1 (MQ=255)                                         aaGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATc                                                                                            >  1:313519/1‑140 (MQ=255)                                         aaGTTCTGCAGCTCGTCACCGTAAGACGGAACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATc                                                                                            <  2:303503/140‑1 (MQ=255)                                             tCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCggg                                                                                        <  1:775090/140‑1 (MQ=255)                                                 cAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTaa                                                                                    <  2:714279/140‑1 (MQ=255)                                                        tCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACg                                                                             <  1:585593/140‑1 (MQ=255)                                                         cACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGa                                                                            <  2:501318/140‑1 (MQ=255)                                                             gTAAGACGGTACCTGCATAAAGTTAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGt                                                                        <  1:705967/140‑1 (MQ=255)                                                             gTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGt                                                                        <  2:313519/140‑1 (MQ=255)                                                              tAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGa                                                                            <  1:400703/135‑1 (MQ=255)                                                              tAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGATCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGa                                                                            >  2:400703/1‑135 (MQ=255)                                                                      tACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCa                                                               >  2:522143/1‑140 (MQ=255)                                                                        ccTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCata                                                             >  1:392006/1‑140 (MQ=255)                                                                           gCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCg                                                          >  1:347056/1‑140 (MQ=255)                                                                              tAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAGGATGACGATGGt                                                                        >  1:965753/1‑123 (MQ=255)                                                                              tAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAGGATGACGATGGt                                                                        <  2:965753/123‑1 (MQ=255)                                                                                aaGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGACTACGATGGTAAATATTCATATCGCCaaa                                                     >  2:589547/1‑140 (MQ=255)                                                                                      gTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGaa                                               >  1:976739/1‑140 (MQ=255)                                                                                          cacCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGccc                                           >  2:361394/1‑140 (MQ=255)                                                                                           acCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCg                                          >  1:414962/1‑140 (MQ=255)                                                                                            cccGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGc                                         <  2:467326/140‑1 (MQ=255)                                                                                            cccGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGc                                         <  1:262347/140‑1 (MQ=255)                                                                                                 agtagtGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAg                                              >  2:484620/1‑130 (MQ=255)                                                                                                 agtagtGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAg                                              <  1:484620/130‑1 (MQ=255)                                                                                                 agtagtGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCttgtt                                    >  2:992020/1‑140 (MQ=255)                                                                                                   tagtGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGa                                  <  1:361394/140‑1 (MQ=255)                                                                                                    agtGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGaa                                 <  2:347056/140‑1 (MQ=255)                                                                                                      tGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATg                               >  2:707358/1‑140 (MQ=255)                                                                                                       gAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGc                              >  2:426520/1‑140 (MQ=255)                                                                                                       gAAATCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGc                              <  1:426520/140‑1 (MQ=255)                                                                                                           tCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGc                              <  1:618500/136‑1 (MQ=255)                                                                                                           tCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGc                              >  2:618500/1‑136 (MQ=255)                                                                                                           tCCTTCATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGtt                          >  2:19952/1‑140 (MQ=255)                                                                                                                cATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATg                               <  2:304873/130‑1 (MQ=255)                                                                                                                cATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATg                               >  1:304873/1‑130 (MQ=255)                                                                                                                cATTTCAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTaagaa                     >  1:252577/1‑140 (MQ=255)                                                                                                                     cAGGATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCAGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGt                >  2:952564/1‑140 (MQ=255)                                                                                                                       ggATCAAAGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTc              >  1:313304/1‑140 (MQ=255)                                                                                                                             aaGGTTCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGtt        >  1:593837/1‑140 (MQ=255)                                                                                                                                 ttCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGc    <  1:992020/140‑1 (MQ=255)                                                                                                                                  tCCCGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCt   <  2:313304/140‑1 (MQ=255)                                                                                                                                   cccGGCTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGTAGCCCGCTTGTTGGATGCTGTTAAGAAGCAGTTCACCGTTTAGCTc  <  2:602890/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       TTCATCGGCAATATGCTGCGCAGATTGCAGCATCAGCTTGAAGTTCTGCAGCTCGTCACCGTAAGACGGTACCTGCATAAAGATAGTGACACCCGGAGTAGTGAAATCCTTCATTTCAGGATCAAAGGTTCCCGGTTTCACCATATTCGCCAGGCTGAATAACGCCGGGCCGCTGCCATCCGGGCTAAGATGACGATGGTAAATATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTC  >  NC_000913/2530349‑2530618

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: ‑