Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,620,207	G→A	*73* (TAG→TAA)	marB →	periplasmic mar operon regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,620,207	0	G	A	100.0%	232.6 / NA	74	*73* (TAG→TAA)	marB	periplasmic mar operon regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/34);  total (40/34)

TCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGC  >  NC_000913/1620079‑1620341                                                                                                                                 |                                                                                                                                       tcttGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTgg                                                                                                                              <  1:59575/139‑1 (MQ=255)  cttGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGc                                                                                                                             >  2:1269923/1‑139 (MQ=255)    tgtgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGccc                                                                                                                           <  2:384551/139‑1 (MQ=255)          aaTGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACAt                                                                                                                     >  1:872462/1‑139 (MQ=255)            tGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACAtct                                                                                                                   >  1:118981/1‑139 (MQ=255)            tGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACAtct                                                                                                                   >  1:465528/1‑139 (MQ=255)               cGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                                <  2:245781/139‑1 (MQ=255)                       tttccccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                        >  1:670952/4‑139 (MQ=255)                       ttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                        >  1:853213/1‑139 (MQ=255)                       ttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGAGCttatt                                                                                                        <  1:1269923/139‑1 (MQ=255)                                cGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCa                                                                                               >  1:60416/1‑139 (MQ=255)                                cGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCa                                                                                               >  1:529531/1‑139 (MQ=255)                                 gCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCAc                                                                                              <  2:670952/139‑1 (MQ=255)                                                gTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACg                                                                               >  2:263385/1‑139 (MQ=255)                                                 tttGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGc                                                                              >  2:1381398/1‑139 (MQ=255)                                                       ttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                        <  1:1381398/139‑1 (MQ=255)                                                        tAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCgga                                                                       <  2:872462/139‑3 (MQ=255)                                                               caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaaa                                                                >  1:899546/1‑139 (MQ=255)                                                                   ggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                            <  2:465528/139‑1 (MQ=255)                                                                            aGTGATAAGTCGGATGTGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCggg                                                                                                               <  2:982756/79‑1 (MQ=255)                                                                            aGTGATAAGTCGGATGTGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCggg                                                                                                               >  1:982756/1‑79 (MQ=255)                                                                              tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGa                                                                                                                        <  2:1488605/68‑1 (MQ=255)                                                                              tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGa                                                                                                                        >  1:1488605/1‑68 (MQ=255)                                                                              tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACtt                                                                                                   <  1:810415/89‑1 (MQ=255)                                                                              tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACtt                                                                                                   >  2:810415/1‑89 (MQ=255)                                                                               gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                        <  1:1133608/99‑1 (MQ=255)                                                                               gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAAATTCCCACCTTTAc                                                                                        >  2:1133608/1‑99 (MQ=255)                                                                                            cgcTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCggg                                                                                                               <  1:1130510/63‑1 (MQ=255)                                                                                            cgcTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCggg                                                                                                               >  2:1130510/1‑63 (MQ=255)                                                                                               tCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                                >  2:310847/1‑59 (MQ=255)                                                                                               tCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                                <  1:310847/59‑1 (MQ=255)                                                                                                   cGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGAc                                                                                                                       >  2:96015/1‑48 (MQ=255)                                                                                                   cGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGAc                                                                                                                       <  1:96015/48‑1 (MQ=255)                                                                                                           attatAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                                          <  1:1347664/117‑1 (MQ=255)                                                                                                           attatAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                                          >  2:1347664/1‑117 (MQ=255)                                                                                                                 aTCAACACGCTATGTAATTTGTTCTGGCCCCGCCATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCgg                                            >  2:1198854/1‑109 (MQ=255)                                                                                                                 aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCgg                                            <  1:1198854/109‑1 (MQ=255)                                                                                                                 aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCa              <  2:118981/139‑1 (MQ=255)                                                                                                                  tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                              >  2:1258073/1‑90 (MQ=255)                                                                                                                  tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                              <  1:1258073/90‑1 (MQ=255)                                                                                                                    aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                                   <  1:442956/83‑1 (MQ=255)                                                                                                                    aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                                   >  2:442956/1‑83 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                                >  1:166026/1‑36 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                                <  2:166026/36‑1 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGc                                                                                                             <  1:300434/39‑1 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGc                                                                                                             >  2:300434/1‑39 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                        <  1:946425/60‑1 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                        >  2:946425/1‑60 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTAc                                                                                >  1:322330/1‑68 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTAc                                                                                <  2:322330/68‑1 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGat                                                       <  1:40072/93‑1 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGat                                                       >  2:40072/1‑93 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                               >  1:192901/1‑101 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                               <  2:192901/101‑1 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc         >  1:86618/1‑139 (MQ=255)                                                                                                                      cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc         >  2:1016066/1‑139 (MQ=255)                                                                                                                      cacGCTATGTAATTTGCTCTGGCCCCGACATCTCGGGGCt                                                                                                            >  1:670047/1‑40 (MQ=255)                                                                                                                      cacGCTAGGTAATTTGTTCTGGCCCCGACATCTCGGGGCt                                                                                                            <  2:670047/40‑1 (MQ=255)                                                                                                                         gCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                        >  2:1400458/1‑57 (MQ=255)                                                                                                                         gCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                        <  1:1400458/57‑1 (MQ=255)                                                                                                                         gCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                            >  2:26293/1‑85 (MQ=255)                                                                                                                         gCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                            <  1:26293/85‑1 (MQ=255)                                                                                                                         gCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                            <  1:1257866/85‑1 (MQ=255)                                                                                                                         gCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                            >  2:1257866/1‑85 (MQ=255)                                                                                                                          cTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                               >  2:596598/1‑97 (MQ=255)                                                                                                                          cTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                               <  1:596598/97‑1 (MQ=255)                                                                                                                            aTGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                              <  1:1378018/80‑1 (MQ=255)                                                                                                                            aTGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                              >  2:1378018/1‑80 (MQ=255)                                                                                                                             tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAAt                                     >  2:642469/1‑104 (MQ=255)                                                                                                                             tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAAt                                     <  1:642469/104‑1 (MQ=255)                                                                                                                             tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGc  >  1:949011/1‑139 (MQ=255)                                                                                                                             tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGc  >  1:280116/1‑139 (MQ=255)                                                                                                                             tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGc  <  2:468204/139‑1 (MQ=255)                                                                                                                              gTAATTTCTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                              >  2:997109/1‑78 (MQ=255)                                                                                                                                 |                                                                                                                                       TCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGC  >  NC_000913/1620079‑1620341

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