Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,637,295	C→T	R21H (CGT→CAT)	der ←	GTPase, multicopy suppressor of ftsJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,637,295	0	C	T	100.0%	164.5 / NA	51	R21H (CGT→CAT)	der	GTPase, multicopy suppressor of ftsJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/26);  total (25/26)

CGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCA  >  NC_000913/2637161‑2637394                                                                                                                                       |                                                                                                    cGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGtt                                                                                                 >  1:257601/1‑139 (MQ=255)   tcttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaa                                                                                               >  1:357954/1‑139 (MQ=255)    cttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaaa                                                                                              <  2:235195/139‑1 (MQ=255)                          ggTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTa                                                                        <  1:837072/139‑1 (MQ=255)                             aTCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                     <  2:257601/139‑1 (MQ=255)                               cAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcg                                                                   <  1:549142/139‑1 (MQ=255)                                       aTAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACaa                                                           >  2:921831/1‑139 (MQ=255)                                       aTAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACaa                                                           >  2:364760/1‑139 (MQ=255)                                        tAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAg                                                          >  1:323801/1‑139 (MQ=255)                                          aaCTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                            >  2:658544/1‑135 (MQ=255)                                          aaCTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                            <  1:658544/135‑1 (MQ=255)                                          aaCTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAgcg                                                        <  2:117555/139‑1 (MQ=255)                                            cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                     <  2:391513/108‑1 (MQ=255)                                            cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                     >  1:391513/1‑108 (MQ=255)                                              cACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                    <  1:688238/139‑1 (MQ=255)                                                 ggCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACcac                                                 >  2:480839/1‑139 (MQ=255)                                                 ggCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACcac                                                 >  1:492443/1‑139 (MQ=255)                                                  gCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACcaca                                                <  2:492443/139‑1 (MQ=255)                                                       cAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAc                                           <  1:921831/139‑1 (MQ=255)                                                                             aCGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTta                     >  2:273611/1‑139 (MQ=255)                                                                               ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                    <  1:310772/106‑1 (MQ=255)                                                                               ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                    >  2:310772/1‑106 (MQ=255)                                                                               ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACca                                                  >  1:315566/1‑108 (MQ=255)                                                                               ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACca                                                  <  2:315566/108‑1 (MQ=255)                                                                               ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                   <  1:273611/139‑1 (MQ=255)                                                                               ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                   >  2:88905/1‑139 (MQ=255)                                                                               ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                   >  1:338852/1‑139 (MQ=255)                                                                                  cACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                     >  1:771319/1‑86 (MQ=255)                                                                                  cACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                     <  2:771319/86‑1 (MQ=255)                                                                                   aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                     >  2:487559/1‑69 (MQ=255)                                                                                   aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                     <  1:487559/69‑1 (MQ=255)                                                                                   aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTCCGTTAGGGCGCCCGACAAGCGCGAc                                                    >  2:664824/1‑102 (MQ=255)                                                                                   aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                            <  2:574548/94‑1 (MQ=255)                                                                                   aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                            >  1:574548/1‑94 (MQ=255)                                                                                   aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                    <  1:664824/102‑1 (MQ=255)                                                                                   aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAatc               >  1:856559/1‑139 (MQ=255)                                                                                      aGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAatcatc            <  2:338852/139‑1 (MQ=255)                                                                                                gggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca  <  2:856559/139‑1 (MQ=255)                                                                                                gggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca  <  1:480839/139‑1 (MQ=255)                                                                                                 ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                    <  2:136637/56‑1 (MQ=255)                                                                                                 ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                    >  1:136637/1‑56 (MQ=255)                                                                                                 ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAcc                                          >  2:683061/1‑98 (MQ=255)                                                                                                 ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAcc                                          <  1:683061/98‑1 (MQ=255)                                                                                                 ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGAc        <  1:21152/132‑1 (MQ=255)                                                                                                 ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGAc        >  2:21152/1‑132 (MQ=255)                                                                                                                 cgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAc                                           <  2:345196/81‑1 (MQ=255)                                                                                                                 cgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAc                                           >  1:345196/1‑81 (MQ=255)                                                                                                                       tcGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca  >  1:332572/1‑116 (MQ=255)                                                                                                                       tcGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca  <  2:332572/116‑1 (MQ=255)                                                                                                                                aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAc                                           <  1:247626/66‑1 (MQ=255)                                                                                                                                aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAc                                           >  2:247626/1‑66 (MQ=255)                                                                                                                                       |                                                                                                    CGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCA  >  NC_000913/2637161‑2637394

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