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%	195.2 / NA	61	R21H (CGT→CAT)	der	GTPase, multicopy suppressor of ftsJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/31);  total (30/31)

TCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAA  >  NC_000913/2637163‑2637420                                                                                                                                     |                                                                                                                              tcttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaa                                                                                                                         >  1:1023271/1‑139 (MQ=255) tcttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaa                                                                                                                         >  2:910491/1‑139 (MQ=255) tcttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaa                                                                                                                         >  1:327765/1‑139 (MQ=255)  cttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaaa                                                                                                                        >  2:157433/1‑139 (MQ=255)    tctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATa                                                                                                                      >  1:207143/1‑139 (MQ=255)            aTCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                                              <  2:1278042/139‑1 (MQ=255)            aTCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                                              <  2:1164191/139‑1 (MQ=255)                      ccgGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGt                                                                                                    >  2:710946/1‑139 (MQ=255)                             cAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcg                                                                                             <  2:327765/139‑1 (MQ=255)                             cAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcg                                                                                             <  2:1023271/139‑1 (MQ=255)                                tACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGccc                                                                                          >  1:572663/1‑139 (MQ=255)                                   aGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcg                                                                                             >  1:101480/1‑133 (MQ=255)                                   aGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcg                                                                                             <  2:101480/133‑1 (MQ=255)                                    gATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                      >  2:942315/1‑139 (MQ=255)                                    gATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                      >  1:304562/1‑139 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAg                                                                                                 <  2:10658/122‑1 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAg                                                                                                 >  1:10658/1‑122 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                <  2:304562/139‑1 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                <  1:710946/139‑1 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                >  1:646336/1‑139 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                >  2:149703/1‑139 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                >  2:135265/1‑139 (MQ=255)                                            cACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                                              <  2:207143/139‑1 (MQ=255)                                            cACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                                              <  1:942315/139‑1 (MQ=255)                                            cACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                                              <  1:157433/139‑1 (MQ=255)                                             aCGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                                              <  1:782222/106‑1 (MQ=255)                                             aCGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                                              >  2:782222/1‑106 (MQ=255)                                             aCGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                      <  2:931714/130‑1 (MQ=255)                                             aCGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                      >  1:931714/1‑130 (MQ=255)                                               ggCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACcac                                                                           <  1:149703/139‑1 (MQ=255)                                               ggCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCACATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcg                                                                                             >  2:897090/1‑121 (MQ=255)                                               ggCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCACATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcg                                                                                             <  1:897090/121‑1 (MQ=255)                                                                     gTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCa                                                     <  2:572663/139‑1 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                             >  2:766251/1‑139 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                             <  2:128361/139‑1 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGATGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                             >  2:915104/1‑139 (MQ=255)                                                                                 aagAGTCAGACCCGGGAAATCGGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAtt                                                                                                             >  2:777580/3‑71 (MQ=255)                                                                                 aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAtt                                                                                                             <  1:777580/71‑1 (MQ=255)                                                                                 aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                          >  1:974690/1‑74 (MQ=255)                                                                                 aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                          <  2:974690/74‑1 (MQ=255)                                                                                 aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAatc                                         >  1:502132/1‑139 (MQ=255)                                                                                  cGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAatca                                        >  2:171161/1‑139 (MQ=255)                                                                                       cAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGa                                   <  2:646336/139‑1 (MQ=255)                                                                                            ccGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacg                              <  1:171161/139‑1 (MQ=255)                                                                                            ccGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacg                              <  1:135265/139‑1 (MQ=255)                                                                                            ccGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacg                              <  1:915104/139‑1 (MQ=255)                                                                                              gggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca                            <  2:502132/139‑1 (MQ=255)                                                                                               ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgcac                           >  2:415962/1‑139 (MQ=255)                                                                                                    tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaaa                                                                                                                        <  1:676517/41‑1 (MQ=255)                                                                                                    tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaaa                                                                                                                        >  2:676517/1‑41 (MQ=255)                                                                                                    tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAt                                                    >  2:406564/1‑109 (MQ=255)                                                                                                    tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAt                                                    <  1:406564/109‑1 (MQ=255)                                                                                                               cgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTgg                                                                                                                >  1:302880/1‑38 (MQ=255)                                                                                                               cgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTgg                                                                                                                <  2:302880/38‑1 (MQ=255)                                                                                                               cgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGc                                                         <  1:951231/93‑1 (MQ=255)                                                                                                               cgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGc                                                         >  2:951231/1‑93 (MQ=255)                                                                                                                        ggTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTaaaaa  <  1:415962/139‑1 (MQ=255)                                                                                                                              aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                             <  2:711005/90‑1 (MQ=255)                                                                                                                              aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                             >  1:711005/1‑90 (MQ=255)                                                                                                                              aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgcacg                          <  2:1266224/109‑1 (MQ=255)                                                                                                                              aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgcacg                          >  1:1266224/1‑109 (MQ=255)                                                                                                                                     |                                                                                                                              TCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAA  >  NC_000913/2637163‑2637420

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