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

TCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAAC  >  NC_000913/2637163‑2637421                                                                                                                                      |                                                                                                                               tcttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTaa                                                                                                                          >  1:554414/1‑139 (MQ=255) tcttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTaa                                                                                                                          >  1:258029/1‑139 (MQ=255) tcttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTaa                                                                                                                          >  1:1327406/1‑139 (MQ=255)      tgtgCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAAc                                                                                                                     >  1:1218109/1‑139 (MQ=255)             tCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGAtt                                                                                                              >  2:875347/1‑139 (MQ=255)                aTCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                           <  2:1340294/139‑1 (MQ=255)                            tCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgc                                                                                               >  1:1508533/1‑139 (MQ=255)                                    gATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                       >  1:1220180/1‑139 (MQ=255)                                    gATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                       >  1:1410698/1‑139 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                 >  2:906490/1‑139 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                 <  2:554414/139‑1 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                 >  1:334155/1‑139 (MQ=255)                                             aCGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAcc                                                                              >  1:1416210/1‑139 (MQ=255)                                             aCGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAcc                                                                              >  1:118290/1‑139 (MQ=255)                                              cGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACca                                                                             <  1:507707/139‑1 (MQ=255)                                               ggCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACcac                                                                            >  1:1179796/1‑139 (MQ=255)                                               ggCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACcac                                                                            <  1:335302/139‑1 (MQ=255)                                               ggCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACcac                                                                            >  2:335302/1‑139 (MQ=255)                                                 ccTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAg                                                                          <  2:1179796/139‑1 (MQ=255)                                                      aaTTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAcc                                                                     <  1:906490/139‑1 (MQ=255)                                                             gCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTa                                                              >  1:1114790/1‑139 (MQ=255)                                                              cACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAAc                                                                                                                     >  2:1519410/1‑83 (MQ=255)                                                              cACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAAc                                                                                                                     <  1:1519410/83‑1 (MQ=255)                                                                 gACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGc                                                          <  2:334155/139‑1 (MQ=255)                                                                  aCCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAAc                                                                                                                     <  2:530123/79‑1 (MQ=255)                                                                  aCCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAAc                                                                                                                     >  1:530123/1‑79 (MQ=255)                                                                           aCGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                                                >  1:1171054/1‑75 (MQ=255)                                                                           aCGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                                                <  2:1171054/75‑1 (MQ=255)                                                                           aCGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTta                                                <  2:1508533/139‑1 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                 >  1:792420/1‑104 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                 <  2:792420/104‑1 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                                               >  1:166063/1‑106 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                                               <  2:166063/106‑1 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                              >  1:1324636/1‑139 (MQ=255)                                                                                      tCAGACCCGGGAAACCCGCAACCCGCGCGTCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGCCAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAg                                     <  2:1217098/139‑1 (MQ=255)                                                                                            ccGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacg                               <  2:1114790/139‑1 (MQ=255)                                                                                              gggAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca                             <  2:258029/139‑1 (MQ=255)                                                                                               ggAAATCCGCAACCAGCGCATCTCGGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacg                               >  1:1271343/1‑137 (MQ=255)                                                                                               ggAAATCCGCAACCAGCGCATCTCGGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacg                               <  2:1271343/137‑1 (MQ=255)                                                                                               ggAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAg                                                                          >  2:1021482/1‑93 (MQ=255)                                                                                               ggAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAg                                                                          <  1:1021482/93‑1 (MQ=255)                                                                                                               cgcATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAg                                                                                                  <  1:1363085/53‑1 (MQ=255)                                                                                                               cgcATCTC‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAg                                                                                                  >  2:1363085/1‑53 (MQ=255)                                                                                                                    ctc‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcg                                                                                              <  1:650848/52‑1 (MQ=255)                                                                                                                    ctc‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcg                                                                                              >  2:650848/1‑52 (MQ=255)                                                                                                                    ctc‑GGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTa       <  2:1220180/139‑1 (MQ=255)                                                                                                                         ggTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAcc                                                                     <  2:1015878/73‑1 (MQ=255)                                                                                                                         ggTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAcc                                                                     >  1:1015878/1‑73 (MQ=255)                                                                                                                          gTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                                                <  2:645017/45‑1 (MQ=255)                                                                                                                          gTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                                                >  1:645017/1‑45 (MQ=255)                                                                                                                          gTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                       <  2:980707/54‑1 (MQ=255)                                                                                                                          gTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                       >  1:980707/1‑54 (MQ=255)                                                                                                                          gTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgcacg                           >  1:806516/1‑114 (MQ=255)                                                                                                                          gTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgcacg                           <  2:806516/114‑1 (MQ=255)                                                                                                                          gTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAAc  >  1:755914/1‑139 (MQ=255)                                                                                                                             cGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                              <  2:617850/92‑1 (MQ=255)                                                                                                                             cGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                              >  1:617850/1‑92 (MQ=255)                                                                                                                               aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca                             >  2:635245/1‑107 (MQ=255)                                                                                                                               aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca                             <  1:635245/107‑1 (MQ=255)                                                                                                                                      |                                                                                                                               TCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTC‑GGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAAC  >  NC_000913/2637163‑2637421

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