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

CTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGC  >  NC_000913/1620077‑1620341                                                                                                                                   |                                                                                                                                       cttcttGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCt                                                                                                                                >  1:1466986/1‑139 (MQ=255)        tgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCg                                                                                                                         >  2:1087474/1‑139 (MQ=255)                    ggTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGc                                                                                                             <  1:598107/139‑1 (MQ=255)                          tCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                        >  1:501642/1‑138 (MQ=255)                          tCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                        <  2:501642/138‑1 (MQ=255)                                aTCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                                                 >  2:1573944/1‑139 (MQ=255)                                    cAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTccc                                                                                                <  2:1660911/136‑1 (MQ=255)                                    cAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTccc                                                                                                >  1:1660911/1‑136 (MQ=255)                                            cccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                        <  2:3899/136‑1 (MQ=255)                                            cccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                        >  1:3899/1‑136 (MQ=255)                                                   tttGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCg                                                                                      <  2:1579833/131‑1 (MQ=255)                                                   tttGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCg                                                                                      >  1:1579833/1‑131 (MQ=255)                                                   tttGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGc                                                                              <  2:1466986/139‑1 (MQ=255)                                                          tAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCa                                                                       <  1:1573944/139‑1 (MQ=255)                                                           aaaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                      >  2:573076/1‑139 (MQ=255)                                                                  aTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                        >  1:255042/1‑114 (MQ=255)                                                                  aTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                        <  2:255042/114‑1 (MQ=255)                                                                          tGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGat                                                       >  2:1207023/1‑139 (MQ=255)                                                                                tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                                                 <  1:37652/91‑1 (MQ=255)                                                                                tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                                                 >  2:37652/1‑91 (MQ=255)                                                                                         ggATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                                <  2:558715/67‑1 (MQ=255)                                                                                         ggATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                                >  1:558715/1‑67 (MQ=255)                                                                                          gATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGG‑CCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCtt                                                                                   <  2:1318981/94‑1 (MQ=255)                                                                                          gATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGG‑CCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCtt                                                                                   >  1:1318981/1‑94 (MQ=255)                                                                                                cTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACa                                                                                                                      >  1:390137/1‑54 (MQ=255)                                                                                                cTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACa                                                                                                                      <  2:390137/54‑1 (MQ=255)                                                                                                cTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                                          >  2:1559776/1‑130 (MQ=255)                                                                                                cTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                                          <  1:1559776/130‑1 (MQ=255)                                                                                                  cGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcac                               <  1:1207023/139‑1 (MQ=255)                                                                                                   ggCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                        <  1:782550/97‑1 (MQ=255)                                                                                                   ggCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                        >  2:782550/1‑97 (MQ=255)                                                                                                     cGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTa                                                                                         <  1:732408/78‑1 (MQ=255)                                                                                                     cGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTa                                                                                         >  2:732408/1‑78 (MQ=255)                                                                                                                    tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAg             >  2:413257/1‑139 (MQ=255)                                                                                                                      aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                        >  2:1169124/1‑78 (MQ=255)                                                                                                                      aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                        <  1:1169124/78‑1 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                        >  2:522341/1‑44 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                        >  1:235711/1‑44 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                        <  1:522341/44‑1 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                        <  2:235711/44‑1 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                        >  2:84863/1‑76 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                        <  1:84863/76‑1 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                      <  1:684432/78‑1 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                                      >  2:684432/1‑78 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCa              >  2:1249175/1‑134 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCa              <  1:1249175/134‑1 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc         >  2:170530/1‑139 (MQ=255)                                                                                                                        cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc         >  2:1670001/1‑139 (MQ=255)                                                                                                                              aTGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                              >  1:1382074/1‑80 (MQ=255)                                                                                                                              aTGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                              <  2:1382074/80‑1 (MQ=255)                                                                                                                               tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCg                                                                                      >  2:1449753/1‑55 (MQ=255)                                                                                                                               tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCg                                                                                      <  1:1449753/55‑1 (MQ=255)                                                                                                                               tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCtt                                                                                   <  1:158058/58‑1 (MQ=255)                                                                                                                               tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCtt                                                                                   >  2:158058/1‑58 (MQ=255)                                                                                                                               tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGc  >  1:876196/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                       CTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGC  >  NC_000913/1620077‑1620341

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