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

TCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAAC  >  NC_000913/2637163‑2637421                                                                                                                                     |                                                                                                                               tcttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaa                                                                                                                          >  1:948246/1‑139 (MQ=255)    tctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATa                                                                                                                       >  2:1538807/1‑139 (MQ=255)          ccATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTgg                                                                                                                 >  2:129929/1‑139 (MQ=255)                aTCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                           <  2:948246/139‑1 (MQ=255)                      ccgGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGt                                                                                                     >  2:441056/1‑139 (MQ=255)                         gTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAg                                                                                                  >  1:836801/1‑139 (MQ=255)                           aTCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                                                <  1:474158/139‑1 (MQ=255)                             cAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcg                                                                                              <  1:1852843/139‑1 (MQ=255)                                     aTAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACaa                                                                                      >  2:1752300/1‑139 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAgcg                                                                                   <  2:1836630/137‑1 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAgcg                                                                                   >  1:1836630/1‑137 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                 <  1:441056/139‑1 (MQ=255)                                          cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                 <  1:129929/139‑1 (MQ=255)                                            cACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                                               <  2:836801/139‑1 (MQ=255)                                             aCGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAcc                                                                              >  2:1531276/1‑139 (MQ=255)                                              cGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACca                                                                             >  2:348340/1‑139 (MQ=255)                                                     cAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcgc                                                                                             >  2:1770760/1‑116 (MQ=255)                                                     cAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGgcgc                                                                                             <  1:1770760/116‑1 (MQ=255)                                                           ccGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATgg                                                                                                                              <  2:320675/77‑1 (MQ=255)                                                           ccGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATgg                                                                                                                              >  1:320675/1‑77 (MQ=255)                                                           ccGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCa                                                                    <  2:1190419/135‑1 (MQ=255)                                                           ccGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCa                                                                    >  1:1190419/1‑135 (MQ=255)                                                           ccGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGtt                                                                >  1:89333/1‑139 (MQ=255)                                                                      tACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAt                                                     >  2:949885/1‑139 (MQ=255)                                                                            cGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAAc                                                                                                                     >  2:658023/1‑69 (MQ=255)                                                                            cGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAAc                                                                                                                     <  1:658023/69‑1 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                           >  2:1289022/1‑78 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                           <  2:865539/78‑1 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                           >  1:865539/1‑78 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                           <  1:1289022/78‑1 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                       >  1:708509/1‑98 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                       <  2:708509/98‑1 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                 <  1:331389/104‑1 (MQ=255)                                                                             ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                 >  2:331389/1‑104 (MQ=255)                                                                                 aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                           >  2:311579/1‑74 (MQ=255)                                                                                 aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                           <  2:280766/74‑1 (MQ=255)                                                                                 aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                           <  1:311579/74‑1 (MQ=255)                                                                                 aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                           >  1:280766/1‑74 (MQ=255)                                                                                 aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAatc                                          >  1:1711510/1‑139 (MQ=255)                                                                                  cGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAatca                                         >  1:620534/1‑139 (MQ=255)                                                                                     gTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAAtcatca                                      >  1:493781/1‑139 (MQ=255)                                                                                         gACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                              <  2:1221053/127‑1 (MQ=255)                                                                                         gACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                              >  1:1221053/1‑127 (MQ=255)                                                                                            ccGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacg                               <  1:949885/139‑1 (MQ=255)                                                                                            ccGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacg                               <  2:620534/139‑1 (MQ=255)                                                                                            ccGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacg                               <  2:89333/139‑1 (MQ=255)                                                                                            ccGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacg                               <  1:348340/139‑1 (MQ=255)                                                                                              gggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca                             <  1:1538807/139‑1 (MQ=255)                                                                                               ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACca                                                                             >  1:886413/1‑90 (MQ=255)                                                                                               ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACca                                                                             <  2:886413/90‑1 (MQ=255)                                                                                               ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGt                                                                        <  2:1456588/95‑1 (MQ=255)                                                                                               ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGt                                                                        >  1:1456588/1‑95 (MQ=255)                                                                                               ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGa                                    >  2:42746/1‑131 (MQ=255)                                                                                               ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGa                                    <  1:42746/131‑1 (MQ=255)                                                                                                gAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgcacg                           <  2:1711510/139‑1 (MQ=255)                                                                                                gAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgcacg                           <  1:1531276/139‑1 (MQ=255)                                                                                                    tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTa                                                                                                   >  1:862416/1‑63 (MQ=255)                                                                                                    tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTa                                                                                                   <  2:862416/63‑1 (MQ=255)                                                                                                               cgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAgcg                                                                                   >  1:827701/1‑68 (MQ=255)                                                                                                               cgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAgcg                                                                                   <  2:827701/68‑1 (MQ=255)                                                                                                               cgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGcgtcgt            >  1:1415635/1‑139 (MQ=255)                                                                                                                         gTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAAc  >  1:317380/1‑139 (MQ=255)                                                                                                                              aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                                                >  1:991945/1‑40 (MQ=255)                                                                                                                              aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                                                <  2:991945/40‑1 (MQ=255)                                                                                                                              aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgcacgc                          >  1:17947/1‑110 (MQ=255)                                                                                                                              aGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgcacgc                          <  2:17947/110‑1 (MQ=255)                                                                                                                                     |                                                                                                                               TCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAAC  >  NC_000913/2637163‑2637421

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