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

ACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCC  >  NC_000913/1620076‑1620334                                                                                                                                    |                                                                                                                                aCTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTc                                                                                                                          <  1:464132/139‑1 (MQ=255)      ttGTGTCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGcc                                                                                                                     >  1:896058/1‑139 (MQ=255)         tgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGtt                                                                                                                           >  2:895468/1‑130 (MQ=255)         tgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGtt                                                                                                                           <  1:895468/130‑1 (MQ=255)         tgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCg                                                                                                                  >  1:1825013/1‑139 (MQ=255)                gTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACAtctc                                                                                                           <  1:1310406/139‑1 (MQ=255)                  cGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                         <  2:372196/139‑1 (MQ=255)                           tCCCCCATCGCAGGAACACCTACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatta                                                                                                >  2:110454/1‑139 (MQ=255)                                 aTCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                                          >  1:293512/1‑139 (MQ=255)                                  tCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGAc                                                                                                                <  2:617431/116‑1 (MQ=255)                                  tCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGAc                                                                                                                >  1:617431/1‑116 (MQ=255)                                          acacCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                 <  1:794123/139‑1 (MQ=255)                                             cccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGc                                                                              >  2:1490517/1‑139 (MQ=255)                                              ccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCg                                                                                                                  <  1:157440/102‑1 (MQ=255)                                              ccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCg                                                                                                                  >  2:157440/1‑102 (MQ=255)                                                   gTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACg                                                                        >  1:1492332/1‑139 (MQ=255)                                                      tGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatta                                                                                                >  2:707991/1‑112 (MQ=255)                                                      tGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatta                                                                                                <  1:707991/112‑1 (MQ=255)                                                               tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                            <  1:966607/139‑1 (MQ=255)                                                                  caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaaa                                                         >  2:932670/1‑139 (MQ=255)                                                                   aTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCTACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAt                                                        >  2:510572/1‑139 (MQ=255)                                                                    tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                       <  2:1518301/139‑1 (MQ=255)                                                                    tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                       <  2:1825013/139‑1 (MQ=255)                                                                    tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                       <  2:896058/139‑1 (MQ=255)                                                                    tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                       <  1:110454/139‑1 (MQ=255)                                                                      ggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCTACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                     <  1:510572/139‑1 (MQ=255)                                                                                 tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                 >  2:418449/1‑100 (MQ=255)                                                                                 tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                 <  1:418449/100‑1 (MQ=255)                                                                                   aTAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                        <  2:293512/139‑1 (MQ=255)                                                                                   aTAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                        <  1:932670/139‑1 (MQ=255)                                                                                             tGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                 >  1:1503144/1‑72 (MQ=255)                                                                                             tGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                 <  2:1503144/72‑1 (MQ=255)                                                                                                  tCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgca                         >  1:874761/1‑139 (MQ=255)                                                                                                           ccTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGa                                                                                                                 >  1:763584/1‑42 (MQ=255)                                                                                                           ccTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGa                                                                                                                 <  2:763584/42‑1 (MQ=255)                                                                                                                tataATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAAt           >  2:1283471/1‑139 (MQ=255)                                                                                                                    aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGc                                                                                                      >  2:376937/1‑44 (MQ=255)                                                                                                                    aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGc                                                                                                      <  1:376937/44‑1 (MQ=255)                                                                                                                    aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                            <  2:961666/86‑1 (MQ=255)                                                                                                                    aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                            >  1:961666/1‑86 (MQ=255)                                                                                                                    aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCa       >  1:506829/1‑139 (MQ=255)                                                                                                                    aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCa       <  1:1283471/139‑1 (MQ=255)                                                                                                                    aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCa       <  2:874761/139‑1 (MQ=255)                                                                                                                     tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                            >  2:963520/1‑85 (MQ=255)                                                                                                                     tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                            <  1:963520/85‑1 (MQ=255)                                                                                                                     tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaaa                                                         <  2:917083/88‑1 (MQ=255)                                                                                                                     tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaaa                                                         >  1:917083/1‑88 (MQ=255)                                                                                                                     tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAg      >  1:224940/1‑139 (MQ=255)                                                                                                                       aaCACGCTATGTAATTTGTTCTGGCCCCGACAtct                                                                                                            >  1:1136984/1‑35 (MQ=255)                                                                                                                       aaCACGCTATGTAATTTGTTCTGGCCCCGACAtct                                                                                                            <  2:1136984/35‑1 (MQ=255)                                                                                                                       aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCtt                                                                            <  2:716930/67‑1 (MQ=255)                                                                                                                       aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCtt                                                                            >  1:716930/1‑67 (MQ=255)                                                                                                                       aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccg                      <  1:1499272/121‑1 (MQ=255)                                                                                                                       aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccg                      >  2:1499272/1‑121 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCgggg                                                                                                       <  2:1600655/38‑1 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCgggg                                                                                                       >  1:1600655/1‑38 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                 >  1:977230/1‑44 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                                 <  2:977230/44‑1 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                 >  2:407487/1‑60 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                 <  1:407487/60‑1 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                               <  2:556685/78‑1 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCaa                                                               >  1:556685/1‑78 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGc                                         <  2:251624/100‑1 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGc                                         >  1:251624/1‑100 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc  >  2:1300763/1‑139 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc  >  2:1250237/1‑139 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc  >  1:1003924/1‑139 (MQ=255)                                                                                                                         cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGAccc  >  1:624838/1‑139 (MQ=255)                                                                                                                            gCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGc                                                                       >  1:923272/1‑67 (MQ=255)                                                                                                                            gCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGc                                                                       <  2:923272/67‑1 (MQ=255)                                                                                                                                    |                                                                                                                                ACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCC  >  NC_000913/1620076‑1620334

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