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

TGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCA  >  NC_000913/1620084‑1620329                                                                                                                            |                                                                                                                           tgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCg                                                                                                             >  1:1162132/1‑139 (MQ=255) tgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCg                                                                                                             >  2:868507/1‑139 (MQ=255)    cAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACa                                                                                                          <  2:1162132/139‑1 (MQ=255)        gTCGTGGTTGTTCCCGCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACAtctc                                                                                                      <  2:488600/139‑1 (MQ=255)          cGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                    <  1:868507/139‑1 (MQ=255)            tGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgggg                                                                                                  >  2:422582/1‑139 (MQ=255)             ggTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGc                                                                                                 <  1:422582/139‑1 (MQ=255)              gttgttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                    <  2:171810/135‑1 (MQ=255)              gttgttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                    >  1:171810/1‑135 (MQ=255)                  ttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                            >  2:854458/1‑139 (MQ=255)                           cGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCa                                                                                   >  2:1170033/1‑139 (MQ=255)                           cGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCa                                                                                   >  2:903365/1‑139 (MQ=255)                                gAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCttt                                                                              >  1:1300621/1‑139 (MQ=255)                                  acacCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                            <  2:61667/139‑1 (MQ=255)                                     cccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGtt                                                                                                                      <  1:1298278/94‑1 (MQ=255)                                     cccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGtt                                                                                                                      >  2:1298278/1‑94 (MQ=255)                                        accGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCttt                                                                      >  2:459563/1‑139 (MQ=255)                                                     aaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAg                                                         <  2:435559/139‑1 (MQ=255)                                                            tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                  <  1:459563/139‑1 (MQ=255)                                                                         tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGccc                                                                                                               >  1:151541/1‑65 (MQ=255)                                                                         tGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGccc                                                                                                               <  2:151541/65‑1 (MQ=255)                                                                           aTAAGTCGGATCCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                   <  1:854458/139‑1 (MQ=255)                                                                               gTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCgg                                                                                                    >  1:727849/1‑70 (MQ=255)                                                                               gTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTCCTGGCCCCGACATCTCgg                                                                                                    <  2:727849/70‑1 (MQ=255)                                                                                tCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                  <  1:458395/119‑1 (MQ=255)                                                                                tCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                  >  2:458395/1‑119 (MQ=255)                                                                                                        tataATCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAAt      >  1:329742/1‑139 (MQ=255)                                                                                                            aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                <  1:988573/93‑1 (MQ=255)                                                                                                            aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                >  2:988573/1‑93 (MQ=255)                                                                                                            aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCa  <  1:1170033/139‑1 (MQ=255)                                                                                                            aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCa  >  1:617713/1‑139 (MQ=255)                                                                                                            aTCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCa  <  1:903365/139‑1 (MQ=255)                                                                                                             tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                   >  2:978057/1‑105 (MQ=255)                                                                                                             tCAACACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                   <  1:978057/105‑1 (MQ=255)                                                                                                               aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCCCCTTTAc                                                                            <  2:71031/62‑1 (MQ=255)                                                                                                               aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTa                                                                             >  1:478130/1‑61 (MQ=255)                                                                                                               aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTa                                                                             <  2:478130/61‑1 (MQ=255)                                                                                                               aaCACGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                            >  1:71031/1‑62 (MQ=255)                                                                                                                 cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCa                                                           <  2:945374/77‑1 (MQ=255)                                                                                                                 cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCa                                                           >  1:945374/1‑77 (MQ=255)                                                                                                                 cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                >  1:1285776/1‑88 (MQ=255)                                                                                                                 cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACa                                                <  2:1285776/88‑1 (MQ=255)                                                                                                                 cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAg                                     <  2:61375/99‑1 (MQ=255)                                                                                                                 cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAg                                     >  1:61375/1‑99 (MQ=255)                                                                                                                 cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                   >  2:323781/1‑101 (MQ=255)                                                                                                                 cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGcc                                   <  1:323781/101‑1 (MQ=255)                                                                                                                 cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                              <  1:773407/106‑1 (MQ=255)                                                                                                                 cacGCTATGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                              >  2:773407/1‑106 (MQ=255)                                                                                                                    gCTATGTAATTTGTTCTGG‑CCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                  <  2:995563/82‑1 (MQ=255)                                                                                                                    gCTATGTAATTTGTTCTGG‑CCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                                  >  1:995563/1‑82 (MQ=255)                                                                                                                        tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTATGCCACCGCAAGCCAAATa                                                  <  1:188762/79‑1 (MQ=255)                                                                                                                        tGTAATTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTATGCCACCGCAAGCCAAATa                                                  >  2:188762/1‑79 (MQ=255)                                                                                                                            |                                                                                                                           TGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCA  >  NC_000913/1620084‑1620329

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