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%	231.4 / NA	71	*73* (TAG→TAA)	marB	periplasmic mar operon regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/33);  total (38/33)

TGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAA  >  NC_000913/1620082‑1620339                                                                                                                              |                                                                                                                                     tgtgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGccc                                                                                                                         >  1:553886/1‑139 (MQ=255) tgtgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGccc                                                                                                                         >  2:543828/1‑139 (MQ=255)  gtgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGcccc                                                                                                                        <  2:1443743/139‑1 (MQ=255)   tgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCg                                                                                                                       >  1:1486518/1‑139 (MQ=255)   tgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCg                                                                                                                       >  2:1606399/1‑139 (MQ=255)         tGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACAtct                                                                                                                 <  1:1019137/139‑1 (MQ=255)         tGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACAtct                                                                                                                 <  2:1187451/139‑1 (MQ=255)            cGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                              <  2:1486518/139‑1 (MQ=255)             gTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCggg                                                                                                             >  1:536657/1‑139 (MQ=255)              tGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                              >  2:1146841/1‑137 (MQ=255)              tGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                              <  1:1146841/137‑1 (MQ=255)                   gttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGACGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttat                                                                                                       >  1:1151157/1‑139 (MQ=255)                     tCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGtt                                                                                                                                <  2:616206/112‑1 (MQ=255)                     tCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGtt                                                                                                                                >  1:616206/1‑112 (MQ=255)                     tCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatta                                                                                                     >  2:654042/1‑139 (MQ=255)                      cccccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTaa                                                                                                    <  1:654042/139‑1 (MQ=255)                             cGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGA‑‑‑CTCGGGGCTTATTAACTTCCCACCt                                                                                          >  2:601774/1‑139 (MQ=255)                              gCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCAc                                                                                            >  1:197410/1‑139 (MQ=255)                                 ggAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCtt                                                                                         <  1:1606399/139‑1 (MQ=255)                                  gAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCttt                                                                                        <  2:1151157/139‑1 (MQ=255)                                    acacCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                      <  2:536657/139‑1 (MQ=255)                                       cccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCTGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                              <  1:182019/112‑1 (MQ=255)                                       cccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                              >  2:182019/1‑112 (MQ=255)                                       cccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTa                                                                                       <  1:503654/135‑1 (MQ=255)                                       cccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTa                                                                                       >  2:503654/1‑135 (MQ=255)                                        ccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCa                                                                                             <  1:1398493/128‑1 (MQ=255)                                        ccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCa                                                                                             >  2:1398493/1‑128 (MQ=255)                                            cGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTAc                                                                              >  2:1450864/1‑139 (MQ=255)                                                   tttAAATCACATGGGTACTGGCATTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTT‑CCACCTTTACCGCTTTACGCCACCGc                                                                      <  2:103244/139‑1 (MQ=255)                                                         tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                                 <  1:1450864/139‑1 (MQ=255)                                                         tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                                 <  2:543211/139‑1 (MQ=255)                                                            caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaaa                                                              <  1:327026/139‑1 (MQ=255)                                                             attggTACTGGGAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCTACACGCTATGTAATTTGtt                                                                                                                                >  1:933226/4‑72 (MQ=255)                                                             aTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGtt                                                                                                                                <  2:933226/72‑1 (MQ=255)                                                                ggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                          <  1:573321/139‑1 (MQ=255)                                                                 gTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAt                                                         >  2:404776/1‑139 (MQ=255)                                                                           tGATAAGTCGGATGCACTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttat                                                                                                       <  2:320168/83‑1 (MQ=255)                                                                           tGATAAGTCGGATGCACTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttat                                                                                                       >  1:320168/1‑83 (MQ=255)                                                                            gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                                 >  2:1331949/1‑120 (MQ=255)                                                                            gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                                 <  1:1331949/120‑1 (MQ=255)                                                                            gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGc                                              >  1:822127/1‑139 (MQ=255)                                                                            gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGc                                              >  1:716802/1‑139 (MQ=255)                                                                             aTAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                             <  2:677989/139‑1 (MQ=255)                                                                                 gTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGt                                         >  2:864042/1‑139 (MQ=255)                                                                                  tCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                                        <  2:822127/139‑1 (MQ=255)                                                                                  tCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                                        <  1:864042/139‑1 (MQ=255)                                                                                     gATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATa                                     >  2:1360542/1‑139 (MQ=255)                                                                                             cGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcac                             <  2:553886/139‑1 (MQ=255)                                                                                                cGTGCCCTATTATAATCAACACGCTATGTAATTTGtt                                                                                                                                >  1:1613862/1‑37 (MQ=255)                                                                                                cGTGCCCTATTATAATCAACACGCTATGTAATTTGtt                                                                                                                                <  2:1613862/37‑1 (MQ=255)                                                                                                cGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc                          >  1:672987/1‑139 (MQ=255)                                                                                                    cccTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCt                      >  2:739399/1‑139 (MQ=255)                                                                                                            taATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTg              >  2:565143/1‑139 (MQ=255)                                                                                                               tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                             <  1:1319564/105‑1 (MQ=255)                                                                                                               tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                             >  2:1319564/1‑105 (MQ=255)                                                                                                                 aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                          <  1:1054848/90‑1 (MQ=255)                                                                                                                 aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                          >  2:1054848/1‑90 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                              <  2:658460/36‑1 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                              >  1:658460/1‑36 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                      <  1:1420527/76‑1 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                      >  2:1420527/1‑76 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                      >  2:898889/1‑76 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                      <  1:898889/76‑1 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                            >  1:136898/1‑86 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                            <  2:136898/86‑1 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc       >  2:507607/1‑139 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc       >  1:126101/1‑139 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc       >  2:779400/1‑139 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCT‑AAAATTGCa            <  2:545655/133‑1 (MQ=255)                                                                                                                   cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCT‑AAAATTGCa            >  1:545655/1‑133 (MQ=255)                                                                                                                        tATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTaa  >  1:787766/1‑139 (MQ=255)                                                                                                                              |                                                                                                                                     TGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAA  >  NC_000913/1620082‑1620339

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